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Section  1 


EXECUTIVE  SUMMARY  AND  RECOMMENDATIONS 


1.1  BACKGROUND 


The  understanding  and  model  validation  of  electron  beam  (e-beam) 
pumped  XeCl  kinetics  has  lagged  that  of  the  other  excimers,  primarily 
because  of  the  lack  of  appropriate  data  with  which  to  check  the  models. 

The  goal  of  this  program  was  to  measure  fundamental  quantities  in  e-beam 
excited  XeCl  mixtures  and  use  this  data  to  help  improve  and  validate 
existing  codes.  As  an  added  benefit,  some  data  was  also  obtained  in  XeF 
and  KrF  laser  mixtures  in  parallel  programs.  This  data,  whose  collection 
was  not  part  of  the  original  statement  of  work  for  this  program,  has  proven 
very  valuable  in  elucidating  the  differences  between  XeCl  and  the  other 
excimers. 

The  work  performed  during  this  prograjn  resulted  in  a  number  of  papers 
being  written,  some  of  which  have  been  published  in  scientific  journals  and 
others  that  are  still  being  reviewed.  Since  these  papers  describe  the 
experiments,  data,  and  analysis  in  great  detail,  the  bulk  of  this  final 
report  are  reprints  of  these  papers.  They  are  contained  in  the  appendix 
after  this  section.  The  remainder  of  this  section  summarizes  the 
accomplishments  of  this  program  and  ends  with  a  discussion  of  the 
consequences  of  the  work  and  recommendations  for  additional  experiments. 


1.2  PROGRAM  ACCOMPLISHMENTS 

During  this  program  the  following  tasks  were  successfully  completed: 

1)  The  electron  beam  for  the  Tahoma  laser,  used  during  all  the 

measurements,  was  upgraded  to  enable  1  /is  e-beam  pulse  durations  and 
energy  loadings  >100  J/2..  The  long  e-beam  pulse  length  is  important 


i 

MI'!' 


for  measuring  certain  kinetics  phenomenon  such  as  halogen  depletion. 


2)  Time-dependent  electron  densities  iri^XeCl^  {JtelQ  and  KrF  were  measured 
as  a  function  of  initial  halogen  and  minor ^afe 'gas  concentrations. 
The  most  noteworthy  result  of  these  measurements  is  that  the  electron 
density  is  higher  than  was  expected. 


3)  Depletion  of  the  HCl  halogen  donor  in  XeCl  was  measured  after 
termination  of  the  e-beam  pulse.  Although  not  a  time-dependent 
measurement,  the  results  demonstrate  that  a  considerable  amount  of 
halogen  is  depleted  during  the  e-beam  pulse.  This  depletion  is  also 
higher  than  expected. 

4)  Time-dependent  xenon  excited  state  densities  in  XeCl  and  XeF  were 
measured  as  a  function  of  initial  halogen  concentration  and  e-beam 
pulse  duration.  The  xenon  excited  state  density  behaves  in  several 
unexpected  ways  that  are  described  in  the  papers. 

5)  The  XeCl  computer  model  was  upgraded  and  compared  with  the  experimental 
data.  In  particular,  it  was  determined  that  careful  calculation  of  the 
electron  energy  distribution  is  necessary  in  order  to  predict  the 
observed  electron  density  behavior.  The  code  comes  close  to  predicting 
the  measured  electron  density,  xenon  excited  state  density,  and  HCl 
depletion;  however,  best  agreement  between  the  code  and  the  data  occurs 
only  when  additional  HCl  depletion  is  included  above  what  the  model 
predicts . 


As  a  result  of  this  work,  the  following  papers  were  written  and  are 
located  in  the  appendix  of  this  report: 

1.  W.D.  Kimura,  D.R.  Guyer,  S.E.  Moody,  J.F.  Seamans,  and  D.H.  Ford, 
"Electron  Density  Measurements  in  Electron-Beam  Pumped  XeCl  Laser 
Mixtures,"  Appl.  Phys.  Lett.,  1569  (1986). 


2.  W.D.  Kimura,  D.R.  Guyer,  S.E.  Moody,  J.F.  Seamans,  and  D.B.  Ford, 
"Electron  Density  Measurements  in  Electron-Beam  Pumped  XeF  and  KrF 
Laser  Mixtures,"  Appl.  Phys.  Lett.,  60  (1987). 

3.  D.R.  Guyer,  W.D.  Kimura,  and  S.E.  Moody,  "HCl  Burnup  and  Electron 
Densities  in  E-Beam  Pumped  XeCl,"  in  Proceedings  of  the  International 
Conference  on  Lasers  *86,  R.W.  McMillan,  Editor  (STS  Press,  McLean, 
1987),  pp.  589-593. 

4.  F.  Kannari  and  W.D.  Kimura,  "Kinetics  Issues  of  Electron  Beam  Pumped 
KrF  Lasers,"  in  Proceedings  of  the  International  Conference  on  Lasers 
’86,  R.W.  McMillan,  Editor  (STS  Press,  McLean,  1987),  pp.  446-452. 

5.  F.  Kannari,  W.D.  Kimura,  J.F.  Seamans,  and  D.R.  Guyer,  "Xenon  Excited 
State  Density  Measurements  in  Electron-Beam  Pumped  XeCl  Laser 
Mixtures,"  Appl.  Phys.  Lett.,  1986  (1987). 

6.  F.  Kannari,  W.D.  Kimura,  J.F.  Seamans,  and  D.R.  Guyer,  "Xenon  Excited 
State  Densities  in  Electron-Beam  Pumped  XeCl  and  XeF,"  submitted  to 
Journal  of  Applied  Physics. 

7.  F.  Kannari  and  W.D.  Kimura,  "High-Energy  Electron  Distribution  in 
Electron-Beam  Excited  Ar/Kr  and  Ne/Xe  Mixtures,"  to  be  published  in 
Journal  of  Applied  Physics,  April  15,  1988. 

8.  F.  Kannari  and  W.D.  Kimura,  "Low-Energy  Electron  Distribution  in 
Electron-Beam  Excited  XeCl  Laser  Mixtures,"  submitted  to  Journal  of 
Applied  Physics. 

9.  W.D.  Kimura,  F.  Kannari,  J.F.  Seamans,  and  D.R.  Guyer,  "Diagnostic 
Techniques  for  Studying  Excimer  Laser  Kinetics,"  Excimer  Laser 
Technology  and  Advanced  Scientific  Applications,  R.  Sauerbrey  and 
J.H.  Tillotson,  Editors,  to  be  published  in  Proceedings  of  SPIE. 
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Besides  the  preceding  publications,  the  following  papers  were 

presented  at  various  conferences  during  this  program: 

1.  W,D.  Kimura,  D.R.  Guyer,  S.E.  Moody,  J.F.  Seamans,  and  D.H.  Ford, 
"Electron  Density  Measurements  of  XeCl,  XeF,  and  KrF  Laser  Mixtures," 
39th  Gaseous  Electronics  Conference,  Madison,  WI,  Oct.  7-10,  1986, 

Paper  N-1. 

2.  E.T.  Salesky  and  W.D.  Kimura,  "Model  Comparisons  of  Electron  Density 
Measurements  in  KrF,  XeF,  and  XeCl,"  39th  Gaseous  Electronics 
Conference,  Madison,  WI,  Oct.  7-10,  1986,  Paper  N-2. 

3.  D.R.  Guyer,  W.D.  Kimura,  and  S.E.  Moody,  "HCl  Burnup  and  Electron 
Densities  in  E-Beam  Pumped  XeCl,"  International  Conference  on  Lasers 
’86,  Orlando,  FL,  Nov.  3-7,  1986,  Paper  HF.8. 

4.  F.  Kannari,  S.E.  Moody,  and  W.D.  Kimura,  "Kinetics  Issues  of  Electron 
Beam  Pumped  KrF  Lasers,"  International  Conference  on  Lasers  ’86, 
Orlando,  FL,  Nov.  3-7,  1986,  Paper  WK.7. 

5.  E.T.  Salesky,  W.D.  Kimura,  and  F.  Kannari,  "The  Role  of  Electron 
Dissociative  Attachment  and  Impact  Dissociation  of  HCl  in  XeCl  Lasers," 
Conference  on  Lasers  and  Electro-Optics  (CLEO) ,  Baltimore,  MD,  Apr.  27- 
May  1,  1987,  Paper  WR4. 

t 

6.  D.R.  Guyer,  F.  Kannari,  W.D.  Kimura,  and  S.E.  Moody,  "Implications  of 
Recent  Kinetics  Measurements  in  XeCl,"  Atomic,  Molecular,  and  Optical 
Physics  Meeting  of  the  American  Physical  Society,  Cambridge,  MA,  May 
18-20,  1987,  poster  paper. 

7.  F.  Kannari,  J.F.  Seamans,  D.R.  Guyer,  and  W.D.  Kimura,  "Xenon  Excited 
State  Density  Measurements  in  XeCl,"  International  Conference  on  Lasers 
’87,  Lake  Tahoe,  NV,  Dec.  7-11,  1987,  Paper  HK.3. 
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^  8.  W.D.  Kimura,  F.  Kannari,  J.F.  Seamans,  and  D.R.  Guyer,  "Diagnostic 

Techniques  for  Studying  Excimer  Laser  Kinetics,"  0-E/LASE  *88,  Excimer 

[Laser  Technology  and  Advanced  Scientific  Applications,  Los  Angeles,  CA, 
Jan.  10-15,  1988,  Paper  894-09. 

9.  W.D.  Kimura,  F.  Kannari,  and  J.J.  Ewing,  "Recent  Developments  in 
'  Understanding  E-Beam  Pumped  XeCl  Lasers,"  Conference  on  Lasers  and 

r  Electro -Optics  (CLEO) ,  Anaheim,  CA,  Apr.  25-29,  1988,  Paper  TuHl. 

W  Spectra  Technology  also  participated  in  the  Workshop  on  XeCl  Kinetics 

r  held  December  19,  1985  at  SRI  International. 
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1.3  DISCUSSION  AND  RECOMMENDATIONS 

\  ' 

As  a  result  of  the  measurements  performed  during  this  program  and  the 
refinements  made  to  the  computer  model,  our  understanding  of^ XeCl  kinetics 
has  increased  dramatically.  Qualitatively,  one  might  state  that  our 
understanding  of  XeCl  is  now  comparable  to  the  other  major  excimers. 

Much  of  the  mystery  associated  with  the  XeCl  kinetics  has  been 
lifted.  For  example,  significant  vibrational  excitation  of  HCl  appears  to 
be  occurring  as  indirectly  indicated  by  the  behavior  of  the  electron 
density.  The  exact  character  of  the  electron  energy  distribution  and  its 
affect  on  the  kinetics  is  now  better  understood.  The  distribution  is  not 
entirely  Maxwellian,  as  was  generally  assumed  by  models  before  this  work. 
Instead,  only  the  low  energy  distribution  (the  distribution  below  the 
ionization  energy  of  xenon)  is  Maxwellian;  the  high  energy  distribution  is 
non-Maxwellian.  The  problems  with  assuming  a  Maxwellian  distribution  can 
be  avoided  by  using  a  Boltzmann  solution  for  the  electron  energy 
distribution;  however,  we  have  found  that  the  high  and  low  energy  parts  of 
the  distribution  can  be  solved  separately,  thereby  simplifying  the 
calculation  process. 
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As  with  ajiy  research  investigation,  occasionally  one  finds  new 
mysteries.  In  our  case,  the  data  all  indicate  that  more  HCl  is  being  lost 
than  can  be  predicted  by  the  model.  Put  another  way,  the  model  predicts 
quite  well  the  time-dependent  electron  and  xenon  excited  state  density 
behaviors,  and  the  total  HCL  depletion,  but  only  when  additional  HCl 
destruction  to  neutrals  is  artificially  included  in  the  model.  Depletion 
of  HCl  through  dissociative  attachment  and  impact  dissociation  is  already 
in  the  model;  hence,  this  additional  destruction  channel  corresponds  to  a 
process  still  unknown.  It  is  important  to  understand  this  HCl  loss  channel 
because  scaling  of  these  lasers  to  long  pulse  lengths  directly  depends  on 
the  HCl  burnup  rate. 

In  light  of  this  new  mystery,  one  recommendation  regarding  further 
work  in  this  area  are  experiments  to  determine  how  the  additional  HCl  is 
being  lost.  Time-dependent  measurements  of  the  HCl(v  =0,  1,  and  2) 
concentrations  during  the  e-beam  pulse  are  needed.  Performing  these 
measurements  as  a  function  of  buffer  gas  may  also  provide  clues  regarding 
the  role  of  neon  and  xenon  is  this  process. 

As  another  recommendation,  it  was  identified  during  the  Workshop  on 
XeCl  Kinetics  that  measurements  of  the  electron  density  and  temperature  are 
critically  needed.  The  electron  density  has  now  been  measured;  to  measure 
the  electron  temperature  is  much  more  difficult.  However,  with  the  ability 
to  measure  the  electron  density,  it  may  be  possible  to  determine  the 
electron  temperature  by  measuring  phenomena  that  are  dependent  on  both 
electron  density  and  temperature.  For  example,  free-free  absorption  by 
inverse  bremsstrahlung  is  a  function  of  the  electron  velocity 
distribution.  This  distribution  is  dependent  on  the  electron  density  and 
temperature;  knowing  the  density,  one  can  deduce  the  temperature. 

Ultimately,  the  interest  is  in  improving  the  performance  of  XeCl 
lasers,  in  particular,  scaling  them  to  longer  pulse  lengths.  As  mentioned 
earlier,  control  of  halogen  depletion  is  crucial  for  achieving  long  pulse 
lengths.  Therefore,  as  a  final  recommendation,  either  experiments  and/or 


tests  of  the  model  should  be  performed  to  investigate  methods  for  reducing 
the  total  HCl  depletion. 


Electron  density  measurements  of  electron  beam  pumped  XeCI  laser 
mixtures 


W.  D.  Kimura,  Dean  R.  Guyer,  S.  E.  Moody,  J.  F,  Seamans,  and  D.  H.  Ford 

Spi’clru  Tcchnoloiiv.  Imurporuicd.  Sorlhup  M'uy.  Bdlcvue,  Washington 
(Received  2  September  1^)86;  accepted  for  publication  13  October  1986) 


Time-dependent  electron  density  measurements  of  electron  beam  pumped  XeCl  laser  mixtures 
(nonlasing)  are  performed  using  a  CO-  (  10.6/<m)  quadrature  interferometer.  The  electron 
beam  pulse  length  is  -  5(X)  ns  ( full  width  at  half-maximum)  and  delivers  ?  10  A/cm  '  at  the 
foil.  Electron  densities  are  measured  as  a  function  of  the  halogen  (  HCl)  and  xenon 
concentration.  For  a  98 . 3%Ne/l  .s/y  Xe/O.lb*^  UCl  mixture  at  30(X)  Torr.  the  electron 


density  peaks  near  the  beginning  of  the  pulse  at  7  x  10'^  cm  ',  decreases  to  3  ■  10' 


then  gradually  increases  to  :r6  x  10”  cm  '  at  the  end  of  the  pulse.  As  the  initial  halogen 
concentration  is  reduced,  the  electron  density  increases,  rising  dramatically  near  the  end  of  the 
pulse.  This  increase  in  the  electron  density  later  in  the  pulse  is  accompanied  by  a  decrease  of 
the  XeCI*  sidelight  fluorescence. 


Ill  electron  beam  (e-beam)  pumped  excimer  lasers,  the 
primaiN  driving  meehamsm  is  e-beam  generation  of  ions, 
excited  slates,  and  secondary  electrons  The  secondaries  are 
imiiKed  with  the  excimer  formation  through  many  different 
channels.'  Electrons  are  also  important  with  regard  to 
quenching  I'f  the  excimer  molecule  and  of  excited  states 
le.idmg  to  the  excimer  formation.  Due  to  the  governing  role 
of  electrons  in  the  kinetics  of  excimer  lasers,  there  is  a  great 
detil  of  interest  in  determimng  the  magnitude  and  temporal 
beliav  tor  of  the  electron  deiisitv  in  these  devices. 

1  here  have  been  some  measurements  of  the  electron 
density  in  discharge  excimer  lasers'  '  and  in  e-beam  pumped 
pure  rare  gases.''  These  measurements  cannot  be  easily  ap- 
plietl  to  c  -beam  pumped  rare  gas  h.ilide  mixtures  due  to  the 
different  pumping  meehaiiisins  which  exist  m  discharge  la¬ 
sers  and  the  importance  of  electron  attachment.  Recently, 
the  electron  density  in  I'-beam  pumped  Ar/HCI  and  Ar/F- 
mixt  ures  lias  been  deter  minevi.'  fins  let  ter  presents  the  first 
measurements  of  electron  densities  m  I'-beam  pumped  XeCI 


ras  mixtures. 


.■\s  shown  sehemalieallv  in  Fie.  1.  the  electron  densitv  is 


measured  using  a  CO-  (  10  6//m  )  quadrature  interferometer 
in  which  the  c-beam  pumped  gas  chamber  is  located  in  one 
arm  of  the  interferometer  This  gas  chamber  is  part  of  the 
Tahoma  laser  system.'  fhe  laser  mirrors  have  been  remov  ed 
and  the  fused  silica  windows  on  the  chamber  replaced  with 
BaF  ,  w  indow  s.  The  CO  -  laser  beam  has  a  diameter  of  ~  I  cm 
within  the  laser  chamber.  As  is  shown  in  Fig.  1,  the  beam 
passes  twice  through  the  chamber  approximately  3  cm  from 
the  foil  and  centered  around  the  nominal  lasing  axis.  Liquid- 
nitrogen-cooled  FlgCdTe  detectors,  with  1.5  mm  -  1.5  mm 
active  areas,  monitor  the  transmitted  energy  in  the  probe 
and  reference  beams. 

Quadrature  detection  offers  several  advantages  over 
conventional  interferometry.  '  The  primary  advantage  is  the 
elimination  of  ambiguous  interpretation  of  the  detector  sig¬ 
nals  when  either  passes  through  an  extremum.  This  is  be¬ 
cause  when  one  detector  is  at  .in  extremum,  the  other  detec¬ 
tor  will  be  at  its  maximum  differential  sensitivity  and 
provieies  enough  information  to  determine  the  sign  of  the 
phase  shift.  The  output  from  the  detectors  is  recorded  on  a 
dual-beam  oscilloscope  for  later  analysis. 
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The  index  of  refraction  of  a  weakly  ionized  plasma  is 
given  by 


(  1  -  )'  ~\ 


where  is  the  pl.isma  fretjuenev  and  <l>  is  the  optical  fre¬ 
quency.  Since  i;):,  -  14^1  l/i,  -  .VIX  ■  10’/;  ,  where  c 
IS  the  electron  charge,  /n  is  the  electron  mass,  and  is  the 
electron  density  in  cm  ',  the  index  change  due  to  electron 
density  is  A/;,.  -  .‘'.0,^  ■  10  "/;  The  I'-beam  dimensions 

are  10  cm  >  }5  cm  (the  c  beam  veas  baffled  down  to  half  its 
nominal  length  of  70  cm  );  hence,  the  total  optical  path  of  the 
laser  plasma  is  2  ■  .T'i  70  cm  t  herefore,  a  shift  of  one 
fringe  on  our  system  corresponds  to  an  electron  density  of 
n,  =  .VOX  lO'-*  cm  ‘ 

Infrared  interferometry  to  measure  electron  densities  is 
a  well  established  technique;  however,  special  care  was  taken 
to  eliminate  possible  errors  in  the  measurement  The  detec¬ 
tors  are  placed  a  considerable  distance  (  ■  1  K  m )  from  the 
laser  chamber  to  eliminate  detection  of  infrared  emission 
from  the  XeCl  mixture.  ‘  Tare  shots  w  ere  taken  w  ith  the  ref¬ 
erence  arm  of  the  interferometer  blocked  to  ensure  that  no 
stray  emission  was  being  detected,  and  that  no  loss  of  signal 
was  occurring  due  to  absorption  or  transient  refractive  index 
effects  '  "  The  issue  of  refractive  index  changes  causeii  by 
resonances  occurring  nearby  the  probe  wavelength  has  been 
considered  by  others’  and  determined  not  to  be  .i  significant 
problem  This  agrees  with  mir  own  estimaiion  that  fringe 
shifts  caused  by  nearby  transitions,  such  as  the  '>  and  I  I  ' 
//m  xenon  lines,"  contribute,  at  most,  a  fraction  of  a  fringe 
shift  In  addition,  some  data  were  taken  at  IDS  f/nt  and 
yielded  essentially  identical  results  as  the  10  0//m  data  Re¬ 
actions  occurring  within  the  mixture  can  .iKo  induce  an  in¬ 
dex  change,  hut  even  for  "worse  case"  estimates,  they  wouhl 
contribute  only  a  fraction  of  a  fringe  shift  The  ret  urn  to  zero 
of  the  detector  signal  within  a  few  microsecoiuls  .ilier  termi- 
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luition  of  the  pump  pulse  indicates  that  acoustic  waves  in  the 
gas  chamber  do  luil  contribute  to  the  detected  index  change. 

The  detector  .^-dB  bandwidths  are  s  W)  VIHz  This 
causes  some  loss  of  signal  amplitude  during  periods  of  rapid 
electron  density  changes;  however,  it  does  not  appreciably 
affect  the  accuracy  of  the  results  The  largest  uncertainty  m 
the  data  is  caused  by  the  lower  bandw  idth  limit  of  the  detec¬ 
tor  amplifiers  (5  MHz)  We  estimate  the  uncertainty  to  he 
approximately  *  0,4  fringes  (  1 .2  •  10'‘cm  ') 

In  all  the  gas  mixtures,  the  halogen  donor  is  HCI  and  the 
diluent  IS  neon.  Total  gas  pressure  is  ftXK)  I  orr  (  .i  '  ama- 
gats)  and  the  initial  temperature  is  2‘>4  k  I'owerifepositum 
in  the  gas  IS  approximately  ISOkW/cm'l  10.\  cm  atthe 
foil ).  and  IS  measured  using  a  pressure  rise  technique  '  The 
e-beam  source  uses  a  carbon  felt,  cold  cathode  at  .VK)  kV  and 
::  2.'  k.\  peak  Current 

figure  2  shows  typical  time-dependent  electron  density 
results  for  varunis  initial  HCI  concentrations  (the  .Xe  con¬ 
centration  IS  consiaiu  at  I  .‘'‘"r  )  A  normalized  ty  pical  e- 
beam  curreni  p.iKe  shape  is  draw  n  at  the  bottom  of  the  fig¬ 
ure  During  the  tirst  50  ns  of  the  e-beam  pulse,  the  electron 
density  rises  sh.irplv,  then  gradually  decreases,  with  the 
amount  of  decrease  dependent  on  the  initial  halogen  concen¬ 
tration  Since  the  dissociative  attachment  rate  constants  for 
HCI  increase  dramatically  with  vibrational  excitation.  "  " 
this  behavior  can  he  attributed  to  the  proifuction  of  HCI 
(  c  -  0)  during  the  pumping  process  for  a  typical  laser  mix¬ 
ture  of  08  .  ,VxNe/l  .Xe/0  lb*';  HCI,  the  electrr'n  den¬ 
sity  vanes  between  and  7  ■  10"  cm 

A  surprising  result  is  the  magnitude  of  the  electron  den¬ 
sity  It  IS  roughly  4  times  larger  than  that  predicted  by  a 
computer  moilel  developeil  at  Spectra  lechnology.  Inc  ' 

I  his  h.is  .1  n umbel  of  important  implic.it ions,  some  c'f  w  hich 
■ire  discussetl  elsewhere’'  "  .md  will  luU  be  examined  here 

Wlien  lasmg  ,\eCl  on  the  Tahoina  device,  lumimal  HCI 
concenlrations  range  from  0 to  0  Ib'  f  f  rom  the  elec¬ 
tron  ilensiiy  results,  it  is  clear  why  leaner  HCI  concentra- 
tionsdo  not  perform  as  well  .-\l  le.m  HCI  concentrations,  the 
electron  ilensiiy  reaches  very  high  ifensities  i  -  lO"  cm  ’). 
which  implies  loss  of  HCI  fuel  .md  mcreaseil  electrr'n 
quenching  1  his  effect  can  be  seen  more  rle.irly  m  fig  .1 
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F'lG  4  Tspical  cleclnm  densities  of  XeCl  mixtures  (OOSG-  HCT.  balance 
Ne  lo  3(XX)  T<'rr )  for  varuius  initial  Xe  concentrations.  Data  points  are  not 
shoNMi,  hov^e^er.  a  representative  error  bar  is  indicated 


which  shows  the  XeCl*  sidelight  fluorescence  for  the  same 
mixtures  given  in  Fig.  2.  Although  the  electron  density  re¬ 
mains  "under  control"  in  the  rich  HCl  mixtures  (>0.16% 
HCl ).  Fig.  3  demonstrates  that  the  overall  XeCl*  population 
decreases  with  HCl  concentrations  above  ~0.08%. 

Figure  4  shows  the  electron  density  for  various  xenon 
concentrations,  with  the  HCl  concentration  fixed  at  0.08%. 
There  is  little  change  with  xenon  concentration  over  the 
range  of  0. 5-3.0%. 
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Electron  density  measurements  of  electron-beam-pumped  XeF  and  KrF  laser 
mixtures 
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A  CO.  quadrature  inlerferonieier  is  used  to  measure  the  temporal  evolution  of  the  electron 
density  in  XeF  and  KrF  electron-beam-)  umped  laser  mixtures  (nonlasing).  Measurements  are 
obtained  using  a  dOO-kV  electron  beam  which  delivers  S  10  A/cm’  at  the  foil  for  pulse 
durations  ^  650  ns  (FWHM).  For  typical  XeF  mixtures  (99.35%  Ne/0.5%  Xe/0. 15%  F,  at 
2888  Torr)  the  electron  density  varies  between  3  and  4.5  x  10'^  cm  ’  during  the  pulse.  The 
electron  density  is  found  to  increase  as  the  initial  F,  concentration  is  decreased.  Both  10%  Kr 
in  argon  diluent  and  argon-free  KrF  mixtures  are  tested.  The  electron  density  for  these 
mixtures  has  the  same  magnitude  and  behavior  as  the  typical  XeF  mixture. 


Several  computer  models  of  electron-beam  (e-beam )- 
pumped  XeF  and  KrF  lasers  exist.  All  the  models  rely  on 
verification  using  experimental  data  on  the  laser  perfor¬ 
mance,  for  example,  the  laser  output  flux  and  efficiency.  The 
models  have  been  very  useful  in  improving  the  understand¬ 
ing  of  the  complex  plasma  chemistry  which  occurs  in  these 
lasers.  One  key  area  is  the  role  of  the  primary  and  secondary 
electrons  which  are  very  important  in  the  formation  of  the 
excimers,  both  through  ionic  and  metastable  channels^  and 
in  quenching  of  excited  states.'’  Despite  the  importance  of 
the  electrons,  none  of  the  codes  has  been  validated  with  re¬ 
gard  to  their  predictions  of  the  electron  densities,  primarily 
because  of  the  lack  of  experimental  data. 

This  letter  reports  the  first  measurements  of  electron 
densities  in  e-beam-pumped  XeF  and  KrF  ( 10%  Kr  and 
argon-free)  laser  mixtures  under  nonlasing  conditions.  Elec¬ 
tron  density  measurements  have  been  previously  reported  in 
discharge  pumped  XeCr  "  and  XeF^ ;  and,  recently,  the  elec¬ 
tron  density  in  e-beam-pumped  Ar/F,  and  Ar/HCI  has  been 
measured."’  "  An  earlier  letter’’  gives  the  electron  density 
results  for  e-beam-pumped  XeCl  laser  mixtures. 

The  experimental  system  is  described  in  more  detail 
elsewhere.’’  Briefly,  the  gas  chamber  of  the  Tahoma  laser 
system ’’is  placed  in  one  arm  ofaCO,  ( lO.b/rm)  quadrature 
interferometer.’''  The  electron  density  is  determined  from 
the  refractive  index  change  of  the  plasma  due  to  the  presence 
of  free  electrons.  In  this  particular  experiment,  one  fringe 
shift  corresponds  to  an  electron  density  of  3.0 X  lO’"*  cm  '. 

The  electron  beam  source  is  a  carbon  felt,  cold  cathode 
at  300  kV  with  a  peak  emission  current  of  ~23  kA.  The 
current  density  at  the  foil  (0.13  mm  aluminized  Kapton)  is 
^  10  A/cm'.  A  pressure  rise  technique’'  is  used  to  calibrate 
the  power  deposition  in  the  gas  mixtures.  In  the  XeF  gas 
mixtures,  the  deposition  is  —180  kW/cm'.  In  the  KrF  gas 
mixtures,  the  deposition  is  kept  constant  at  120  kW/cm' 
by  adjusting  the  total  pressure  of  the  gas  mixtures  to  com¬ 
pensate  for  the  different  electron  stopping  powers  of  Ar  and 
Kr.  All  measurements  are  performed  at  294  K. 

Figure  1  shows  typical  time-dependent  electron  density 
results  for  XeF  as  a  function  of  the  F,  concentration.  The 
XeF  mixture  is  0.5%  Xe  in  Ne  diluent  to  a  total  gas  pressure 
of  2888  Torr.  For  reference.  Fig.  2  shows  a  normalized  typi¬ 


cal  trace  of  the  excimer  sidelight  fluorescence  for  the  0. 1 5% 
F,  XeF  mixture  ( the  same  trace  is  also  representative  of  the 
KrF*  sidelight  fluorescence  for  the  data  given  in  Fig.  3 ) .  The 
top  of  Fig.  3  gives  a  normalized  typical  pulse  shapie  of  the  e- 
beam  current.  The  e-beam  pulse  length  does  vary  slightly 
shot  to  shot;  some  of  the  data  curves  (e  g.,  the  0.15%  and 
0.20%  Fj  mixtures  in  Fig.  1 )  represent  data  for  pulse  lengths 
shorter  than  shown  in  Fig.  2  and  the  top  of  Fig.  3. 

As  seen  in  Fig.  1,  as  the  initial  halogen  concentration  is 
reduced,  the  electron  density  tends  to  increase  more  rapidly 
towards  the  end  of  the  pulse.  This  is  not  surprising,  since  as 
the  F,  is  depleted  (burned  up),  due  to  the  dissociative  at¬ 
tachment  with  electrons,  less  F,  is  available  to  control  the 
electron  density.  Without  the  attacher,  the  electron  density 
appears  to  approach  the  higher  densities  found  in  e-beam- 
pumped  pure  rare  gas  mixtures. Dissociative  recombina¬ 
tion  of  the  electrons  with  dimer  ions  does  not  help  control 
the  electron  density  because  the  dimer  ion  number  density 
tends  to  be  too  small’  ’  for  the  conditions  of  this  experiment. 
The  rapid  increase  in  the  electron  density  as  the  halogen  is 
depleted  has  been  predicted  in  e-beam  sustained  discharge 
models,''  but  the  increase  tends  to  be  at  least  an  order  of 

20 

Ne/0  5%  Xe/Fp 


FIG- 1  Typical  electron  densities  of  XeF  mixtures  (0  Xe.  balance  Ne  to 
2888  Torr  I  for  various  initial  F  concentrations  Data  points  are  not  shown; 
however,  a  representative  error  bar  is  indicated 
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FIG  2.  Normalized  waveform  of  typical  XeF*  and  KrI  *  sidelight  fluores¬ 
cence  emission  for  the  data  given  in  Figs  I  and  .V 

magnitude  larger  than  the  results  in  Fig.  1  Since  the  mecha¬ 
nism  for  electron  population  growth  m  discharge  lasers 
(metastable  ionization  leading  to  electron  avalanche  and 
arcing)  is  very  different  from  that  in  e-beam-pumped  sys¬ 
tems,  care  must  be  taken  when  comparing  the  two  pumping 
schemes. 

From  past  experience'^  on  the  Tahoma  laser  device,  for 
these  pulse  lengths,  the  optimum  laser  performance  on  XeF 
occurs  at  Fi  concentrations  of  0.10-0. 159(-.  At  higher  F, 
concentrations,  although  the  electron  density  may  decrease, 
there  is  a  significant  increase  in  F,  quenching  and  F,  absorp¬ 
tion  which  reduces  the  laser  efficiency.  The  electron  density 
profiles  in  Fig.  1  do  not  display  a  peak  in  the  density  at  the 
beginning  of  the  pulse  as  is  seen  in  the  XeCl  data.  '• 

The  electron  density  results  for  the  KrF  mixtures  are 
given  in  Fig.  3.  The  mixtures  are  8^.7%  Ar/10%  Kr/0.27% 
F,  at  1007  Torr,  and  99.6%  Kr/0.4%  F,  at  665  Torn  For 
both  mixtures,  the  F,  partial  pressure  is  2.7  Torr,  corre¬ 
sponding  to  a  density  of  8.9  x  10'"  cm  '.  The  electron  den¬ 
sity  behavior  is  very  similar  between  the  two  KrF  mixtures. 
The  slight  increase  in  density  towards  the  end  of  the  pulse  is 
apparently  due  to  F_,  burn-up. 

To  investigate  the  effects  of  high-energy  loading,  mea¬ 
surements  are  performed  on  the  10%  Kr  and  Kr/F,  mix¬ 
tures  for  e-beam  pulse  lengths  of  650  ns  ( FWHM ).  The  re¬ 
sults  are  depicted  in  Fig.  4.  The  F_,  concentration  has  been 
increased  in  the  10%  Kr  mixture  to  5  Torr  (1.6x10'^ 
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FIG  3.  Typical  electron  densities  of  KrF  mixtures  ( 8*)  7^  \x/\Q*7f  Kr/ 
()  27^r  F.,al  1007  Torr.  and  99  Kr/0  4^£  F  -at  66?  Torr )  The  top  trace 
IS  the  normalized  e-beam  current  pulse.  Data  pt)ints  are  not  shown;  how¬ 
ever.  a  representative  error  bar  is  indicated 
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FIG,  4  Typical  electron  densities  of  KrF  mixtures  (  89, 5G  Ar/lO'^r  Kr/ 
()-5n;-  1007  Torr.  and  9n4^r  Kr/0  6G-  F,  at  665  Torr ),  The  top  trace  is 

the  normalized  e-beam  current  pulse.  Data  pHunls  are  not  shown;  however.  ^ 
a  representative  error  har  is  indicated 


cm  '),  and  in  the  Kr/F-  mixture  to  4  Torr  (1.3  a  10'’ 
cm  '),  At  this  higher  F,  concentration,  the  electron  density 
for  the  Ar/Kr/F-  mixture  does  not  increase  after  300  ns  as 
seen  in  Fig.  3,  but  remains  essentially  constant.  On  the  other 
hand,  the  Kr/F-  mixture  still  displays  a  gradual  increase  in 
the  electron  density  towards  the  end  of  the  pulse.  Although 
the  initial  F-  concentration  in  the  Kr/F .  mixture  is  slightly 
less  than  in  the  Ar/Kr/F_,  mixture,  this  tendency  for  the 
electron  density  to  be  higher  in  the  Kr/F,  mixtures  is  consis¬ 
tent  with  model  predictions’  and  long  pulse  extraction  mea¬ 
surements.'"  Kr  ‘  is  formed  more  efficiently  in  Kr/F,-  than 
Ar*  is  in  Ar/Kr/F-,''  which  means  more  electrons  are 
freed  per  unit  of  power  deposition.  Higher  gas  healing  effects 
in  the  Kr/F-  mixtures  also  slightly  increase  the  dissociative 
attachment  rate  constant  of  F-, '''  which  leads  to  a  faster  de¬ 
pletion  of  the  halogen  fuel  and  an  increase  in  the  electron 
density. 

The  KrF  mixtures  given  in  Figs.  3  and  4  represent  opti¬ 
mum  mixtures  on  the  Tahoma  laser.  For  the  optimum  XeF 
and  KrF  laser  mixtures,  the  electron  density  varies  between 
3  and  4.5  x  10'^  cm  *.  The  density  can  be  estimated  analyti¬ 
cally  from  the  expression 

Py(eW,) 

n,  = - ,  1 ) 

A„[F-1 

where  is  the  electron  density,  P,  is  the  input  power  deposi¬ 
tion.  (efV, )  is  the  energy  loss  per  ion  pair  produced,’"  is 
the  F;  attachment  rate  constant,  and  [  F,  ]  is  the  F,  number 
density.  Using  a  typical  F-  attachment  rate  constant  of 
2x  10  “'em  Vs,'"  and  values  for  (elF', )  of  36  and  26  eV  for 
Ne  and  Ar,  respectively,’"  we  find  that  electron  densities  of 
1-2  X  lO'"  cm  '  are  predicted  which  are  2-3  times  smaller 
than  measured. 

Although  Eq.  ( 1 )  may  predict  lower  electron  densities, 
preliminary  comparisons’'  with  a  computer  model’  yield 
good  agreement  in  overall  magnitude  and  temporal  behavior 
with  the  data.  The  model  uses  a  time-dependent  Boltzmann 
equation  solution”  and  indicates  that  the  secondary  electron 
energy  distribution  is  significantly  non-Maxwellian  for  the 
pump  conditions  of  this  experiment.  It  also  predicts  that  the 
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F,  dissociative  attachment  rate  constant  is  lower  than  that 
assumed  in  Eq.  ( I ),  which  leads  to  higher  electron  densities. 

The  reasons  the  electron  distribution  is  non-Maxwellian 
are  complex  because  the  distribution  is  affected  by  a  variety 
of  interrelated  kinetic  processes.  For  example,  the  electron 
dissociative  attachment  rate  constant  of  F,  is  directly  depen¬ 
dent  on  gas  temperature'''  and  inversely  dependent  on  elec¬ 
tron  temperature.'''-'  Flowever,  Rosenberg  cr  o/. show 
that  the  dissociative  attachment  rate  constant  is  also  a  func¬ 
tion  of  the  F.  concentration  and  c-beam  current  density,  and 
that  adding  F,  has  the  net  effect  of  increasing  the  average 
electron  temperature,  caused  by  removal  of  low-energy  elec¬ 
trons  through  attachment.  This  affects  the  electron  energy 
distribution  and,  at  the  concentrations  and  densities  used  in 
this  work,  may  lead  to  a  lower  attachment  rate  constant. 
Other  implications  of  the  electron  density  results  are  report¬ 
ed  elsewhere.-'* 
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Abstract 


Two  separate  topics  are  discussed  in  this  paper.  The  first  examines  the  kinetic 
processes  associated  with  the  secondary  electrons  in  electron-beam  excited  rare-gas  halide 
lasers.  The  results  of  experimental  measurements  of  the  secondary  electron  density  in  KrF 
laser  mixtures  suggest  that  the  assumption  of  a  Maxwellian  electron  energy  distribution, 
used  in  some  computer  models,  may  not  be  valid  under  certain  circumstances.  Calculations 
of  the  high  energy  distribution  of  the  secondary  electrons  produced  by  the  electron  beam 
also  predict  a  non-Maxwell ian  distribution.  Moreover,  the  calculated  W-values 
[eV/electron-ion  pair]  and  yields  of  rare  gas  metastables  show  a  dependence  on  the  laser 
mixture  composition.  The  second  topic  of  this  paper  deals  with  time-resolved  calculations 
of  multiplexed  short-pulse  amplification  in  KrF.  The  calculations  include  nonsaturable  and 
saturable  absorption,  amplified  spontaneous  emission  effects,  and  end  mirrors  for  a  double 
pass  amplification  arrangement.  It  is  found  that  a  train  of  short  pulses,  with  a  short 
time  interval  between  the  pulses,  can  extract  energy  from  the  gain  media  with  little  cost 
to  extraction  efficiency  and  is  comparable  to  the  performance  of  a  contiguous  train  of 
short  pulses.  This  suggests  that  a  novel,  simpler  angular-multiplexing  system  can  be 
designed  for  high  energy,  short  pulse  KrF  laser  systems. 


Introduction 


Efforts  are  currently  in  progress  to  scale  up  electron-beam  (e-beam)  pumped  KrF  laser  to 


multi-10  kJ  levels,  to  serve  as  drivers  in  advanced  inertial  confinement  fusion.  Kinetics 
models  of  the  KrF  laser  are  very  importajit  tools  for  predicting  the  output  performances  of 
such  systems.  However,  there  are  several  important  kinetics  issues,  which  affect  the  basic 
design  of  these  models,  that  are  still  not  well  understood.  Some  of  these  issues  are:  1) 
The  secondary  electron  energy  distribution;  2)  the  yields  of  primary  ionization  and 
excitation;  3)  vibrational  relaxation  in  KrF*(B)  and  (C)  states;  4)  mixing  between  KrF  (B) 
and  (C)  ;  5)  Kr2F  formation  channels  in  atmospheric  pressure  laser  mixtures;  6) 
nonsaturable  absorptions;  7)  the  temperature  dependences  of  the  absorbers;  and  8) 
photoassociative  KrF  formation.  Further  development  of  the  models  requires  additional 
fundamental  kinetics  measurements,  such  as  measurements  of  the  secondary  electron  density 
and  the  electron  mean  temperature,  the  fluorescence  yield  of  the  KrF*(C)  state, 
time-resolved  saturated  absorption  at  various  gas  temperatures,  and  so  on. 


One  of  these  measurements,  the  secondary  electron  density,  has  been  measured 
recently.  The  results  for  e-beam  pumped  KrF  laser  mixtures^  indicate  that  the  density 

is  significantly  higher  than  that  calculated  by  some  models.  This  seems  to  imply  that  the 
treatment  of  the  secondary  electron  processes,  which  generally  has  assumed  a  Maxwellian 
distribution  for  the  electron  energy  distribution,  is  invalid.*^  In  light  of  these 
measurements,  the  first  half  of  this  paper  examines  the  importance  of  the  secondary 
electron  processes  with  regard  to  the  kinetics  models.  Presented  are  the  calculated 
results  of  the  e-beam  source  energy  distribution  in  the  gas,  and  new  results  pertaining  to 
the  W-value  (the  energy  expended  per  ion  produced)  and  metastable  yields  in  the  gas 
mixtures . 


The  second  half  of  this  paper  deals  with  another  important  topic  related  to  KrF  laser 
fusion  drivers,  and  that  is  multiplexed  short  pulse  amplification.  Recent  measurements 
have  shown  that  a  train  of  short  pulses,  with  an  appropriate  time  interval  between  the 


pulses,  caul  extract  the  same  energy  as  that  by  long  pulses  consisting  of  a  train  of 
contiguous  short  pulses  without  intervening  time  spaces.  This  is  because  of  repumping 
occurring  during  the  time  spaces.  This  dynamic  characteristic  of  KrF*  permits  the  design 
of  a  novel,  simpler  multiplexed  short  pulse  amplifier  system.  By  reducing  the  number  of 
multiplexed  beams  needed,  while  maintaining  a  similar  output  energy  as  the  concatenated 
multiplexed  system,  the  complexity  and  cost  of  the  system  is  reduced.  In  the  latter  half 
of  this  paper,  the  results  of  a  computer  simulation  of  this  amplification  concept  are 
described.  Examined  are  extraction  efficiency  and  amplified  spontaneous  emission  (ASE) 
effects  as  a  function  of  the  time  space  between  the  pulses. 

Electron  Density  Measurements  and  E-Beam  Source  Calculation 


A  CO2  quadrature  interferometer  is  used  to  measure  the  temporal  evolution  of  the 
electron  density  in  e-beam  pumped  KrF  laser  mixtures.  The  experimental  system  has  been 
described  extensively  elsewhere . ^ ^  Figure  1  shows  the  results  for  mixtures  of 
89. 5X  Ar/lOX  Kr/0.5X  F2  at  1007  torr,  and  99. 4X  Kr/0.6X  F2  at  665  torr .  The  pumping  pulse 
width  of  the  electron  beam  is  «700  ns  with  an  average  excitation  rate  of  oi20  kW/cm"’. 
Measured  electron  densities  are  ^3  x  10^  cm  for  both  mixtures  at  early  times  during  the 
pumping . 

The  secondary  electron  density  can  be  estimated  by  using  the  excitation  rate,  a  rate 
constant  for  F2  dissociative  attachment,  and  a  W-value .  Assuming  a  secondary  electron 
temperature  of  1  eV  and  a  corresponding  attachment  rate  of  4.5  x  10  cm"^/sec,  which  are 
the  same  assumptions  used  by  Trainer  and  Jacob  in  their  paper  on  the  electron  quenching  of 
KrF*,  an  electron  density  of  5  x  10^  cm  is  calculated.  A  detailed  model  calculation, 
using  a  Maxwellian  assumption  for  the  energy  distribution  of  secondary  electrons  in 
conjunction  with  the  energy  balance  equation,  gives  an  electron  density  of  1  x  10^“^  cm  ^ . 
Hence,  the  magnitude  of  the  measured  electron  density  is  roughly  3  to  6  times  larger  than 
that  predicted  by  the  simple  analysis  and  the  model  calculation. 


E-BEAM  CURRENT 


89  7%  Ar/10%  Kr/O  27%  F, 


Figure  1.  Typical  electron  densities  of 
KrF  mixtures 

(89. 5X  Ar/lOX  Kr/O.SX  Fg  at 
1007  torr,  and 
99. 4X  Kr/0.6X  Fg  at  665 
torr) .  The  top  trace  is  the 
normalized  e-beam  current 
pulse . 
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To  understand  the  reasons  for  these  differences,  it  is  helpful  to  review  the  energy  flow 
channels  for  the  electrons  in  an  e-beam  pumped  laser  mixture.  First,  the  primary  electrons 
from  the  e-beam  directly  produce  ions,  secondary  electrons,  and  excited  molecules.  This 
process  is  usually  accounted  for  in  the  models  by  using  a  stopping  power  constant  for  the 
gases.  The  hot  secondary  electrons  are  cooled  by  collisions  with  gas  molecules  and  produce 
more  ions  and  metastables.  This  process  is  usually  modeled  by  using  W-values  and  the  yield 
of  metastables  relative  to  the  ion  formation.'^  W-value  is  defined  as  the  energy  loss  per 


electron-ion  pair  produced,  and  is  generally  higher  (typically  a  factor  of  2)  than  the 
ionization  potential  of  the  atom. 

The  primary  electrons  from  the  e-beam  auid  the  hot  secondary  electrons  are  generally 
treated  collectively  as  the  e-beam  source.  Because  the  preceding  two  processes  are  very 
fast  (usually  <100  ps  for  typical  e-beam  voltages  and  laser  mixtures) ,  the  formation  of  the 
ions  and  metastables  by  the  e-beam  source  has  been  treated  as  an  nstantaneous  process  by 
simply  using  gross  macroscopic  values  for  the  stopping  power  constant  atnd  W-values . 

The  cold  secondary  electrons  redistribute  their  energy  through  slow  electron  reactions 
such  as  dissociative  attachment  with  the  halogen,  excitation  and  ionization  of  metastables, 
recombination  with  positive  ions,  elastic  collisions  with  neutral  gases,  superelastic 
collisions,  and  electron-electron  collisions.  This  redistribution  can  be  calculated  by 
solving  the  Boltzmann  transport  equation,  as  is  usually  done  for  discharge  pumped  lasers. 
However,  in  e-beam  pumped  models,  the  distribution  of  the  cooled  secondary  electrons  is 
many  times  either  ignored  or  assumed  to  be  Maxwellian. 

To  illustrate  how  non— Maxwellian  the  electron  distributions  can  become  in  e-beam  pumped 
gases.  Fig.  2  shows  the  energy  distribution  of  the  e-beam  source  calculated  for  a  laser 
mixture  of  380  torr  Ar/380  torr^Kr.  The  method  used  for  this  calculation  is  based  on  that 
described  by  Elliott  and  Green.'*  The  calculations  in  this  paper  also  include  excitation  of 
the  ground  state  rare-gases,  as  well  as  their  ionizations. 


Figure  2.  Energy  distribution  of  the 

e-beam  source  in  a  mixture  of 
50%  Ar/50*  Kr  at  760  torr. 


_  1  1  i  1.  .  J 

0  10  20  30  40  50  60 

ENERGY  (eV) 

In  Fig.  2,  the  e-beam  source  distribution  is  truncated  at  a  first  excitation  threshold 
of  Kr(53)  at  10  eV.  Note,  the  distribution  below  this  first  excitation  threshold  is  very 
close  to  being  Maxwellian.  However,  the  distribution  above  this  threshold  is  no  longer 
Maxwellian.  Although  the  slow  electron  reactions  of  the  low  energy  electrons  change  the 
distribution,  especially  in  the  low  energy  part  of  the  distribution  below  the  excitation 
threshold,  the  high  energy  tail  of  the  distribution  remains  relatively  unchanged  as  it  is 
driven  by  the  e-beam  source  because  the  effects  of  slow  electron  reactions  are  negligible 
for  high  energy  electrons. 

Figure  3  shows  how  the  effective  W-value  for  the  Ar  ion  changes  as  a  function  of  the  Kr 
concentration  in  an  Ar/Kr  mixture.  The  effective  W-value  is  the  W-value,  calculated  from 
the  total  energy  deposition,  times  the  fraction  of  energy  deposited  in  Ar .  Also  plotted  in 
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Figure  3.  Plots  of  the  effective 

W-values  for  the  Ar  ion  and 
yields  of  Ax  metastables 
relative  to  the  Ar  ion 
formation  in  Ar/Kr  mixtures  as 
a  function  of  Kr 
concentration. 


Fig.  3  axe  the  Ar  metastable  yields,  both  for  individual  states  and  for  the  total  yield, 
relative  to  the  ion  formation.  It  can  be  seen  that  the  W-value  is  not  a  unique  constant, 
but  depends  on  the  mixture  composition.  The  total  yield  of  the  metastables  is 
approximately  0.45,  which  is  much  higher  than  the  amount  typically  used  in  the  models  of 
0.3.  Figure  4  shows  the  same  calculations  for  Kr  in  an  Ar/Kr  mixture. 

A  Boltzmann  code  calculation  is  necessary  to  calculate  the  low  energy  distribution  of 
the  secondary  electrons  and  to  obtain  accurate  rate  constants  for  the  secondary  electron 
reactions.  This  is  especially  true  with  regard  to  dissociative  attachment  by  fluorine 
molecules,  which  can  affect  the  low  energy  part  of  the  electron  distribution.  A  detailed 
Boltzmaoin  calculation  for  the  low-energy  electron  distribution,  in  conjunction  with  a 
proper  calculation  of  the  e-beam  source  distribution  for  the  high  energy  electrons,  will 
significantly  improve  the  models. 
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Figure  4.  Plots  of  the  effective 

W-values  for  the  Kr  ion  and 
yields  of  Kr  metastables 
relative  to  the  Kr  ion 
formation  in  Ar/Kr  mixtures  as 
a  function  of  Kr. 


Multiplexed  Short  Pulse  Amplification 

This  last  section  of  the  paper  discusses  calculations  related  to  short  pulse 
amplification  in  e-beam  pumped  KrF  amplifiers.  Figure  5  is  a  schematic  comparison  of  the 
transient  gain  dynamics  of  a  cw  amplifier,  which  corresponds  to  the  usual  concatenated 
multiplexed  amplifier,  and  the  gain  dynamics  of  a  multiplexed  short  pulse  amplifier  with 
intervening  time  periods  between  pulses.  With  cw  amplification,  the  gain  is  depleted  and 
kept  constant  during  the  energy  extraction.  Hence,  the  number  of  beams  required  for  this 
scheme  is  simply  calculated  by  dividing  the  pumping  pulse  length  by  the  input  laser  pulse 
length.  The  KrF  laser  medium  is  not  an  energy  storage  one;  the  lifetime  of  KrF*  is 
typically  3  to  5  ns .  Therefore,  if  the  input  laser  pulses  to  the  amplifier  are  arranged  in 
a  series  of  short  pulses  with  an  appropriate  time  interval  between  the  pulses,  the  gain 
medium  is  repumped  during  the  interval,  and  the  following  short  pulse  can  see  the  highly 
recovered  gain  and  extract  the  extra  energy  stored  during  the  interval.  In  this  scheme, 
the  number  of  the  beams,  ,  needed  is  calculated  using  the  following  equation. 
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The  gain  recovery  time  is  defined  by  the  balance  between  the  pumping  and  the 
deactivation  lifetime  of  KrF  .  The  gain  depletion  time,  r,  i  ^.  •  ,  which  corresponds  to 

the  time  that  the  input  pulse  can  extract  the  extra  energy,  is  obtained  by. 
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where  Tp  is  the  deactivation  lifetime  of  KrF*,  ^  is  the  input  beam  intensity,  and  ^  is  the 
saturation  intensity.  The  pulse  length  required  for  laser  fusion  is  comparable  to  the 
deactivation  lifetime  of  KrF*;  this  means  the  transient  dynamics  of  the  KrF*  medium  can  be 
utilized  to  reduce  the  number  of  optics  and  laser  beams. 
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Figure  5.  Scbenatics  of  the  transient  gain  dynamics  of  a  cw  amplifier  and  a  multiplexed  short  pulse 
amplifier  with  intervening  time  periods  between  pulses. 


Model  calculations  are  performed  and  compared  with  the  results  of  multiplexed  short 
pulse  experiments  reported  by  Shaw  et  The  code^^  used  for  this  calculation  contains  the  full  KrF 

kinetics,  including  vibrational  relaxation  and  KrF* (B) /(C)  mixing,  but  not  a  detailed  secondary  electron 
distribution  calculation;  however,  such  a  calculation  is  not  necessary  when  evaluating  the  performance  of 
this  novel  amplification  scheme.  Time-dependent  ASE  calculations  are  also  available  from  the  model.  The 
model  results  show  that  a  series  of  5-ns  pulses,  with  a  7-ns  interval  between  the  pulses,  can  obtain  an 
energy  gain  of  «1.7  tines  larger  than  that  obtained  by  a  series  of  12-ns  pulses  without  tine  intervals.  A 
peak  input  intensity  of  0.1  MV/cm^  is  assumed.  This  is  in  good  agreement  with  the  experimental  data. 


Next  analyzed  is  the  amplification  scheme  being  considered  for  the  Large  Aperture  Module  (LAM)  KrF 
amplifier  of  the  Aurora  system  at  Los  Alamos  National  Laboratory  (LANL) .  The  gain  length  is  2  n  with  20  cm 

unpumped  regions  at  the  ends,  the  excitation  rate  is  240  kW/cn'*,  and  the  laser  mixture  is  665  torr  Kr/Fg. 

The  input  pulses  to  the  LAM  are  amplified  by  a  double  pass  arrangement.  For  the  case  when  the  two-way 

propagation  time  of  the  laser  pulse  in  the  amplifier  is  longer  than  the  pulse  width,  the  gain  depletion 

caused  by  the  forward  and  backward  pulses  has  to  be  properly  distributed  in  the  amplifier  model  to  account 
for  depleted  regions  of  gas  in  the  amplifier.  The  gain  recovery  time  is  dependent  on  location  within  the 
amplifier  because  when  and  where  the  pulses  pass  through  a  region  of  gas  depends  on  the  pulse  length  and 
interpulse  time  interval.  The  net  effect  is  that  the  recovery  time  is  generally  shorter  than  the  interval 
between  the  pulses. 


Figure  6  shows  a  plot  of  the  extraction  efficiency  and  energy  gain,  with  respect  to  cw  amplification,  as 
a  function  of  the  tine  interval  between  pulses.  For  this  case,  the  pulse  width  of  the  input  pulse  is  5  ns. 

In  theje  calculations,  ASE  effects  have  been  neglected.  The  effects  of  ASE  are  calculated  separately  and  are 
given  later.  The  extraction  efficiency  decreases  with  increasing  time  interval.  At  a  5  ns  interval,  the 
reduction  of  the  efficiency  is  20%.  However,  with  a  5  ns  interval,  the  number  of  beams  needed  for 
multiplexing  is  5CWI  of  that  needed  for  the  normal  concatenated  multiplexed  system,  and  the  energy  gain  is 
three  times  higher. 
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Figure  6.  Plots  of  extraction  efficiency 
and  energy  gain  for 
multiplexed  short  pulse 
amplification  relative  to  cw 
amplification.  The  pulse 
width  of  the  input  laser  is 
5  ns . 


Figure  7  shows  the  effects  of  ASE  on  the  extraction  efficiency.  In  these  calculations,  the  gain  length 
is  shorten  to  40  cm  to  conserve  on  computation  time.  However,  the  gain-length  product  is  the  same  as  the  2- 
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Figure  7.  Decrease  of  extraction  efficiency  for  multiplexed  short  pulse  amplification  and  cw  amplification 
due  to  the  effects  of  ASE. 


device.  Since  the  propagation  time  in  the  40  cm  gain  length  is  shorter  than  the  input  pulse  width  of  5  ns, 
the  gain  recovery  time  is  longer  than  for  a  2-m  gain  length,  even  with  the  same  input  pulse  width  and 
interval.  Therefore,  Fig.  7  represents  a  worse  case  example  of  the  affects  of  ASE  build-up.  The  ASE 
significantly  degrades  the  extraction  efficiency  for  both  cw  and  short  pulse  amplifications.  The  percentage 
of  degradation  is  about  the  same  if  a  higher  input  intensity  is  assumed  during  short  pulse  amplification. 

Conclusion 

The  results  of  secondary  electron  density  measurements  and  e-beam  source  calculations  reveal  that  the 
energy  distribution  of  the  secondary  electrons  in  e  beam  pumped  laser  gases  is  not  Haxwellian.  Detailed 
Boltzmann  calculations  for  the  low  energy  electron  distribution  and  detailed  calculations  of  the  e-beam 
source  for  the  high  energy  distribution  are  essential  in  order  to  obtain  the  correct  rate  constants  for 
electron  reactions.  These  calculations  are  especially  important  when  determining  the  dissociative  attachment 
rate  constant  for  F2  and  predicting  the  number  density  of  the  secondary  electrons. 

The  transient  dynamics  of  the  KrF  gain  medium  permit  a  novel  multiplex  short  pulse  amplification  scheme 
to  be  used,  which  reduces  the  number  of  optics  and  laser  beams  needed.  Model  predictions  agree  well  with  the 
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experimental  results,  and  demonstrate  higher  energy  gain  than  the  usual  concatenated  multiplexing  schemes. 

The  extraction  efficiency,  during  the  condition  that  reduces  the  number  of  beams  by  501{,  is  approximately  30% 
loeer  than  the  usual  multiplexing  scheme  This  multiplexed  short  pulse  amplification  scheme  is  applicable 
even  for  the  amplification  of  picosecond  short  pulses. 
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Abstract 

Presented  are  the  results  of  electron  density  and  HCl  depletion  measurements  in 
electron-beam  excited  XeCl  laser  mixtures.  The  results  are  used  to  infer  the  importance 
of  different  kinetic  processes. 

Introduction 

There  is  a  body  of  experimental  evidence  suggesting  that  the  achievable  energy  density 
in  e-beam  pumped  XeCl  lasers  is  limited  by  cumulative  kinetic  processes.  Several 
laboratories  have  observed  limiting  efficiencies  of  the  laser  output  at  a  threshold  of  is5 
J/2..  There  is  currently  no  consensus  in  the  field  concerning  the  mechanisms  causing 

this  effect.  Interpretation  of  these  observations  is  difficult  since  there  is  little 
experimental  data  on  the  concentration  changes  of  key  species  during  e-beam  excitation. 

As  part  of  a  continuing  program  to  characterize  and  understand  the  reaction  kinetics 
occurring  in  these  systems,  we  have  performed  measurements  to  determine  the  time-dependent 
electron  density  and  extent  of  halogen  consumption  in  e-beam  pumped  XeCl .  This  paper 
reports  the  first  measurements  of  HCl  consumption  in  e-beam  pumped  XeCl  laser  mixtures 
under  non-lasing  conditions.  The  results  of  the  electron  density  measurements  have  also 
been  reported  el sewhere . ^ ^ 

Experimental  Results  and  Discussion 

The  electron  beam  source  used  in  both  studies  is  a  carbon  felt,  cold  cathode  source 
operating  at  ^300  kV  and  delivers  >10  A/cm*  at  the  foil.  To  better  characterize  the  gain 
length  in  the  gas  mix,  the  e-beam  was  baffled  down  prior  to  the  foil  with  metal  plates. 
Power  deposition  in  the  gas  is  ol80  kW/cm*  and  is  measured  using  a  pressure  rise 
technique.^  For  all  gas  mixtures  presented  in  this  study,  the  halogen  donor  is  HCl,  the 
diluent  is  neon,  and  the  xenon  concentration  is  1.5%.  The  total  gas  pressure  is  3000  torr 
and  the  initial  temperature  300  K. 

Electron  Densities 

The  experimental  setup  used  to  determine  electron  densities  is  shown  schematically  in 
Fig.  1.  The  gas  chamber  of  the  Tahoma  laser  system^  is  placed  in  one  arm  of  a  10.6  jim 
quadrature  interferometer.^  The  laser  mirrors  are  removed  and  the  fused  silica  windows 
replaced  with  BaF2 .  The  COg  probe  laser  passes  ®!3  cm  from  the  foil,  centered  on  the 
nominal  lasing  axis.  As  the  index  of  refraction  of  the  plasma  is  proportional  to  electron 
number  density,  the  measured  fringe  shifts  during  the  pulse  directly  yield  the  temporal 
changes  in  the  electron  density. 

Figure  2  shows  the  temporal  evolution  of  electron  density  for  a  variety  of  initial  HCl 
concentrations.  The  dashed  line  at  the  bottom  of  the  figure  depicts  a  normalized  typical- 
current  wavefo’^m.  During  the  first  50  nsec  of  the  pulse,  the  electron  density  rises 
sharply  then  gradually  decreases.  Since  the  dissociative  attachment  rate  constants  for 
HCl  increase  dramatically  with  vibrational  excitation,  this  behavior  can  be  attributed  to 
the  production  of  v i brat i onal 1 y~ exc i ted  HCl  during  the  pumping  process. 


V.A-V- 


tf*.  ■ 


.  •  »  ■ 


POLARIZER 


z 

til 


<  2 
til  fC 


Figure  2.  Electron  densities  of  XeCl 

mixtures  for  various  initial 
HCl  concentrations.  The 
dashed  curve  is  a  normalized 
trace  of  a  typical  electron 
beam  current  waveform.  Data 
points  are  not  shown;  however, 
a  representative  error  bar  is 
shown . 


This  initial  peak  in  electron  density  is  also  observed  in  model  predictions  for  e-beam 
discharges  under  these  conditions . ^ The  surprising  result  is  the  magnitude  of  the 
electron  density,  which  is  2-4  times  greater  than  the  model  predictions.  Peak  laser 
performance  on  Tahoma,  for  this  pulse  duration,  is  obtained  with  HCl  concentrations  of 
The  reason  leaner  mixtures  do  not  perform  as  well  is  evident  from  the  electron 
density  results.  At  lower  HCl  concentrations,  the  electron  density  is  found  to  increase 
rapidly  at  late  times  during  the  pulse,  implying  that  the  HCl  fuel  has  been  exhausted. 
Halogen  burnup  at  low  HCl  concentrations  is  also  evident  in  the  XeCl*  sidelight 
fluorescence  shown  in  Fig.  3. 

HCl  Depletion 

The  experimental  setup  is  schematically  represented  in  Fig.  4.  FlCl(v,J)  densities  are 
measured  directly  by  infrared  absorption  on  the  first  vibrational  overtone  of  v"=0,l,  and 
2.  All  measurements  were  made  on  the  P(3)  rotational  transition  and  rotational 
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Figure  3.  XeCl  sidelight  fluorescence 
waveforms  for  the  same 
mixtures  shown  in  Fig.  2. 


200  400 

TIME  (NSEC) 


REFERENCE 


RATIOMETER 


PROBE 


Nd:YAG  PUMPED 
DYE  LASER 


TAHOMA 
LASER  CHAMBER 


Hg  RAMAN  CELL 


50%  BEAM 
SPLITTER 


ttttttt  t 

ELECTRON  BEAM 


Figure  4.  Schematic  of  the  experimental  set  up  for  the  halogen  depletion  studies. 

equilibration  at  a  final  gas  temperature  of  “400  K  Is  assumed.  To  generate  tunable  light 
around  1.8  jim,  the  output  of  a  Nd ; YAG  pumped  dye  laser  was  Raman-shifted  in  a  high 
pressure  H2  cell.  The  desired  Stokes  component  was  isolated  and  directed  through  the 
nominal  active  zone  of  the  Tahoma  laser  chamber.  A  portion  of  the  probe  beam  was  split 
off  prior  to  the  chamber  for  normalization  purposes.  The  transmitted  power  was  measured 
with  a  pyroelectric  detector  and  ratioed  to  the  reference  power  to  yield  the  relative 
amount  of  absorption . 

From  measurements  with  static  gas  fills  of  known  HCl  concentrations,  a  significant 
deviation  from  a  simple  exponential  absorption  was  observed.  Figure  5  shows  the  measured 
HCl  number  density  as  calculated  assuming  an  exponential  absorption  and  using  published 
overtone  absorption  cross  sections.^  This  behavior  was  reproducible  throughout  the 
measurements  and  is  attributed  to  saturation  effects.  A  polynomial  fit  to  these  data 
points  was  made  and  used  to  correct  the  measurements  taken  during  e-beam  excitation. 

The  presence  of  a  strong  absorption  during  the  e-beam  pulse  prevented  the  time- 
dependent  loss  rate  of  HCl  from  being  determined.  This  absorption  is  structureless  within 
our  resolution  (“0.3  cm  and  was  found  to  be  strongest  in  the  absence  of  MCI.  This 
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Figure  5.  Measured  HCl  number  densities 
as  determined  assuming  an 
exponential  absorption.  The 
staight  line  shows  the 
expected  behavior  and  the 
curve  is  a  polynomial  fit  to 
the  data  points. 
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signal  decayed  quite  rapidly  after  termination  of  the  e— beam,  allowing  the  residual  HCl 
density  to  be  measured  2  /isec  after  the  excitation  pulse.  The  identity  of  this  absorbing 
species  could  not  be  accurately  determined,  but  the  lack  of  spectral  structure,  the 
strength  of  the  absorption,  and  the  short  lifetime  suggest  an  assignment  of  Xcg- 

Figure  6  gives  the  results  of  the  HCl (v=0)  depletion  measurements.  For  initial 
concentrations  below  SO.IR  HCl  the  data  show  nearly  complete  removal  of  the  HCl.  This  is 
consistent  with  the  electron  density  measurements;  for  these  lower  HCl  concentrations,  the 
HCl  is  depleted  at  some  point  during  the  pulse  and  this  absence  of  an  attacher  allows  the 
electron  density  to  increase.  At  higher  initial  HCl  concentrations,  the  fractional  burnup 
does  not  level  off  but  continues  to  increase.  At  the  highest  initial  concentrations 
studied  (0.2055  HCl),  the  fractional  loss  of  HCl  (v=0)  is  found  to  be  61*855.  Higher 
vibrational  states  were  observed  only  at  the  higher  HCl  densities.  At  0.2055  initial  HCl, 
the  population  in  HCl(v=l)  shows  a  maximum  of  2.455  of  the  HCl  added.  HCl (v=2)  was  never 
observed,  allowing  an  upper  limit  of  0.255  to  be  determined  from  our  noise  level  for  this 
transition . 


Figure  6.  Measured  HCl (v=0)  loss  after 
firing  the  electron  beam  for 
various  initial  HCl 
concentrations.  The  straight 
line  shows  the  expected 
behavior  for  10055  consumption 
of  the  HCl  present. 
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Quantitatively,  these  results  represent  a  minimum  cumulative  loss  of  HCl.  Since  the 
probe  pulse  is  delayed  by  2  ;isec,  some  degree  of  HCl  regeneration  will  have  occurred  via 
three-body  recombination  and  associative  detachment  reactions .  The  residual  populations 
in  HCl(v=l,2)  also  must  be  regarded  as  a  minimum  since  vibrational  relaxation  by 
superelastic  collisions  with  electrons  is  expected  to  be  significant. 

Concluding  Remarks 

The  kinetic  processes,  involved  in  e-beam  pumped  XeCl  lasers,  represent  a  rich 
chemistry  that  at  present  is  still  incompletely  understood.  The  data  presented  here  show 
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both  a  significantly  higher  electron  density  and  a  greater  HCl  loss  than  expected  from 
model  predictions  for  this  system.  These  results  suggest  that  electron  associative 
detachment  cannot,  by  itself,  account  for  the  depletion  of  HCl  and  that  other  loss 
channels  must  be  considered.  A  likely  candidate  is  direct  electron  impact  dissociation  of 
HCl  leading  to  neutral  atomic  products.  Further  analysis  of  the  data  is  currently  in 
progress . 
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Xenon  excited  state  density  measurements  in  electron  beam  pumped  XeCI 
laser  mixtures 

F.  Kannari,  W.  D.  Kimura,  J.  F.  Seamans,  and  Dean  R.  Guyer 
Spectra  Technology,  Inc.,  2755  Sorthup  fVay,  Bellevue,  IVashinglon  98004-1495 

(Received  1 1  August  1987;  accepted  for  publication  12  October  1987) 

Time-dependent  density  measurements  of  the  lowest  xenon  excited  states  (5/>^6r)  in  electron 
beam  (e-beam)  pumped  XeCl  laser  mixtures  (nonlasing)  are  performed  using  pulsed  hook 
interferometry.  The  e-beam  pulse  length  is  =;0.45  /rs  (full  width  at  half-maximum)  with  an 
average  excitation  rate  of  =;250  kW/cm  ’.  Density  differences  (AJV*)  of  transitions  at  823.2 
nm  (6s[3/2]5-6p[3/2]j),  828.0  nm  (&s[3/2]?-6p[  l/2]o),  and  840.9  nm  (6ls[3/2]5- 
6p[3/2] , )  are  obtained  for  various  HCI  and  Xe  concentrations.  For  a  98.3%  Ne/1.5% 

Xe/0.16%  HCI  mixture  at  3000  Torr,  A/V*  (823.2  nm)  and  AA'*( 828.0  nm)  are  relatively 
constant  at  =3  X  10‘*  and  =:  1.7x  10‘^  cm~’,  respectively.  At  lower  initial  HCI 
concentrations,  the  AN  *  ( 823.2  nm )  density  starts  out  similar  to  the  0. 1 6%  HCI  case,  but  tends 
to  increase  dramatically  during  the  e-beam  pulse.  For  a  0.04%  HCI  mixture,  a  AN  *  ( 840.9  nm ) 
density  of  =;4.5  X  10”  cm“^  is  measured  at  the  end  of  the  e-beam  pulse.  Preliminary 
comparisons  of  the  data  with  a  computer  model  show  good  a,^reement  for  HCI  concentrations 
>0.16%,  but  disagreements  at  leaner  concentrations. 


The  formation  of  XeCl  )  in  electron  beam  (e- 

beam )  pumped  lasers  has  been  the  topic  of  considerable  de¬ 
bate.  '  ’  This  letter  describes,  to  our  knowledge,  the  first  time- 
dependent  measurements  of  xenon  excited  state  (Xe*) 
densities  in  e-beam  pumped  XeCl  mixtures  (nonlasing). 
These  results  together  with  earlier  measurements  of  the 
time-dependent  electron  density’  and  the  HCI  depletion"*  in 
XeCl  were  performed  to  help  improve  our  understanding 
and  modeling  of  XeCl  formation  and  laser  kinetics. 

Excited  states  of  xenon  are  important  for  the  neutral 
channel  formation  of  XeCl.  The  first  excited  levels  of  xenon, 
Xe*(5p’[’P?,,  ]6i),  have  very  small  XeCl  formation  rate 
constants  with  vibrational  ground-state  HCI’ ;  however,  the 
rate  constants  increase  dramatically  with  vibrational  excita¬ 
tion  of  HC1.‘ 

The  Xe*  density  is  measured  using  hook  interferome¬ 
try^  in  which  the  e-beam  pumped  laser  medium  is  placed  in 
one  branch  of  a  laser  interferometer  ( see  Fig.  1 ) .  An  excimer 
pumped  broadband  dye  laser  output  passes  twice  through 
the  gas  chamber  3  cm  from  the  foil  and  centered  around 
the  nominal  lasing  axis.  The  output  of  the  interferometer  is 
directed  into  a  spectrometer  whose  output  is  an  interference 
pattern  that  contains  “hooks”  caused  by  the  anomalous  dis¬ 
persion  around  atomic  absorption  lines.  Interference  pat¬ 
terns  are  detected  by  a  vidicon  camera  and  recorded  on  vid¬ 
eotape  for  later  analysis.  The  e-beam,  excitation  rate  is  250 
kW/cm’ ,  which  is  higher  than  during  the  earlier  electron 
density  measurements.’ 

The  wavelength  separation  A  between  the  hooks  in  the 
interference  pattern,  centered  around  a  transition  wave¬ 
length  /l,y  (where  i  and  j  are  the  lower  and  upper  levels, 
respectively),  is  given  by^* 


Ilf, 

ttK  ' 


(1) 


where  N,^  are  the  population  densities,  g,j  are  the  statistical 
weights  of  each  level,  r,,  is  the  classical  electron  radius,  /  is  the 
length  of  the  medium,  is  the  oscillator  strength  of  the 


transition,  and  AT  is  a  constant  that  can  be  determined  from 
an  undistorted  interference  pattern.  A  knowledge  of  a  par¬ 
ticular  transition  oscillator  strength  allows  the  population 
difference,  AN*  =  iV,  —  (g,/gy)A^y,  to  be  determined  from 
the  hook  spectra. 

Hook  interferometry  offers  several  advantages.  It  has  a 
relatively  large  dynamic  range,  the  results  are  insensitive  to 
line  shape  perturbations  such  as  pressure  broadening,  and 
complicated  geometrical  or  detector  calibrations  are  unnec¬ 
essary. 

Figure  2  is  an  energy-level  diagram  showing  the  transi¬ 
tions  examined  during  the  experiment:  823.2  nm  (6is[  3/2)5 
-6/>[3/2]2),  828.0nm  (6s[3/2]?-6;>[  l/2]o),  and  840.9  nm 
(6ij(3/2]5-6/>[3/2]  i).  Oscillator  strengths  for  these  lines 
are  0.23,  0.12,  and  0.012,  respectively.’  The  total  optical 
path  of  the  excited  laser  medium  is  70  cm.  The  value  of  AT  in 
Eq.  ( 1 )  is  typically  =:  1 .5  x  10’ .  Thus,  on  our  system  a  hook 
separation  of  1  nm  at  823.2  nm  corresponds  to  a  AN  *  density 
of 7x10”  cm-’. 

A  typical  hook  interference  pattern  is  shown  in  Fig.  3. 
As  can  be  seen,  not  all  the  transitions  are  well  separated  from 
one  another.  In  our  case,  two  relatively  strong  transitions  are 
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FIG  1.  Schemalic  of  xenon  excited  state  density  measurement  system 
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FI(i.  2  Friergy-leve!  diagram  df  xenon  showing  transitions  examined 

observed  at  826.6  and  822. 1  nm  near  the  828.0  and  823.2  nm 
transitions.  The  transition  at  826.6  nm  corresponds  to  Xe 
(6.s'(  1/2  |"-6/7'[  1/2] , ).  The  transition  at  822.1  nm  (not 
shown  in  Fig.  3)  has  not  been  identified.  These  lines  distort 
the  hook  profiles  and  position  on  the  short  wavelength  side 
ofthe  828.0and  823.2  nm  transitions,  while  theeffectson  the 
long  wavelength  side  ofthe  transitions  are  negligible. 

For  an  isolated  transition,  the  hook  separation  is  sym¬ 
metric  about  the  transition  wavelength.  The  quadratic  Stark 
effect  can  displace  the  hooks  towards  long  wavelengths,  '  but 
because  the  electron  density  is  relatively  low  '  (  ^  10'' 
cm  '  ),  this  effect  can  be  ignored.  Therefore,  during  analy¬ 
sis  of  the  data,  it  is  assumed  the  total  hook  separation  can  be 
calculated  by  doubling  the  distance  between  the  hook  on  the 
long  wavelength  side  and  line  center  (see  Fig.  3). 

In  ty  pical  laser  kinetics  models,  Xe**  is  a  generic  highly 
excited  state  representing  Xe(6/7)  states  above  Xe(6,v).  Cal¬ 
culations  indicate  that  the  total  population  of  Xe**  is  signifi¬ 
cant.  Since  cross  sections  for  stepwise  excitation  of  Xe*  are 
estimated  to  be  large,"'  the  two-body  relaxation  rates  of 
Xe**  by  neutral  molecules  are  small,' '  and  the  electron  den¬ 
sity  is  appreciable,'  thermalization  among  these  levels  at 
near  the  electron  temperature  of  ~  1.5  eV  (Ref  12)  is  ex¬ 
pected.  Within  this  assumption,  the  upper  state  populations 
of  the  probed  transitions  are  estimated  to  be  ^  l()-369f  of 
the  lower  state.  Measurements  of  Xe(6.s'[  1/2  |.)  in  XeCl, 
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FIG.  4.  Time-dej>cndenl  A.V*  densities  of  a  30(X)  Torr  mixture  of  3G 
Ne/F5^  Xc/0. 16^  HCl  at  294  K.  Closed  and  open  circles  eorresp<ind  m 
A.V  •  ( 823.2  nm)  and  A.V  •  ( 828,0  nm  ).  respect  nets .  The  dashed  curse  is  a 
normalized  trace  of  the  e-bcam  current  waveform. 


XeF,  and  Ne/Xe  mixtures  also  suggest  significant  upper 
state  densities."  Due  to  the  uncertainty  in  the  upper  state 
population,  all  measurements  herein  are  presented  as  den¬ 
sity  differences  (A.V*)  rather  than  level  densities. 

Figure  4  shows  typical  experimental  results  obtained  at 
the  828.0  and  823.2  nm  transitions,  for  an  initial  HCl  con¬ 
centration  of  0. 1 6fT.  The  estimated  measurement  uncertain¬ 
ty  associated  with  all  the  data  points  presented  in  this  letter  is 
+  lOi^f .  As  can  be  seen,  the  density  differences  are  relatively 
constant  during  the  e-beam  pumping. 

Xenon  excited  states  are  produced  primarily  by  recom¬ 
bination  reactions  of  dimer  ions  ( Xe,'  ,  NeXe  '  )  rather  than 
direct  excitation  of  ground-state  Xe  by  electrons.  Cl  from 
HCl  controls  the  Xe*  population  by  intercepting  Xe ' ,  the 
precursor  to  the  dimer  ions,  to  form  XeCI.  Hence,  as  long  as 
dissociative  attachment  is  the  dominant  electron  loss  pro¬ 
cess.  the  Xe*  density  decreases  with  increasing  initial  HCl 
concentration.  In  fact,  experimental  results  for  0,32'T-  HCl 
(not  shown)  are  ~47^  of  A.V*  at  0.16'^  HCl.  However, 
Figs.  5  and  6  show  time-dependent  A.V  ♦  densities  for  rela- 
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FIG.  6.  Timc-depcndcnt  A.V  •  densities  for  an  initial  IICI  concentration  of 
0  ()4^r  The  other  gas  parameters  and  symbols  are  the  same  as  in  Figs  4 
and  5 


lively  low  initial  HCl  concentrations  of  0.08  and  0,04'^,  re¬ 
spectively.  In  these  figures,  data  for  A.V  *  ( 840.4  nm  )  are  also 
shown,  particularly  for  densities  >  10'“  cm  'where  broad¬ 
ening  of  the  transition  linewidth  obscures  the  hooks  for  the 
828.2  and  828.0  nm  lines.  Two  interesting  Xe*  density  beha¬ 
viors  are  observed.  First,  although  higher  Xe*  densities  are 
expected  for  lower  HCl  concentrations,  densities  measured 
during  the  first  200  ns  are  nearly  the  same  as  given  in  Fig.  4. 
Second,  at  low  HCl  concentrations,  the  A.V*  density  dra¬ 
matically  increases  at  a  certain  point  during  the  e-beam 
pulse.  As  the  initial  HCl  concentration  is  reduced,  the  A.V  * 
density  begins  to  increase  at  earlier  times.  A  peak 
A.V  *(  840.9  nm )  density  of  4.5  X  lO'^  cm  '  is  measured  for 
the  0.04%  HCl  mixture.  The  abrupt  increase  in  A.V  *  densi¬ 
ties  is  likely  caused  by  exhaustion  of  HCl  during  the  pump 
pulse,  which  leads  to  an  abrupt  increase  in  the  secondary 
electron  density '  and  results  in  an  increase  in  electron-dimer 
ion  recombination. 

Measurements  also  show  that  the  Xe*  density  increases 
with  increasing  initial  xenon  concentration.  At  150  ns  into 
thee-beam  pulse,  A.V  *(  823.2  nm )  densitiesof  2.0  «  10''*  and 
4  5  ■  10'^  cm  'are obtained  for0,5%  and  3.0%  Xe,  respec- 
tivelv,  with  HCl  fixed  at  0.16%. 
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Preliminary  comparisons  with  a  computer  model 
yield  good  agreement  with  the  data  for  initial  HCl  concen¬ 
trations  ;^-0. 16%.  At  leaner  concentrations  the  mtxlel  does 
not  predict  the  rapid  increase  in  Xe*  density  seen  during  the 
pulse.  This  may  be  related  to  the  model’s  underprediction  of 
HCl  depletion, ■*  The  model  also  does  not  predict  the  k'ng 
Xe*  density  lifetime  observed  after  the  end  of  the  e-beam 
pulse  for  lean  HCl  mixtures.  This  lifetime  appears  correlated 
with  the  electron  density, '  and  is  sensitive  to  the  NeXe*  and 
NeXe  '  formation  rates,  which  may  mean  that  an  improved 
estimation  of  these  rates''  is  needed.  Hence,  the  Xe*  density 
measurements  demonstrate  that  our  understanding  of  Xe* 
formation  is  gtwd  except  under  conditions  when  a  signifi¬ 
cant  amount  of  halogen  is  lost.  Additional  discussion  regard¬ 
ing  the  implications  of  these  measurements,  as  well  as  .Xe* 
and  Xe**  measurements  in  XeF',  will  be  presented  in  another 
paper. ' ' 
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ABSTRACT 


Sararal  aoral  diacaoatic  tachaiquaa  hara  baaa  aaad  to  aaalyaa  tka  laaar  kiaaties 
occurriag  ia  alactroa-baaa  punpad  axciaar  laaara.  Oaocribad  ara  a  quadratura  iatarfaroaatar 
for  aaaauriae  tiaa-dapaadeat  alaetroa  daaaitiaa,  aa  ovartoaa  abaorptioa  tackaiqua  for 
aaaauring  tka  HCl  concaatratioa ,  a  broadband  dya  laaar  aystaa  for  aaaauriag  tka  apactrally 
raaolvad  gain  in  XaF,  and  a  kook  intarfaroaatar  for  aaaanring  tiaa-dapandant  axcitad  atata 
daaaitiaa.  Tkaae  tackniquaa  can  alao  ba  appliad  to  otkar  non-axciaar  laaara. 

1.  INTRODUCTION 

Furtkar  advancaaant  of  axciaar  laaar  tacknology  raliaa  on  tka  furtkar  adraaceaant  of  our 
undaratanding  of  tka  laaar  kinatica.  Modala  of  tka  laaar  parforaaaca,  wkick  ara  uacd  for 
axaapla  to  acala  tkaaa  daricaa  to  larga  aiaaa,  ara  only  aa  good  aa  our  undaratanding.  Our 
coafidanca  in  tkaaa  aodala  ia  only  aa  good  aa  tka  data  arailabla  to  ralidata  tka  nodala. 
Henca ,  tkara  ia  a  aaad  to  apply  appropriata  diagnoatic  aaaauraaanta  on  tkaaa  laaara  to 
ravaal  fundaaantal  ckaractariatica  of  tka  laaar  aadiua. 

Uaual  diagnoatica  tkat  aaaaura  paraaetara  auck  aa  tka  laaar  output  aaargy,  afficiancy, 
pulaa  ahape,  and  lasing  apectrun,  only  prorida  aacroacopic  information  ragarding  tka  laaar 
aadiua  ckaractariatica .  Tka  data  cannot  aaaily  ravaal  dataila  about  tka  gaa  kinatica 
bacauae  tkay  rapraaent  aaaantially  tka  nat  raault  of  many  complax  procasaaa  occurring 
aimultanaoualy .  Gain,  abaorption,  and  aidalight  fluoraacanca  maaauramanta  ravaal  mora 
dataila,  but  they  are  atill  tka  nat  affect  of  a  complicated  molecular  ckamiatry  occurring 
within  the  laser  plasma.  Tkarafore,  to  understand  tka  kinatica  on  a  fundaaantal  level,  it 
in  necansary  to  measure  specific  quantities  on  a  molecular  level.  It  is  at  this  molecular 
level  tkat  tka  kinetics  models  ara  baaed,  and  it  ia  at  this  laval  tkat  critical  data  ia 
lacking. 

This  paper  dincunnan  diagnoatic  tackniquaa  used  at  Spectra  Tacknology,  Inc.  (STI)  on  an 
alectron-baam  (a-baam)  pumped  axcimer  laser^  to  measure  various  ckaractariatica  of  the  laaar 
active  medium.  Tka  first  diagnostic  ia  an  infrared  quadratura  inter farometar  to  measure  the 
electron  density  in  XaCl,^  XeF,  and  KrF;^  the  second  is  a  Raman-nkif tad  Nd:YAG  pumped  dye 
laser  to  maaaura  tka  HCl  depletion  in  XaCl;*  tka  third  is  aa  axcimer  pumped  broadband  dya 
laser  aystam  to  measure  apactrally  rasolved  small  signal  gain  in  XeF;^  ud  the  fourth  ia  a 
kook  intarfarometar  to  determine  the  Xa  axcitad  atata  danaitiaa  in  XaCl°''  and  XaF.~ 

8.  DESCRIPTION  OF  B-BEAM  PUMPED  LASER 


Tka  laser  (Takoma)  during  all  these  meaauramenta  baa  been  used  axtanaivaly  in  the  past 
for  lasing  experiments  on  XeCl ,  XeF,  KrF,  and  other  axcimers.  Details  of  tka  laser 
aystam  can  ba  found  in  the  aforementioned  rafarancaa.  Bowevar,  all  tka  maasuranants 
described  in  tkia  paper  ara  performed  under  nonlasing  conditions  (tkat  is  tka  normal  laser 
mirrors  ara  removed).  Tka  maasuramenta  can  alao  ba  performed  under  lasing  conditions,  with 
tba  various  probe  laser  beams  introduced  off-axis  to  the  active  volume,  and,  indeed,  this 
would  ba  vary  useful  to  perform  at  some  future  time. 

The  nominal  gain  length  is  70  cm,  but  it  can  ba  shorten  by  blocking  tka  ends  of  the  a- 
beam  with  baffles.  This  ability  to  shorten  the  gain  length  and  tkaraby  tuna  the  sensitivity 


.  .-.v.  ■r. 


of  til*  ■•Murenent  i*  loportaat  b*c*u**  occMiosally  th*  cliajic**  induced  on  thm  prob*  Inner 
benne  «•»"  be  beyond  the  dynaaic  rnn(*  of  th*  detectorn. 


The  *-b*na  pule*  length  ie  rnrinbl*  froa  «100  ne  to  1  pe.  Generally,  th*  tiae-dependence 
of  th*  quantity  of  interest  relative  to  th*  *-b*aa  pule*  shape  is  desired.  However,  rather 
than  perform  all  measurements  at  the  longest  s-beaa  puls*  length,  aeanuremente  acs  .al  io  don* 
at  MO. 6  ps.  This  is  because  the  optimum  halogen  concentration  for  lasing  tends  to  increase 
with  energy  loading  (pump  rate  times  *-beaa  pulse  length);  however,  higher  halogen 
concentration  also  increases  the  amount  of  halogen  cfuenching.  Therefore,  for  a  fixed  pump 
rat*  and  halogen  concentration,  it  is  necessary  to  vary  the  s-beaa  puls*  length  in  order  to 
investigate  other  operating  regimes  of  interest. 

Timing  of  th*  various  lasers  and  detectors  msed  in  th*  diagnostics  is  controlled  by 
separate  digital  delay  generators  triggered  by  a  coaputer  program  that  primes  and  fires  the 
laser,  and  acquires  the  data.  A  video  frame  store  unit  in  th*  computer  system  captures  the 
output  from  the  vidicon  camera  used  during  some  of  the  experiments,  and  provide*  a 
convenient  method  for  storing  data  for  later  analysis. 


3.  ELECTRON  DENSITY  MEASUREMENT  USING  gUADRATURE  INTERFEROHETRY 


Figure  1  is  a  schematic  of  th*  quadrature  interferometer  system  used  to  measure  th*  time- 
dependent  electron  density.  The  quadrature  interferometer  unes  a  single  mode,  cw  CO^ 

laser  to  probe  the  active  volume  of  the  laser  chamber  positioned  in  one  branch  of  a 
Michelson  interferometer.  Double  use  of  th*  interferometer  central  beamsplitter  together 


with  a  X/4  plate  provide  two  output  beams  that  are  00*  out  of  phase  with  on*  another.  Th* 
"cat’s  eye"  retro-reflectors  permit  the  incident  and  reflected  beasM  to  be  offset  from  on* 
another,  while  at  the  same  tine  not  disturbing  th*  polarisation. 


CO; 

LASER 


Figure 


1. 


Schematic  of 
quadrature 
interferometer 
system  for 
measuring  the 
tine-dependent 
slectron  density. 


Liquid  nitrogen-cooled  HgCdT*  detectors  measure  the  10.0  pa  signals.  Du*  to  infrared 
emission  from  the  laser  chamber,  th*  detectors  are  placed  *18  a  away  to  eliminate  any  false 
signals.  (Narrowband  filters  centered  at  10.6  pa  would  have  also  solved  this  problem.) 

The  electron  density  is  determined  by  measuring  the  fringe  shifts  occurring  during^the  e- 
beam  excitation  and  calculating  the  density  using  the  plasma  dispersion  relationship.  For 
this  experiment,  because  of  th*  fairly  high  electron  density  and,  therefore,  large  number  of 
fringe  shifts,  the  *-b*am  is  baffled  to  35  cm  in  length  so  that  th*  00^  laser  beam  travels 
through  70  cm  of  active  medium.  One  fringe  shift  then  corresponds  to  an  electron  density  of 
3  X  10^^  em^. 

An  example  of  th*  two  detector  output*  is  shown  in  Fig.  S.  Note  th*  rapid  fluctuations 
of  the  signal  during  the  beginning,  indicating  a  rapid  increase  of  electron  density.  The 
signal  varies  more  slowly  after  MIOO  ns,  indicating  a  more  constant  density,  and  then  ends 
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Figura  2.  Oacilloacope  traces  from  the 
HgCdTe  detectors.  The  top 
trace  is  00*  out  of  phase 
with  the  bottom  one.  The 
gas  mixture  is  08.45(  Ne/l.BX 
Xb/0.08%  HCl  at  3000  Torr 
total  pressure  and  an  e-beam 
deposition  rate  of 
•180  kW/cm~®. 


with  a  relatiTely  steadp  rariation  indicating,  in  this  case,  a  steady  increase  of  density. 
Although  all  the  phase  information  can  be  ascertained  from  only  one  detector,  the  second 
detector  eliminates  the  ambiguity  that  can  arise  when  the  first  detector  signal  passes 
through  an  extremum.  At  an  extremum,  the  sign  of  the  phase  shift  is  unknown;  bowerer,  by 
examining  the  other  detector,  which  will  be  at  its  maximum  differential  sensitivity,  the 
direction  of  the  phase  shift  is  known.  An  example  of  the  electron  density  results  for  XeCl 
are  given  in  Pig.  3. 
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Figure  3.  Typical  electron 
deneities  of  XeCl 
mixtures  (1.5X  Xe, 
balance  Ne  to  3000 
Torr)  with  the  e- 
beam  deposition 
S180  kW/cm®.  The 
dashed  curve  is  a 
normalised  trace 
of  a  typical 
electron  beam 
current  waveform. 
Data  points  axe 
not  shown; 
however,  a 
representative 
error  bar  is 
indicated . 


4.  HCl  DEPLETION  MEASUREMENT  USING  OVERTONE  ABSORPTION 

Balogcn  depletion  in  excimer  lasers  during  the  pump  pulse  is  a  critical  issue.  To 
measure  the  HCl  concentration,  the  technique  chosen  is  direct  absorption  measurements  on  th 
first  vibrational  overtone  of  v"b0,  1,  and  2  at  *1.8  /m.  This  is  because  HCl  absorbs  too 
strongly  on  the  v’^l  '*  v*b0  transition  to  permit  convenient  measurement  over  the  length  of 
the  active  volume.  All  measurements  are  made  on  the  P(3)  rotational  transition  and 
rotational  equilibration  at  a  final  gas  temperature  of  *400  K  is  assumed. 

Figure  4  is  schematic  of  the  halogen  depletion  measurement  system.  The  laser  medium  in 
probed  by  a  Nd:YAG  pumped  dye  laser  whose  output  is  Baman-shif ted  to  1.8  pm  in  a  high 
pressure  H2  cell.  Ratioing  the  signals  from  the  reference  and  probe  pyroelectric  detectors 
gives  the  relative  amount  of  absorption.  For  these  measurements,  the  gain  length  in  70  cm. 


An  nxcimer  pumped  dye  Imaer  (TUQ  dye) ,  operating  broadband  (*50  A  wide)  and  centered  at 
around  350  nm,  ia  ueed  to  obtain  tbe  full  wpectrum,  waall  signal  gain  of  e-beam  excited  XeF 
during  a  single  shot.  To  our  knowledge,  this  is  the  first  time  tbe  spectrally  resolred  gain 
spectrum  baa  been  measured  in  XeF.  Although,  as  mentioned  earlier,  gain  measurements  are 
not  usually  simply  related  to  basic  gas  kinetic  quantities,  being  able  to  examine  the  entire 
gain  spectrum  does  reweal  information  regarding  tbe  degree  of  wibrational  and  rotational 
coupling  in  tbe  XeF  manifold.  Such  information  is  particularly  important  when  attempting  to 
understand  tbe  narrowband  extraction  characteristics  of  XeF. 

A  schematic  of  the  gain  measurement  system  is  shown  in  Fig.  6.  The  dye  laser  output  is 
split  into  a  probe  and  reference  beam.  Within  the  laser  chamber,  the  probe  intensity 
averaged  over  the  entire  XeF  gain  band  is  <60  kW/cm^.  The  dye  laser  pulse  length  is  <6  ns; 
hence,  by  varying  the  time  delay  of  the  dye  laser  with  respect  to  the  e-beam,  the  temporal 
dependence  of  the  gain  during  the  e-beam  pumping  can  be  obtained . 
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Figure  6. 


SPECTROMETER 


SchsBatic  of  (sin 
■•ssursBcnt  ■ystea 
for  obtsinisg 
•psctrslly 
rssolTsd  (sin  in 
XsF. 


The  probs  snd  reference  besaa  sre  eent  through  cslibrsted  aeutrsl  density  filters  and  are 
▼erticslly  displaced  across  the  entrance  slit  to  a  1  b  spectroaeter .  Hence,  at  the  output 
of  the  spectroneter  (used  without  an  exit  slit)  are  two  bands  of  light,  one  aboTs  the  other, 
that  correspond  to  the  dye  laser  spectruB.  The  probe  besa  is  saplified  by  the  excited  XeF 
BediuB;  taking  the  ratio  of  the  probe  to  the  reference  then  yields  the  gain. 

A  UV  widicon  esaers  detects  both  the  probe  and  reference  spectra  siBultaneously .  This  is 
shown  scheastically  in  Fig.  7.  A  coaputer-controlled  wideo  frsae  store  unit  captures  the 
frsae  with  the  spectra  and  stores  the  data  in  eoeputer  BSBory.  By  aligning  the  eaaera  with 
its  raster  lines  oriented  along  the  dispersion  direction  of  the  spectra,  the  raster  lines 
within  a  band  of  light  can  be  aTeraged  pixel-to-pixel,  thereby  iaproring  the  signal-to-noise 
of  the  data.  This  aethod  of  data  collection  proxides  a  convenient  Beans  to  acquire  and 
analyse  spectral  data. 
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Figure  7.  Data  acquisition 
systSB  for 
collecting  output 
froB  the  gain 
spectrooeter  shoim 
in  Fig.  6. 


A  typical  exsaple  of  the  small  signal  gain  Beasured  is  given  in  Fig.  8.  The  two  dips  in 
the  gain  spectrua  at  351.62  and  352.05  na  correspond  to  neon  absorption  lines^^  and  are  used 
to  provide  an  accurate  absolute  wavelength  calibration  of  the  spectra. 

e.  XENON  EXCITED  STATE  MEASUREMENTS  USING  HOOK  INTERFEROMETRY 
12 

Book  interferooetry  is  a  Bethod  by  which  excited  state  population  densities  c».n  be 
aeasured.  It  has  been  used  on  lasers  before^^  but,  to  our  knowledge,  this  is  the  first  time 
it  has  been  applied  to  e-beaa  pumped  excimer  lasers.”* 


Figure  8.  Typica.1  amell  aigom.!  gain 
■pectruB  of  XeF  for  aa  •- 
beiia  deposition  of 
*180  kW/ca^  and  a  gas 
mixture  of  60. 4X  Ne/0.6X 
Xe/O.lX  F2  at  3.5  aaagats . 


A  schematic  of  the  hook  interferometer  system  is  shown  in  Fig.  0.  A  short  pulse, 
broadband  dye  laser,  centered  at  *>830  nm,  (LDS  821  dye  in  DMSO  *  lOX  methanol)  feeds  a 
kichelson  interferometer  in  which  the  laser  chamber  is  in  one  branch.  The  output  from  the 
interferometer  is  analyzed  by  a  spectrometer  whose  output  is  detected  by  a  widicon  camera. 
Due  to  anomalous  dispersion  occurring  around  an  atomic  resonance  transition,  *hooks"  form  in 
the  interference  pattern  seen  by  the  camera.  The  separation  of  these  hooks  is  directly 


Schematic  of  book 
interferometer  system  for 
measuring  xenon  excited 
state  densities. 


proportional  to  the  excited  state  population  difference  of  the  transition.  If  the 
transition  oscillator  strength  and  the  active  medium  length  are  known,  the  absolute  value  of 
the  excited  state  population  difference  can  be  determined. 


Figure  6. 


SPECTROMETER 


An  example  of  the  hook  spectra  obtained  during  the  experiment  for  XeCl  is  given  in 
Fig.  10.  Visible  are  the  xenon  transition  lines  at  823.2  nm  (6b [3/2)2  -  8p[3/2]2)  and 
828.0  nm  (68[3/2]j  -  6p[l/2]^).  (The  line  at  826.6  nm  corresponds  to  Xe  (6s’[l/2]j  - 
6p’[l/2]^)  and  was  not  examined  during  this  experiment). 


Figure  10.  Typical  hook  spectra 

obtained  from  the  system 
shown  in  Fig.  0.  The  hook 
separation  is  represented  by 
A. 
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Figur*  11  ahowB  txpicAl  xanon  axcitad  atxta  da&aity  rasulta  for  XaCl .  This  type  of 
information  ia  important  bacauaa  of  tha  rola  axcitad  atata  xanon  playa  in  tha  formation  of 
XaCl .  The  Xa*  danaity  ia  alao  affaetad  by  tha  conaumption  of  tha  halosen  gaa. 


Thia  aama  tachniqua  ia  applicabla  to  othar  axcimer  laaara,  auch  aa  KrF,  whara 
maaaurananta  of  tha  krypton  axcitad  atataa  danaitiaa  ia  poaaibla. 


Ne/Xe/HCI  s  08.2/ 1.5/0.32% 
Total  Prasaura  3000  Torr 


n 


400  600  800 

TIME  (ns) 


'igura  11.  Typical  xanon 
axcitad  atata 
danaity  raaulta 
for  XaCl  (08. 2X 
Na/1.5*  Xe/0.32S 
BCl  at  3000 
Torr) .  The  a-beaa 
dapoaition  ia 
a250  kW/cm®;  a 
normal iaad  trace 
of  tha  a-baam 
currant  pulae  ia 
ahown  at  the 
bottom  of  the 
fifura.  Solid  and 
opan  circlaa 
corraapond  to 
danaity 

diffarancaa  of  the 
823.2  and  828.0  nm 
tranaitiona , 
raapactirely . 


7.  CONCLUSIONS 

Measurementa  of  fundamental  quantitiea,  auch  aa  time-dependent  electron  and  excited  atate 
densities,  has  helped  improve  our  understanding  of  the  complex  kinetics  occurring  in  e-beam 
excited  axcimer  lasers.  It  bas  also  provided  important  data  for  validating  models  of  these 
lasers.  Unlike  typical  data  used  to  anchor  the  models,  thia  data  tests  the  models  on  a 
basic  level. 


I 


The  techniques  described  are  applicable  to  other  typas  of  lasers,  for 
excited  axcimer  and  other  gas  lasers. 


sple  discharge 
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ABSTRACT 


Understanding  of  the  neutral  channel  formation  kinetics  in  excimer  laser 

gas  mixtures  has  been  limited  by  the  lack  of  data  on  the  pertinent  excited 

state  populations  in  these  mixtures.  Presented  are  time-dependent 

measurements  of  the  lower  level  xenon  excited  state  densities  in  electron-beam 

(e-beam)  pumped  XeCl  and  XeF  laser  mixtures  (neon  diluent) .  Measurements  are 

obtained  using  book  interferometry  under  nonlasing  conditions  at  an  average 

3 

excitation  rate  of  R)250  kW/cm  and  e-beam  pulse  lengths  of  0.4  and  1  jts.  The 

♦ 

population  differences,  AN  ,  between  four  different  electronic  transitions 
[three  in  the  Xe* (6s) -Xe** (6p)  manifold,  and  one  in  the  Xe** (6s ’) -Xe*** (6p ’ ) 

manifold]  are  examined  as  a  function  of  halogen  concentration.  For  both  XeCl 

* 

and  XeF  at  high  initial  halogen  concentrations  (>4  Torr) ,  the  AN  densities  of 

the  Xe*(6s)  and  Xe**(6s’)  transitions  are  relatively  constant  during  the 

*  14  -3  .  ,  14 

0.4  jis  e-beam  pulse  [For  AN  (6s-6p)  :  ®<4  x  10  cm  for  XeCl,  and  «1.5  x  10 

-3  .  *  .  . 

cm  for  XeF].  At  lower  initial  halogen  concentrations,  the  AN  densities  of 


6s-6p  and  6s ’-6p’  start  at  the  beginning  of  the  pulse  at  approximately  the 

same  densities  as  the  richer  halogen  mixtures,  but  at  a  certain  point  during 
♦ 

the  pulse,  the  AN  densities  abruptly  increase.  This  increase  can  be  >10 

times  for  very  lean  halogen  mixtures  (1-2  Torr) ,  and  occurs  at  earlier  times 

as  the  initial  halogen  concentration  is  reduced.  From  other  measurements, 

this  increase  appears  related  to  the  depletion  of  the  halogen.  The  observed 
* 

lifetime  of  the  Xe  (6s)  densities  is  <^2  fis  for  the  low  initial  halogen 
concentration  mixtures.  Additional  density  data  for  halogen-free  Ne/Xe  and 
Ar/Xe  mixtures  are  also  presented. 


I.  INTRODUCTION 


XeCl  and  XeF  lasers  have  been  operated  with  electron  beam  (e-beam) 

1  2 

pumping  for  more  than  ten  years  ’  ;  however,  many  questions  remain  regarding 
the  chemical  kinetics  and  dynamical  behavior  of  the  laser  media.  In 
particular  for  XeCl  lasers,  the  formation  of  XeCl*  in  e-beam  pumped  lasers  has 

O  A 

been  the  topic  of  considerable  debate^*  .  In  general,  there  has  been  a  lack 
of  data  on  fundamental  quantities  that  are  directly  related  to  the  media 
kinetics.  This  paper  presents  experimental  results  on  the  lower  level  xenon 
excited  state  densities  (more  precisely,  population  differences)  in  e-beam 
pumped  XeCl  and  XeF  laser  mixtures  (nonlasing).  As  will  be  shown,  the  results 
have  demonstrated  some  unexpected  behaviors  and  revealed  gaps  in  our 
understanding  of  the  neutral  channel  kinetics.  The  measurements  have  also 
provided  useful  data  to  help  validate  computer  models  of  these  excimer 
lasers . 

5 

For  XeCl,  reaction  of  the  first  excitation  state  of  xenon  (5p  6s)  with 

HCl(i/=0)  is  nearly  thermoneutral  (i.e.  the  branching  ratio  is  nearly  zero) 

with  respect  to  XeCl*  formation.  Either  higher  xenon  excited  states  or 

vibrationally  hot  HCl  are  necessary  to  efficiently  form  XeCl*  through  the 

harpoon  reaction^.  Recently,  quenching  rate  constants  and  branching  fractions 

for  XeCl*  formation  generated  by  reaction  of  various  xenon  excited  states  with 

0 

halogen  donors  have  been  measured  using  a  tunable  UY  synchrotron  source”  and 
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by  laser  induced  fluorescence  .  These  results  showed  that  HCl(iM))  gives  high 
XeCl*  formation  yield  only  for  the  6p,  5d,  7s,  and  presumably  higher 


excitation  states  of  xenon.  Enhancement  of  XeCl*  formation  yield  from  Xe*(6s) 
with  vibrationally  excited  HCl  has  also  been  observed  .  It  has  been  assumed 
in  the  past  that  the  population  in  higher  xenon  excited  states  is  minimal, 
that  the  role  of  excited  state  kinetics  is  less  important  in  XeCl  than  in 
other  rare-gas  halide  lasers,  such  as  KrF*,  and  that  XeCl*  is  predominantly 
formed  through  Xe*-Cl~  recombination®*^^. 

In  Ne/Xe  mixtures  (without  HCl),  neon  ions  and  excited  states,  produced 
by  the  primary  electrons  of  the  e-beam  and  by  high  energy  secondary  electrons, 
produce  xenon  ions  through  three-body  charge  transfer  reactions  and  Penning 
ionization.  Xenon  excited  states  are  then  produced  through  recombination 
reactions  of  xenon-based  dimer  ions  (Xeg,  NeXeg) .  The  fraction  of  xenon 
excited  states  formed  directly  by  electron  collisions  is  negligible  due  to  the 
low  density  of  xenon  compared  to  neon  gas  in  these  laser  mixtures  (typically 
Xe  is  0.5-1. 5%  of  the  total  mixture).  With  HCl  in  the  mixture,  as  long  as  the 
Cl  concentration  is  high,  the  Xe*-Cl  recombination  reaction  is  fast  and  the 
xenon  excited  state  formation  channels  are  effectively  intercepted  preventing 
the  formation  of  xenon  excited  states.  Hence,  measurement  of  the  time- 
dependent  xenon  excited  state  densities  provides  an  important  check  of  our 
understanding  of  not  only  the  neutral  channel  kinetics,  but  also  the  ion 
channel  kinetics.  For  example,  if  the  measured  xenon  excited  state  densities 
are  reasonably  low,  this  is  indirect  evidence  of  high  Cl~  ions  density; 
however,  the  dominant  source  of  Cl~  is  dissociative  attachment  from 
vibrationally  excited  HCl^^  and  there  are  experimental  results  indicating  that 
the  density  of  vibrationally  excited  HCl  is  low  ~  .  In  conjunction  with  our 
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previous  measurements  on  the  secondary  electron  density  ’ 
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depletion  ,  the  xenon  excited  state  density  measurements  related  herein  can 
significantly  help  improve  our  understanding  of  the  complex  raire-gas  halogen 
kinetics. 

In  this  paper,  the  Racah  notation  will  be  used  to  specify  a  given  state 
of  xenon.  Specific  xenon  excited  states  are  divided  to  three  groups  according 
to  their  energy  levels.  The  symbol  Xe*,  Xe**,  and  Xe***  denotes  Xe(6s), 
Xe(6s*,  6p,  and  5d) ,  and  Xe (higher  than  5d;  states,  respectively. 


II.  EXPERIMENTAL  TECHNIQUE 

1  ft 

Xenon  excited  state  densities  are  measured  using  book  interferometry 

17 

This  technique  hais  been  used  before  in  gas  lasers  ,  but  this  is  the  first 

application  in  e-beam  pumped  excimer  lasers.  The  experimental  apparatus  is 

described  in  detail  elsewhere  ;  a  schematic  is  shown  in  Fig.  1.  The  total 

optical  path  of  the  excited  laser  medium  is  70  cm,  and  the  e-beam  excitation 
3 

rate  is  «250  kW/cm  ,  which  is  higher  thaui  during  the  earlier  electron  density 
13  14 

measurements  ’ 

The  basic  principle  of  hook  interferometry  relies  on  the  anomalous 
dispersion  occurring  around  an  atomic  resonance  transition;  the  strength  of 
that  dispersion  is  directly  proportional  to  the  population  difference  of  that 
transition.  By  placing  the  laser  medium  in  a  Michelson  interferometer  and 
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i 


analyzing  the  interference  pattern  with  a  spectrometer,  the  anomalous 


dispersion  manifests  itself  in  the  form  of  'hooks*  in  the  interference  pattern 


at  the  output  of  the  spectrometer. 


Near  an  isolated  transition  at  a  wavelength  (where  i  and  j  are  the 


lower  and  upper  levels,  respectively),  but  outside  the  core  of  the  spectral 


line,  the  refractive  index  n  is  given  by 


r  f..X..N. r  g.N. 

n(X)-l  =  1  -  ^ 

4ir  X-X..  g.N. 

ij  I  1 


where  r^  is  the  classical  electron  radius,  fj^j  is  the  oscillator  strength  of 


the  transition,  N*  .  are  the  population  densities,  and  g.  .  are  the 

j  1 ( J 


statistical  weights  of  each  level.  The  wavelength  separation.  A,  between  the 


hooks  in  the  interference  pattern  is  directly  proportional  to  the  square  root 


of  the  population  density  difference  between  the  two  states  and  is  given  by: 


AN*  =  N.  1 


r  x^.ef. . 
o  ij  ij 


where  Z  is  the  length  of  the  medium  and  K  is  a  constant  that  can  be  determined 


from  an  undistorted  interference  pattern.  If  the  oscillator  strength  for  a 


particular  transition  is  known,  then  the  population  difference,  AN  ,  can  be 


determined  from  the  hook  spectra. 


An  energy  level  diagram  showing  the  xenon  transitions  examined  during 
the  experiment  is  given  in  Fig.  2.  Four  particularly  strong  transitions  are 


examined:  1)  823.2  nm  (6s[3/2]^  -  6p[3/2]2);  2)  828.0  nm  (6s[3/2]J  - 
6p[l/2]^j):  3)  834.7  nm  (6s’ [1/2]  J  -  6p’[3/2]2);  and  4)  840.0  nm  (6s[3/2]°  - 
6p[3/2]j^).  yeasured  oscillator  strengths  for  these  lines  are  0.23,  0.12, 

0.37,  and  0.012,  respectively  .  In  this  paper,  the  measured  population 
density  differences  corresponding  to  each  of  the  transitions  are  referred  to 
as  AN* (823. 2  nm) ,  AN* (828.0  nm) ,  AN* (834.7  nm) ,  and  AN* (840.9  nm) , 
respectively. 

3 

The  value  of  K  in  Eq.  (2)  is  typically  «1.5  x  10  during  the 

experiments.  Thus,  on  our  system  a  hook  separation  of  1  nm  at  823.2  nm 
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corresponds  to  a  AN  density  of  <^7  x  10  cm  .  Because  the  823.2  nm  and 

828.0  am  lines  are  close  to  one  another  and  they  have  relatively  large 

oscillator  strengths,  the  hook  patterns  about  these  lines  can  overlap  each 

other  at  high  xenon  excited  state  densities.  This  can  be  avoided  in  principle 

by  reducing  the  active  length  of  the  medituB  (i.e.  £.)  ;  however,  this  was  not 

convenient  to  do  during  the  experiment.  Bence,  most  of  the  xenon  excited 

15  -3 

state  densities  greater  than  10  cm  are  measured  during  the  experiment  by 
examining  the  840.9  nm  transition,  because  of  its  smaller  oscillator  strength, 
rather  than  the  823.2  and  828.0  nm  lines. 

Some  of  the  transitions  examined  have  a  spectral  line  nearby  that  can 
interfere  with  the  hooks  formed  on  the  side  of  the  transition  nearest  the 
spectral  line  .  Two  relatively  strong  transitions  are  observed  at  826.6  and 
822.1  nm  near  the  828.0  and  823.2  nm  transitions.  The  transition  at  826.6  nm 
corresponds  to  Xe  (6s’ [1/2]°  -  6p’[l/2]|).  The  transition  at  822.1  nm  has  not 


been  positively  identified,  but  it  is  probably  from  Ne^,  because  it  is  also 
observed  in  pure  neon  mixtures.  These  lines  disturb  the  hook  formation  on  the 
short  wavelength  side  of  the  828.0  and  823.2  nm  transitions,  while  the  effects 
on  the  long  wavelength  side  of  the  transitions  are  negligible. 

For  an  isolated  transition,  the  hook  separation  is  symmetric  about  the 

transition  wavelength.  The  quadratic  Stark  effect  can  displace  the  hooks 

20 

towards  long  wavelengths  ,  but  because  the  electron  density  is  relatively 
13,14  15  -3 

low  (<oiO  cm  ),  this  effect  cam  be  ignored.  Therefore,  during  analysis 

of  the  data,  it  is  assumed  the  total  hook  separation  can  be  calculated  by 
doubling  the  distance  between  the  hook  on  the  long  wavelength  side  and  line 
center.  The  experimental  uncertainty  of  the  data  presented  in  this  paper  is 
«*10fS. 


Because  hook  interferometry  gives  the  density  difference  between  two 
states,  to  obtain  the  absolute  density  of  the  lower  state  requires  the  upper 
state  population  to  be  either  negligibly  small  or  estimable  through  other 
means.  As  explained  in  the  following  section,  we  believe  the  upper  state 
populations  for  our  experimental  conditions  are  signif iczoit;  however, 
estimating  the  populations,  even  by  comparing  the  data  gathered  for  different 
transitions,  is  difficult.  Therefore,  in  order  to  avoid  biasing  the  data  with 
an  assumption  of  the  upper  state  densities,  all  the  data  presented  in  this 
paper  represent  the  population  difference,  AN  ,  between  the  upper  and  lower 
state  of  the  transition  examined. 
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III.  RESULTS  AND  DISCUSSION 


Before  presenting  the  experimental  results,  it  is  helpful  to  review  the 

important  processes  affecting  the  xenon  excited  state  population 

distributions.  In  high  pressure  gas  mixtures  excited  by  an  e-beam  at  a 

relatively  high  excitation  rate,  the  population  distributions  in  excited  state 

manifolds  are  primarily  defined  by  seven  different  processes:  1)  The  initial 

21 

distribution  as  a  result  of  formation  by  dimer  ion  recombinations  ;  2)  two- 
body  collisional  redistributions  with  rare  gases;  3)  three-body  collisional 
deactivations  (e.g.  dimer  molecule  formations);  4)  radiative  deactivations;  5) 
excitations  and  de-excitations  with  secondary  electrons;  6)  quenching  with  the 
halogen  gas;  and  7)  ionization  of  high  level  excited  states. 

If  the  collisional  mixing  by  either  neutral  species  or  electrons  is 

fast,  the  initial  formation  distribution  will  be  quickly  reduced  to  a  steady- 

state  distribution.  The  nature  of  this  distribution  will  depend  on  the 

relative  rates  of  neutral  and  electron  energy  transfer  collisions.  The 

collisional  deactivation  rate  constants  measured  by  Inoue  et  al.  for  Ne  and 

Ar  show  that  two-body  deactivations  redistribute  xenon  densities  between 

excited  states  separated  by  <5000  cm  ^  with  rate  constants  ranging  from  0.2  to 

36  X  10  cm  s  Radiative  lifetimes  measured  by  Inoue  et  al.  range  from 

30  to  150  ns  for  xenon  states  located  above  6p[l/2]j^.  Since  the  secondary 

electron  number  density  is  typically  lO^^-lO^®  cm~^  in  the  excited  gas  mixture 
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and  their  energy  is  »l-2  eV  ,  it  is  reasonable  to  expect  that  stepwise 
excitation  in  the  xenon  excited  state  manifold  will  effectively  compete  with 


two-body  collisional  deactivation.  If  the  distribution  in  the  xenon  excited 

state  manifold  is  characterized  by  a  secondary  electron  temperature  of  1.5  eV, 
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which  is  a  typical  average  electron  temperature  predicted  by  our  model  ’ 
for  our  experimental  conditions,  the  upper  state  density  of  the  probed 
transitions  would  be  23%,  12%,  and  36%  of  the  lower  state  density  for  the 
823.2,  828.0,  and  840.0  nm  transitions,  respectively.  This  represents  an 
upper  limit  for  the  higher  state  populations  since  other  heating  sources  do 
not  exist.  Moderation  of  this  distribution  will  be  governed  by  interaction 
with  the  ^300  K  gas.  Moreover,  ionization  processes  of  higher  Xe  excited 
states  tend  to  decrease  the  population  of  upper  states. 

At  high  pressures,  three-body  relaxations  are  also  important.  Xe***  and 
Xe**  form  unstable  Xe^**  aad  Xe^*  dimers,  respectively,  that  caui  eventually 
predissociate  into  lower  xenon  excited  state  states.  If  the  contribution  of 
these  three-body  collisional  deactivations  is  significant,  the  characteristic 
temperature  defining  the  distribution  between  Xe*  and  Xe**  can  be  lower  than 
the  electron  temperature.  Consequently,  the  fraction  of  upper  state  densities 
can  be  lower  than  those  estimated  in  the  previous  paragraph.  The  issue  of  the 
xenon  manifold  temperature  is  discussed  in  more  detail  later. 

Figures  3  and  4  show  typical  experimental  results  obtained  at  the  828.0, 
823.2,  and  834.7  nm  transitions,  for  XeCl  laser  mixtures  with  an  initial  HCl 
concentration  of  0.32%  and  0.16%,  respectively.  The  XeCl  laser  mixture  is 
1.5%  Xe  in  Ne  diluent  at  a  total  gas  pressure  of  3000  Torr.  A  normalized  e- 
beam  current  pulse  shape  [400  ns  full  width  at  half  maximum  (FWHM)]  is  shown 
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at  the  bottom  of  the  figures.  As  can  be  seen  in  Figs.  3  2uid  4,  the  density 
♦  * 

differences  of  AN  (823.2  nm)  and  AN  (828.0  nm)  are  relatively  constant  during 
the  e-beam  pumping,  and  decay  rapidly  at  the  termination  of  the  e-beam  pulse. 

As  noted  before,  excited  xenon  is  produced  primarily  by  recombination  of 
dimer  ions  with  secondary  electrons  rather  thaui  direct  excitations  of  ground 
state  Xe  by  high  energy  electrons.  At  high  HCl  concentrations,  secondary 
electrons  are  efficiently  removed  by  dissociative  attachment  reactions  with 
vibrationally  excited  HCl.  The  Cl~  ions  produced  reduce  the  dimer  ion 
populations  through  ion-ion  recombination  reactions.  Therefore,  as  long  as 
dissociative  attachment  dominates  the  electron  loss  process,  the  Xe*  density 
will  decrease  with  increasing  initial  HCl  concentration.  In  fact,  as  seen  in 
Figs.  3  and  4,  experimental  results  for  0.3251  HCl  of  AN  (823.2  nm)  are  tt475{  of 
the  0.165(  HCl  mixture.  These  Xe  densities  agree  with  model  predictions 
which  assume  high  dissociative  attachment  rate  constants  for  vibrationally 
excited  HCl  as  reported  by  Allan  and  Wong^^. 

* 

Figures  5-8  show  time-dependent  AN  densities  for  lower  initial  HCl 

concentrations  from  0.12  to  0.045S.  Two  unexpected  xenon  excited  state  density 

behaviors  are  observed.  First,  although  higher  Xe*  densities  are  expected  at 

* 

the  beginning  of  the  pulse  for  lower  HCl  concentrations  and,  in  fact,  the  AN 
measured  during  the  first  200  ns  does  increase  when  the  initial  HCl 

concentration  is  reduced  to  0.099,  the  AN  densities  actually  tend  to  decrease 

* 

as  the  HCl  concentration  is  reduced  to  0.08  and  0.049.  The  AN  (834.7  nm) 
density  of  0.089  HCl  is  also  smaller  than  that  of  the  0.169  HCl  mixture  during 


the  e-beam  pulse.  Second,  at  low  HCl  concentrations,  the  AN  densities 

dramatically  increase  at  a  well  defined  point  during  the  e-beam  pulse.  As  the 

initial  HCl  concentration  is  reduced,  the  AN  density  begins  to  increase  at 

earlier  times.  This  increase  is  so  rapid  that  it  appears  to  indicate  some 

type  of  threshold  effect.  For  the  particular  energy  loading  and  pulse 

duration  used  in  Figs.  3-8,  this  increase  occurs  at  the  very  end  of  the  e-beam 

* 

pulse  for  the  0.12%  HCl  mixture,  and  the  AN  density  decays  rapidly  after 

* 

termination  of  the  e-beam.  However,  for  HCl  concentrations  <0.12%,  the  AN 
density  does  not  decay  immediately  after  termination  of  the  e-beam  pulse,  but 
instead  the  decay  time  becomes  progressively  longer  with  decreasing  initial 
HCl  concentration.  This  will  be  discussed  in  more  detail  later. 

* 

Also  shown  in  Fig.  7  are  results  of  measuring  the  AN  density  of  neon 
excited  state  for  the  (3p[5/2]2  -  3d[7/2]^)  transition  at  837.8  nm.  It  is 
interesting  that  the  neon  excited  state  density  does  not  increase  abruptly 
like  the  xenon  densities,  but  rather  follows  the  e-beam  pumping.  The  abrupt 
increase  of  the  xenon  excited  state  density  appears  to  be  caused  by 
destruction  of  the  reaction  balance  between  ion-ion  and  electron-dimer  ion 
recombination  processes,  as  a  result  of  HCl  consumption.  Since  neon  excited 
states  are  mainly  produced  by  high  energy  electron  excitations,  they  are  not 
affected  by  the  HCl  consumption.  During  the  electron  density 
measurements^  ’  ,  the  electron  densities  also  tended  to  increase  abruptly  at 

a  certain  point  during  the  e-beam  pumping.  This  point  also  tends  to  occur  at 
earlier  times  as  the  initial  HCl  concentration  is  lowered.  The  consumption  of 
HCl  during  the  pump  pulse  increases  the  secondary  electrons  and  reduces  Cl 
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ions  resulting  in  an  increase  in  electron-dimer  ion  recombinations. 

23  . 

The  model  predicts  these  abrupt  increases  only  when  additional 
dissociation  of  HCl  to  neutrals  is  included  at  effective  rates  larger  than 
commonly  accepted.  In  other  words,  the  HCl  burn-up  rate  appears  to  be 
underestimated,  especially  at  low  initial  HCl  concentrations.  Additional  HCl 
destruction  channels  at  low  HCl  concentrations  may  still  not  be  identified  yet 
and  will  be  discussed  later. 

At  low  HCl  concentrations,  the  AN  (840.9  nm)  density  decays  slowly  after 
termination  of  the  e-beam  pulse.  In  high  pressure  pure  rare  gas  mixtures 
excited  by  an  e-beam,  electrons  are  controlled  by  recombinations  with  dimer 
ions.  The  rare  gas  excited  states  formed  by  recombination  processes  produce 
dimer  molecules  through  three-body  collisions,  and  the  inner  energy  of  the 
total  xenon  system  decays  mainly  through  radiative  deactivations  of  dimer 
molecules.  Therefore,  as  long  as  the  three-body  dimer  molecule  formation  is 
fast,  rare  gas  excited  states  cannot  survive  for  very  long  after  the  e-beam 
termination  even  after  the  HCl  is  completely  consumed.  During  the  secondary 
electron  density  measurements,  a  long  electron  density  lifetime  was  also 
observed  at  low  HCl  concentrations  '  .  To  conserve  charge  balance,  the 

secondary  electron  density  must  be  balanced  with  some  positive  ions. 

Therefore,  the  long  xenon  excited  state  lifetime  may  be  caused  by  relatively 
long  deactivation  processes,  and  by  some  feeding  processes  such  as 
recombinations  of  dimer  ions  even  after  the  e-beam  excitation.  Energy 
transfer  processes  from  excited  Cl  atoms  are  another  possible  path  to  form 


xenon  excited  states.  Comparing  Figs.  6  and  7,  it  is  interesting  to  note  how 
sensitive  this  lifetime  is  to  small  changes  in  the  initial  halogen 
concentration . 

To  eliminate  the  effects  of  the  halogen,  Ne/Xe  and  Ar/Xe  mixtures  were 
also  examined  using  the  hook  interferometry  system.  The  results  for  the 
binary  mixtures  are  shown  in  Fig.  9.  The  Xe*  lifetime  is  surprisingly  long, 
and  its  decay  behavior  is  not  purely  exponential. 

To  help  interpret  the  binary  mixtures  results,  Fig.  10  depicts  the  key 
kinetics  pathways  in  an  e-beam  excited  Ne/Xe  gas  mixture.  The  major  dimer  ion 
species  in  Ne/Xe  mixtures  are  believed  to  be  NeXe*  and  Xe^.^'^^  In  neon  gas 
containing  a  trace  of  xenon,  Xej  is  formed  mainly  from  NeXe*  rather  than 
directly  from  Xe*.  Similarly,  the  major  dimer  molecules  in  Ne/Xe  mixtures  are 
thought  to  be  NeXe*  and  Xe^.  However,  the  heteronuclear  species  NeXe"^  and 
NeXe*  are  bound  very  weakly  (<0.05  eV) ,  so  that  the  reverse  reaction  to 
produce  Xe  and  Xe  from  these  heteronuclear  species  should  be  very  fast 
If  these  reverse  reactions  are  fast  enough  that  the  effective  densities  of 
heteronuclear  species  are  negligibly  small,  and  if  the  formation  reactions  of 
molecular  species  of  Xe^  and  Xe^  with  neon  as  a  third- body  species  are 
significantly  slower  than  expected,  viz.  the  reactions: 


then  the  total  rates  of  for  nation  of  Xe^  aod  Xe^  will  be  much  less  than  those 

for  tne  corresponding  rare  gas  dimer  species  in  pure  xenon.  In  fact,  the 

rates  used  in  the  model  are  estimated  or  extrapolated  from  results  for 

collision  species  other  than  neon  .  Consequently,  the  densities  of  Xe- 

baksed  molecular  species  are  much  smaller  than  would  be  expected  considering 

the  total  pressure.  If  the  preceding  assumptions  are  correct,  then  the 

dominant  ion  in  the  Ne/Xe  mixture  is  Xe^,  and  the  electron  loss  mechanism  by 

electron-dimer  recombination  is  slow,  with  the  result  that  the  electron 

density  in  the  Ne/Xe  mixture  remains  high  after  the  termination  of  e-beam 

pumping.  Penning  ionization  processes  can  also  help  feed  the  electron  density 

population.  This  is  one  noteworthy  difference  between  a  pure  rare  gas  and  a 

rare  gas  containing  a  trace  of  xenon.  Therefore,  even  after  termination  of  e- 

beam  pumping,  the  recombination  process  producing  xenon  excited  states 

continues  at  a  low  rate,  some  of  the  excited  species  are  ionized  back  to  Xe**^ 

through  Penning  ionization,  and  the  total  xenon  excited  state  system  decays 

through  the  relatively  slow  Xe^  formation  process.  If  the  electron 

temperature  remains  high  enough,  the  produced  Xe^  is  excited  to  unstable  Xe^*^ 
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states,  which  again  produce  Xe  . 

Figure  11  shows  the  fluorescence  time  history  of  the  Ne/Xe  mixture 
obtained  with  a  VUV  photodiode  (ITT  F4115,  120  nm<X<320  nm) .  The  major 
contribution  is  dimer  emission  of  Xe^.  Surprisingly,  the  fluorescence 
continues  to  increase  after  termination  of  the  e-beam  pulse  and  peaks  at 
^1.25  /is  after  the  pulse.  This  fluorescence  behavior  of  Xe^  is  caused  by  a 
decrease  of  Xe^*  formation  by  secondary  electrons  due  to  the  drop  of  the 
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electron  temperature  after  the  e-beam  ends^  . 


In  Fig.  0,  AN  (823.2  nm,  and  828.0  nm)  measured  for  an  Ar/Xe  mixture  are 

also  shown.  The  decay  time  of  the  xenon  excited  state  system  is  much  faster 

than  in  the  Ne/Xe  mixture.  This  is  because  the  binding  energies  of  ArXe^  and 

ArXe*  are  relatively  large  compared  to  those  of  NeXe*  and  NeXe*,  respectively, 

and  a  larger  amount  of  dimer  molecules  exist  in  Ar-based  mixtures'^  .  Xe^  and 

Xe2  formation  reactions  with  argon  as  a  third-body  species  are  also  faster 
26~28 

than  with  neon^  .  Therefore,  the  recombination  and  deactivation  processes 
for  xenon  excited  state  are  faster  than  those  in  Ne-based  mixtures. 


Figures  12  and  13  shows  time-dependent  AN  densities  obtained  under 
extended  (1  /(s)  e-beam  excitation  for  an  initial  BCl  concentration  of  0.32  and 
0.16%,  respectively.  Sidelight  fluorescence  of  XeCl*(B-X)  is  shown  at  the 

bottom  of  each  figure  with  the  normalized  e-beam  current  pulse.  For  0.32% 

* 

HCl,  the  AN  density  does  not  stay  constant,  but  gradually  increases  towards 

the  end  of  the  e-beam  pulse.  This  gradual  increase  may  be  caused  by  an 

increase  in  the  gas  temperature  during  the  e-beam  pumping  which  lea^ls  to  a 
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decrease  in  the  ion-ion  recombination  rates  .  For  0.16%  HCl,  the  AN  density 
increases  abruptly  at  ^600  ns  into  the  e-beam  pulse.  The  sidelight 

fluorescence  for  the  0.16%  HCl  mixture  begins  to  decrease  at  approximately  the 

*  * 
same  time  as  the  increase  of  the  AN  density.  Since  the  increase  of  AN 

density  also  suggests  a  corresponding  increase  of  the  Xe'**  and  Xe^'** 

densities,  photoionizations  of  these  higher  excited  states  should  increase  the 

absorption  losses  in  a  laser  and  tend  to  reduce  the  laser  output.  This  may  be 

part  of  the  reason  for  the  premature  laser  pulse  termination  observed  in  some 


experiments  . 


Measurements  also  shov  that  the  Xe*  density  increases  with  increasing 

initial  xenon  concentration.  With  0.16!{  BCl  and  at  150  ns  into  the  e-beam 
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pulse,  AN  (823.2  nm)  densities  of  Rf2.0  x  10  and  k<4.5  x  10  cm  are 

obtained  for  0.5%  and  3.0^  Xe,  respectively. 

It  is  interesting  to  examine  whether  the  abrupt  increase  of  xenon 
excited  state  density  is  a  common  behavior  when  using  different  halogen 
species.  Results  of  AN  density  measurements  performed  for  XeF  laser  mixtures 
are  shown  in  Figs.  14-16.  The  XeF  laser  mixture  is  0.5%  Xe  in  Ne  diluent  to  a 
total  gas  pressure  of  2888  Torr. 

♦ 

One  must  be  careful  when  comparing  the  AN  results  for  Fg  with  HCl 

because  of  the  different  dissociative  attachment  rate  constauits  between  ground 

state  HCl  and  F2.  The  rate  constants,  however,  actually  become  comparable  to 

one  another  once  the  HCl  is  vibrationally  excited.  Although  the  0.1%  F2 

mixture  (Fig.  15)  and  the  0.09%  HCl  mixture  (Fig.  6)  start  out  with  similar 

* 

amounts  of  halogen  and  the  rapid  increase  in  AN  density  occurs  approximately 
at  the  same  time,  the  lifetime  after  e-beam  termination  of  the  xenon  excited 
states  is  much  longer  for  the  0.1%  F2  case.  Based  on  the  earlier  discussion 
regarding  the  long  lifetimes  observed  in  the  binary  rare  gas  mixtures,  this 
seems  to  imply  that  more  of  the  F2  has  been  consumed  than  in  the  HCl  case. 

The  estimated  effective  lifetime  of  the  xenon  AN  excited  state  for  the  0.05% 
F2  mixture  is  ^2  jis,  which  agrees  well  with  the  value  for  the  XeCl  laser 


mixture . 


* 

From  the  AN  density  results  for  ECl  and  F2  halogens,  it  is  confirmed 
that  regardless  of  which  halogen  gas  is  used,  the  halogen  burn-up  causes  an 
abrupt  change  in  the  energy  flow  of  the  dimer  ion  recombination  channel  and 
increases  the  xenon  excited  state  density.  In  the  Ne-based  mixture  with  a 
trace  of  xenon,  the  main  ion  species  is  Xe^,  and  recombination  processes  are 
relatively  slow.  Therefore,  the  decay  rates  of  the  secondary  electron  and 
xenon  excited  state  densities  are  very  long. 

The  xenon  excited  state  density  results,  besides  demonstrating  the 
interdependence  on  the  halogen  concentration,  also  provide  an  opportunity  to 
examine  part  of  the  xenon  excited  state  distribution  and  therefore  estimate 
indirectly  the  manifold  temperature.  However,  to  perform  such  an  analysis 
requires  some  assumption  regarding  the  form  of  the  distribution 
(e.g.  Boltzmann  versus  non-Boltzmann).  Because  of  this  and  because  the 
experiment  only  examined  a  limited  portion  of  the  xenon  manifold  and  did  not 
measure  absolute  level  densities,  the  following  temperature  analysis  is 
intended  to  illustrate  only  one  possible  interpretation  of  the  data.  It 
should  be  emphasized  that  the  experiment  was  not  designed  to  measure  the 
manifold  temperature. 

The  analysis  is  further  complicated  by  the  fact  that  at  least  part  of 
the  xenon  manifold  must  be  in  nonequilibrium  as  evidenced  by  the  existence  of 
infrared  xenon  lasers.  Lasing  occurs  between  the  5d  and  6p  xenon  excited 
states  when  mixtures  of  Ar/Xe  or  Ne/Xe  at  >1  atm  are  excited  by  e-beams  or  e- 


.  If,  however,  we  assume  the  levels  examined  in 
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beam  sustained  discharges 
this  experiment  are  roughly  in  equilibrium,  then  one  approach  for  analyzing 
the  data  is  to  assume  a  simple  model  where  the  xenon  excited  state  manifold  is 
characterized  by  a  single  temperature.  It  is  possible  to  estimate  such  a 
temperature  by  assuming  a  Boltzmann  distribution  for  the  entire  xenon  excited 
state  manifold  between  the  6s  and  6p'  levels  (see  Fig.  2).  For  example, 

I 

*  4t  ♦ 

comparing  the  AN  (834.7  nm)  data  to  the  AN  (823.2  nm)  or  AN  (828.0  nm)  results 

given  in  Fig.  4  for  0.165(  HCl  would  indicate  that  the  excitation  temperature 

* 

in  the  Xe  manifold  is  ^0.8  eV.  This  is  comparable  to  the  predicted  electron 
temperature  of  1.5  eV  mentioned  earlier  and  supports  the  hypothesis  that 
electron  mixing  is  significant  in  the  lower  excited  state  manifolds  for  this 
experiment. 

I 

As  another  example  of  the  application  of  this  simple  analysis  to  the 
data,  the  calculated  excitation  temperature  from  the  AN  (840.0  nm)  and 
AN* (834.7  nm)  data  for  a  0.08%  HCl  mixture  is  ^0.0  eV  during  the  e-beam 
pumping.  The  excitation  temperature  decreases  after  the  e-beam  termination  to 
«0.3-0.4  eV.  For  a  Ne/Xe  mixture,  the  excitation  temperature  decreases  from 
0.5  eV  at  400  ns  into  the  e-beam  pulse  to  0.4  eV  at  1.4  /is  after  the  e-beam 
pulse.  Comparing  both  mixtures,  the  halogen  gas  tends  to  increase  the 
excitation  temperature,  but  the  excitation  temperatures  without  any  initial 
halogen  or  after  the  halogen  has  been  depleted  are  very  similar.  This  trend 
is  expected  if  the  excitation  temperature  is  a  positive  function  of  average 
electron  temperature.  With  the  halogen  gas,  low  energy  electrons  are  removed 
through  dissociative  attachment  reactions  and  the  average  electron  temperature 


23 

is  raised;  indeed,  this  is  predicted  by  the  model  . 

If  indeed  the  excitation  temperature  of  Xe*^  state  is  comparable  to  the 
estimated  secondary  electron  temperature,  then  considering  the  large 
degeneracies  of  the  high  excitation  states,  the  total  xenon  excited  state 
densities  above  the  first  excitation  state  (6s)  would  appear  to  be 
significant.  For  example,  assuming  the  entire  xenon  excited  state  manifold  is 
governed  by  a  Boltzmann  distribution,  one  would  predict,  for  a  1  eV  excitation 
temperature,  an  effective  degeneracy  for  Xe**  of  order  50  aind  for  Xe***  a 
value  30  times  greater  than  that.  Relative  to  Xe^,  this  implies  the  Xe'*'*  and 
Xe**'*  densities  are  considerably  larger  than  the  Xe*  density.  However,  due  to 
the  high  ionization  rates  of  the  high  level  xenon  states,  it  is  more  likely 
that  the  excitation  temperature  of  Xe***  is  smaller  than  that  of  Xe**. 
Nonetheless,  even  if  the  densities  of  the  high  level  xenon  states  are  only 
comparable  to  the  Xe*  density,  the  impact  on  the  kinetics  is  significant. 

Xe  can  produce  XeCl  through  reactions  with  ground  state  HCl  at  a  branching 
ratio  of  nearly  one^.  Assuming  the  total  Xe*^  and  Xe***  densities  are  the 
same  as  the  Xe*  density,  and  using  a  Xe*  density  of  5  x  10^^  cm~^  obtained 
from  our  experiment  for  a  0.16%  BCl  mixture,  an  estimated  XeCl*  formation 
efficiency  through  the  harpoon  reactions  with  higher  excited  states  of  xenon 
is  7%,  which  is  approximately  half  of  the  formation  efficiency  estimated  from 
the  small-signal  gain  of  a  XeCl  laser^^.  This  XeCl*  formation  is  also  an 
additional  HCl  burn-up  channel.  Moreover,  since  the  higher  excited  states  of 
xenon  are  easily  ionized  by  electrons  or  photons,  they  are  a  source  of  extra 
secondary  electrons  which  tend  to  accelerate  the  HCl  burn-up. 


The  preceding  analysis  relies  on  the  distribution  being  approximately 
Boltzmann.  If  in  reality  it  is  significantly  non-Boltzmann,  then  it  is  not 
possible,  without  additional  measurements  of  the  higher  xenon  excited  state 
level  densities,  to  infer  the  manifold  temperature  from  our  data. 

IV.  CONCLUSIONS 

The  time-dependent  excited  state  density  measurements  of  xenon  in  e-beam 
pumped  Ne/Xe/HCl  and  Ne/Xe/F2  laser  mixtures  as  a  function  of  halogen 
concentration  have  yielded  important  new  insights  into  the  neutral  and  ionic 
kinetics  occurring  in  these  plasmas. 

For  both  laser  mixtures,  as  long  as  the  halogen  density  is  high, 
secondary  electrons  and  xenon  ions  are  controlled  by  dissociative  attachment 
reactions  with  the  halogen  and  ion-ion  recombination  reactions,  respectively. 
In  addition,  the  amount  of  electron-dimer  ion  recombination  is  relatively 
small,  which  results  in  low  xenon  excited  state  densities.  When  a  certain 
amount  of  halogen  is  consumed,  the  energy  flow  channel  is  very  rapidly 
switched  to  the  dimer  recombination  channel  and  the  excited  state  density  of 
xenon  abruptly  increases.  This  increase  of  the  xenon  excited  state  population 
is  coincident  with  a  similar  increase  of  the  secondary  electrons.  Although 
the  kinetics  model  can  qualitatively  predict  these  trends,  the  HCl  depletion 
is,  thus  far,  underpredicted  especially  for  low  initial  HCl  concentrations.  A 
population  increase  of  the  higher  level  xenon  excited  states  is  also  observed 


as  the  halogen  is  burned  up;  this  presumably  results  in  an  increase  of 
photoionization  at  the  laser  wavelength.  This  may  contribute  to  the  premature 
termination  of  the  XeCl  laser  output  under  high  energy  loading  conditions. 

In  binary  mixtures  of  neon  with  a  trace  of  xenon,  the  major  positive  ion 
species  is  xenon  ion,  and  the  electron  loss  mechanism  by  dissociative 
recombination  is  slow,  with  the  result  that  the  electron  density  in  the  Ne/Xe 
mixture  remains  high  after  the  termination  of  e-beam  pumping.  Penning 
ionization  processes  with  Xe***  also  help  maintain  the  electron  density.  The 
densities  of  xenon-based  molecules  such  as  NeXe^,  Xe^,  NeXe^,  and  Xe^,  are  low 
because  of  the  relatively  small  binding  energy  of  the  heteronuclear  molecules, 
and  the  low  rate  constants  for  three  body  collisional  reactions  with  neon  as  a 
third-body  species.  Therefore,  even  after  the  termination  of  e-beam  pulse, 
the  recombination  processes  producing  xenon  excited  states  still  continue  at  a 
slow  rate,  and  the  total  xenon  excited  state  system  decays  through  the 
relatively  slow  Xe^  formation  process.  A  xenon  excited  state  Xe*(6s)  lifetime 
of  /is  is  observed  for  the  Ne/Xe  mixture  and  also  for  mixtures  with  lean 
halogen  concentrations.  This  is  a  noteworthy  difference  from  pure  rare 
gases. 

Assuming  a  Boltzmann  distribution  for  the  xenon  manifold,  the  excitation 
temperature  during  e-beam  pumping  estimated  from  the  data  for  the  lower  Xe 
manifolds  is  *•!  eV,  which  is  very  close  to  the  predicted  electron  temperature 
of  1.5  eV  by  the  model.  If  most  of  the  distribution  of  the  xenon  excitation 
manifolds  are  characterized  by  this  high  temperature,  a  significant  population 


of  xenon  excited  states  exists  above  the  first  excitation  level  of  Xe*(6s)  due 
to  the  large  degeneracy  of  the  higher  excitation  levels.  The  higher  state 
xenon  excited  states  can  produce  XeCl'*  with  a  high  formation  yield,  and  are 
also  easily  ionized  by  the  electrons  and  photons  at  the  laser  wavelength. 

These  reactions  are  additional  HCl  depletion  mechanisms,  and  will  produce 
extra  electrons  which  can  also  dissociatively  attach  to  HCl. 

A  detail  modeling  of  the  xenon  manifolds  is  needed  to  fully  understand 
the  experimental  results,  and  at  the  same  time  it  must  be  able  to  predict  the 
nonequilibriuffl  distributions  observed  in  xenon  infrared  lasers.  Additional 
xenon  excited  state  density  measurements,  particularly  for  the  high  level 
xenon  manifolds,  would  also  be  very  useful  for  estimating  the  characteristic 
distribution  temperatures  and  to  provide  feedback  to  the  theoretical  model. 
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FIGURE  CAPTIONS 


Fig.  1 

Fig.  2 

Fig.  3 

Fig.  4 


Fig.  5 


Schematic  of  xenon  excited  state  density  measurement  system. 

Energy  level  diagram  for  Xe  from  65000  to  00000  cm~^.  The  5d,  7p, 
and  6d  states  are  not  shown.  The  four  transitions  measured  during 
the  experiment  2ire  indicated. 

♦ 

Time-dependent  AN  densities  of  a  XeCl  mixture  (1.5S  Xe,  balance  Ne 
to  3000  torr,  temperature  294  K)  for  an  initial  HCl  concentration  of 
0.32^.  Solid  and  open  circles  correspond  to  AN  (823.2  nm)  auid 
AN* (828.0  nm) ,  respectively.  The  dashed  curve  is  a  normalized  trace 
of  the  e-beam  current  waveform  (400  ns  FWHM) . 

♦ 

Time-dependent  AN  densities  of  a  XeCl  mixture  for  an  initial  HCl 
concentration  of  0.16^.  Other  gas  parameters  and  plot  parameters  are 

the  same  as  in  Fig.  3,  except  open  triangles  correspond  to  AN  (834.7 

♦ 

nm)  and  solid  squares  correspond  to  AN  of  the  neon  excited  state 

O 

(3p[5/2]2  -  3d [7/2]^)  transition  at  837.8  nm.  The  AN  of  the  neon 
excited  state  is  scaled  by  a  factor  of  10. 

♦ 

Time-dependent  AN  densities  of  a  XeCl  mixture  for  an  initial  HCl 

concentration  of  0.12%.  Other  gas  parameters  are  the  same  as  in 

♦ 

Fig.  3.  Solid  triangles  correspond  to  AN  (840.9  nm) . 


Fig.  6.  Time-dependent  AN  densities  of  a  XeCl  mixture  for  an  initial  HCl 
concentration  of  0.09%.  Other  gas  parameters  are  the  same  as  in 
Fig.  3.  Plot  parameters  are  the  same  as  in  Fig.  5. 

♦ 

Fig.  7.  Time-dependent  AN  densities  of  a  XeCl  mixture  for  an  initial  HCl 
concentration  of  0.08%.  Other  gas  parameters  are  the  same  as  in 
Fig.  3.  Plot  parameters  are  the  same  as  in  Figs.  3-5. 

* 

Fig.  8.  Time-dependent  AN  densities  of  a  XeCl  mixture  for  an  initial  HCl 
concentration  of  0.04%.  Other  gas  parameters  are  the  same  as  in 
Fig.  3.  Plot  parameters  are  the  same  as  in  Figs.  3  and  5. 

* 

Fig.  9.  Time-dependent  AN  densities  of  a  Ne/Xe  mixtures  (1.5  %  Xe,  balance 

Ne  to  3000  torr,  temperature  294  K) .  Plot  parameters  are  the  same  as 
in  Figs.  4  and  5.  Dotted  lines  with  solid  auid  open  circles  are 
AN*' (823. 2  nm)  and  AN*  (828.0  nm)  ,  respectively,  measured  for  an  Ar/Xe 
mixture  (1.5%  Xe,  balanced  Ar  to  1500  torr,  temperature  294  K) . 

Fig.  10.  Block  diagram  showing  key  kinetics  paths  in  the  Ne/Xe  rare  gas 
system. 

Fig.  11.  Fluorescence  time  history  in  the  wavelength  region  of  120<X<320  nm, 
obtained  from  a  Ne/Xe  mixture  (1.5  %  Xe,  balance  Ne  to  3000  torr, 
temperature  294  K) .  The  dashed  curve  is  a  normalised  trace  of  the  e- 
beam  current  waveform  (400  ns  FWHH) . 


Fig.  12.  Time-dependent  AN  densities  of  a  XeCl  mixture  for  an  initial  HCl 

concentration  of  0.32^.  Other  gas  psurameters  and  plot  parameters  are 
the  same  as  in  Fig.  3.  The  dashed  curve  is  a  normalized  trace  of  the 
e-beam  current  waveform  (1  /<s  FWHM) .  The  solid  curve  is  a  normalized 
trace  of  the  XeCl*  sidelight  fluorescence  waveform. 

Fig.  13.  Time-dependent  AN  densities  of  a  XeCl  mixture  for  an  initial  HCl 

concentration  of  0.16X.  Other  gas  parameters  are  the  same  as  in 

Fig.  3.  Plot  parameters  are  the  same  as  in  Figs.  3  and  5.  The  dashed 

curve  is  a  normalized  trace  of  the  electron  beam  current  waveform 

* 

(1  ns  FWHM) .  The  solid  curve  is  a  normalized  trace  of  the  XeCl 
sidelight  fluorescence  waveform. 

♦ 

Fig.  14.  Time-dependent  AN  densities  of  a  XeF  mixture  (0.5!(  Xe,  balance  Ne  to 
2888  torr,  temperature  294  K)  for  initial  F2  concentration  of  0.15%. 
The  dashed  curve  is  a  normalized  trace  of  the  e-beam  current  waveform 
(400  ns  FWHM).  Plot  parameters  are  the  same  as  in  Figs.  3  and  4. 

Fig.  15.  Time-dependent  AN  densities  of  a  XeF  mixture  for  initial  F2 

concentration  of  0.10%.  Other  gas  parameters  are  the  same  as  in 
Fig.  14.  The  dashed  curve  is  a  normalized  trace  of  the  e-beam 
current  waveform  (400  ns  FWHM) .  Plot  parameters  are  the  same  as  in 
Figs.  3-5. 

Fig.  16.  Time-dependent  AN  densities  of  a  XeF  mixture  for  initial  F2 
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ABSTRACT 


The  electron  energy  distribution  in  electron-beam  (e-beam)  excited  Ar/Kr 
and  Ne/Xe  gas  mixtures  is  examined  in  detail.  The  binary  rare  gas  mixtures 
are  similar  to  those  used  in  excimer  lasers.  Cooling  processes  for  the 
secondary  electrons  generated  in  the  gas  mixture  plasma  by  the  e-beam  are 
calculated  using  a  reduced  Boltzmann  equation  in  which  elastic  and  electron- 
electron  collisions  for  electron  energy  distributions  above  the  first 
excitation  threshold  of  the  rare  gas  are  ignored.  During  the  calculations  for 
the  Ar/Kr  and  Ne/Xe  mixtures,  all  electron  related  reactions  and  the 
interaction  between  the  two  different  rare  gases  in  the  mixture  are 
simultaneously  considered.  The  high  energy  secondary  electrons  produce  a 
steady-state  distribution  within  a  very  short  time;  however,  it  is  found  that 
the  distribution  is  not  Maxwellian.  W-values  [eV/electron-ion  pair]  and 


yields  of  rare-gas  excited  states  calculated  from  the  steady-state  high-energy 
electron  distribution  show  a  dependence  on  the  mixture  composition,  especially 
for  mixtures  with  low  concentrations  of  the  minor  rare  gas.  Thiu  implies  that 
the  practice  in  excimer  kinetics  models  of  using  the  W-values  determined  from 
pure  rare  gases  is  not  entirely  accurate. 


I.  INTIODUCTION 

Most  of  the  ions  and  excited  species  in  electron-beam  (e-beam)  generated 
plasmas  are  formed  either  directly  by  the  primary  and  hot  secondary  electrons, 
or  by  dissociative  recombination  of  electrons  with  molecular  ions  or  attaching 
species.  Therefore,  for  e-beam  pumped  gas  lasers,  the  kinetics  role  of  both 
the  electrons  and  the  heavy  particles  are  very  important  with  regard  to 
understanding  the  plasma  physics. 

The  motivation  of  this  paper  arises  from  recent  measurements  of  the 
time-dependent  electron  density  in  e-beam  pumped  rare-gas  halide  lasers^”^. 

O 

The  results  show  considerably  higher  electron  densities  than  that  estimated^ 
using  a  typical  value  for  the  rate  constant  for  dissociative  attachment  by  a 
halogen,  and  using  a  typical  rare  gas  W-value  (the  energy  expended  to  form  an 
electron-ion  pair) .  A  detailed  model  calculation^,  assuming  a  Maxwellian 
distribution  for  the  secondary  electrons  in  conjunction  with  an  energy  balance 
equation,  also  gives  lower  electron  densities.  To  help  understand  these 
differences  requires  a  more  careful  examination  of  the  electron  energy 
distribution.  The  objective  of  this  work  is  to  describe  in  detail  the 
electron  distribution  characteristics  in  these  laser  plasmas. 


To  understand  the  issues  related  to  the  usual  development  of  the  kinetic 
models,  it  is  helpful  to  review  the  energy  flow  channels  for  the  electrons  in 
an  e-beam  pumped  laser  mixture.  First,  the  primary  electrons  from  the  e-beam 
directly  produce  ions,  secondary  electrons,  and  excited  molecules.  This 


process  is  usually  accounted  ior  in  the  models  by  using  a  stopping  power 
constant  for  the  gases. ^  The  hot  secondary  electrons  are  cooled  by  collisions 
with  gas  molecules  and  produce  more  ions  and  excited  states.  This  process  is 
typically  modeled  by  using  W-values  and  the  yield  of  excited  states  relative 
to  the  ion  formation.  The  W-value  is  generally  higher  (typically  a  factor 
of  1.7-1. 8)  than  the  ionization  potential  of  the  atom.  The  primary  electrons 
from  the  e-beam  and  the  hot  secondary  electrons  are  usually  treated 
collectively  as  the  e-beam  source.  Because  the  preceding  two  processes  are 
very  fast  (<100  ps  for  typical  rare-gas  halide  laser  mixtures) ,  the  formation 
of  the  ions  and  excited  states  by  the  e-beam  source  has  been  treated  as  an 
instantaneous  process  by  simply  using  gross  macroscopic  values  for  the 
stopping  power  constant  and  W-values. 

The  cold  secondary  electrons  redistribute  their  energy  through  slow 
electron  reactions  such  as  dissociative  attachment  with  the  halogen, 
excitation  and  ionization  of  excited  states,  recombination  with  positive  ions, 
elastic  collisions  with  neutral  gases,  superelastic  collisions,  and  electron- 
electron  collisions.  This  redistribution  can  be  calculated  by  solving  the 
Boltzmann  transport  equation,  as  is  usually  done  for  discharge  pumped  lasers. 
However,  in  models  of  e-beam  pTimped  lasers,  the  distribution  of  the  secondary 
electrons  is  many  times  either  ignored  or  assumed  to  be  Maxwellian.  This  may 
lead  to  an  inaccurate  estimation  of  the  attachment  reaction  rate  of  the 
halogen  gas,  and  the  ionization  rates  by  the  hot  secondary  electrons  in 
populations  above  the  ionization  threshold  energy. 


This  paper  concentrates  primarily  on  the  high-energy  part  of  the 
electron  distribution  function.  The  high-energy  electron  distribution  is  the 
effective  driving  term  of  the  Boltzmann  equation  for  the  low-energy  electron 
distribution.  Calculations  of  the  high-energy  electron  distributions  have 
been  made  by  Elliot  and  Green  ,  and  Bretagne,  et  al.,  but  these  calculations 
were  limited  to  pure  rare  gases  of  argon  or  xenon.  Since  most  r are-gas  halide 
lasers  use  a  mixture  of  different  rare  gases  with  a  trace  of  halogen,  the 
calculations  need  to  include  the  distributions  of  inelastic  collision  cross- 
sections  for  each  of  the  rare  gases  as  a  function  of  energy.  Because 
competition  among  the  inelastic  collision  processes  (to  be  explained  later) 
between  different  rare  gases  can  be  accounted  for  only  through  a  simultaneous 
calculation,  the  electron  energy  distributions  and  the  W-values  of  mixtures 
cannot  be  estimated  simply  from  the  results  of  pure  rare  gases.  This  paper 
presents  calculations  for  mixtures  of  Ar/Kr  and  Ne/Xe,  which  are  typical  rare 
gas  combinations  for  e-beam  excited  KrF  (248  nm) ,  and  XeCl  (308  nm)  or  XeF 
(353  nm)  lasers,  respectively. 

Section  II  gives  a  brief  description  of  our  model.  Section  III  presents 
the  calculated  results  of  the  high-energy  electron  distributions  for  pure  rare 
gases  of  Ar,  Kr,  Ne,  Xe,  and  for  mixtures  of  Ar/Kr  and  Ne/Xe.  W-values  are 
calculated  for  mixtures  and  compared  with  those  of  pure  rare  gases.  The 
yields  of  excited  species  and  their  dependence  on  the  mixture  are  calculated. 
Also  presented  are  the  formation  rates  of  various  ions  and  excited  species 
produced  either  by  direct  e-beam  processes  or  by  the  steady-state  high-energy 


electron  distribution.  Conclusions  are  given  in  Section  TV. 


II.  MODELING  AND  DEEIVATION  OF  RATE  CONSTANTS 

Our  model  development  begins  by  using  the  same  basic  calculation 
procedure  as  presented  by  Bretagne,  et  al.^  The  distribution  of  electron 
energy  is  derived  from  a  reduced  Boltzmann  equation  with  no  elastic  and 
electron-electron  collisions: 

8f(e,t)  ran 

in 

where  f(e,t)  is  the  electron  energy  distribution  function  of  energy  e  at  a 
time  t.  The  first  term  on  the  right-hand  side  corresponds  to  the  production 
and  removal  of  electrons  with  energy  e  by  inelastic  collisions  (ionization, 
and  excitation),  recombination,  and  superelastic  de-excitation;  and  the  second 
term  corresponds  to  the  production  rate  of  secondary  electrons  with  energy  e 
by  primary  electrons  in  the  e-beam  with  energy  S(ep,e)  is  the 

corresponding  differential  ionization  cross-section,  and  I(t)  is  a  value  which 
is  proportional  to  the  current  density  of  the  e-beam. 

The  effects  of  electron-electron  collisions  can  be  neglected  as  long  as 
the  collisional  rates  of  the  inelastic  reactions  are  much  higher  than  the 

electron-electron  collision  rate.  The  value  ♦(€)  is  defined  by  the  ratio  of 
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the  electron-electron  collision  rate  to  the  inelastic  collision  rate: 
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where  is  the  secondary  electron  density,  n^  is  the  rare  gas  density  '  ®in 

9 

the  inelastic  collision  cross-section,  and  A  is  the  Spitzer  parameter.  The 

-3 

value  of  ♦(£)  is  less  than  10  for  electrons  whose  energies  are  higher  than 
the  first  excitation  threshold  of  the  major  gas  constituent.  Therefore,  as 
long  as  the  minimum  energy  considered  is  higher  than  the  first  excitation 
threshold,  the  influence  of  electron-electron  collisions  can  be  neglected. 
However,  if  the  concentration  of  the  minor  rare  gas  is  extremely  low,  then  the 
formation  of  its  ion  and  excited  states  can  be  affected  by  electron-electron 
collisions.  Hence,  the  electron  cooling  processes  involving  inelastic 
collisions  with  the  minor  rare  gas  must  be  treated  as  a  slow  process  if  the 
rare  gas  concentration  is  low  (e.g.  <1?{)  .  This  will  be  discussed  later.  In 
this  work,  the  energy  distribution  above  the  first  excitation  threshold  of  the 
major  rare  gas  is  calculated  as  a  steady-state  high-energy  electron 
distribution  function. 


Similarly,  the  elastic  collision  term  can  be  neglected  as  seen  by 
evaluating  the  function  0(6): 


V 


where  m  /M  is  the  electron  to  rare  gas  mass  ratio  and  Q  ,  is  the  elastic 
collision  cross-section.  Contributions  of  inelastic  collisions  involving 
excited  molecules  are  ignored,  due  to  the  relatively  low  density  of  these 
molecules,  during  the  calculation  of  the  high-energy  secondary  electron 
distribution. 

The  energy  distribution  of  the  secondary  electrons  is  calculated  using 
differential  ionization  cross-sections.  Peterson  and  Allen^^,  Green  and 
Sawada  ,  and  Bretagne,  et  al.  have  proposed  empirical  formulas  for  M-shell 
ionization  of  argon  gas.  For  other  rare  gases,  the  formulas  given  by  Green 
and  Sawada^^  are  available.  The  results  of  calculating  the  total  ionization 
cross-sections  and  energy  losses  for  various  incident  electron  energies  using 
the  different  empirical  formulas  for  argon  gas  are  shown  in  Table  I.  Although 
the  three  different  formulas  give  similar  values  for  the  total  ionization 
cross-sections  over  a  wide  range  of  incident  electron  energies,  the  energy 
loss  values  calculated  by  Peterson  and  Allen’s  formula  are  significantly 
larger  than  the  others.  This  is  because  their  formula  predicts  larger 
differential  ionization  cross-sections  for  the  high  energy  secondary  electrons 
which  leads  to  a  greater  energy  loss.  For  comparison,  the  energy  loss  values 
calculated  using  Bethe’s  stopping  power  theory  ,  which  provides  a  macroscopic 
calculation  of  energy  loss  values  for  high  energy  electrons,  are  also  listed 
in  Table  I.  Even  the  values  calculated  from  Peterson  and  Allen’s  formula  are 
lower  than  those  obtained  from  Bethe’s  theory. 


MME  (MS  Hm.IK  (IMH)  KIHETICS(U>  SKCTM  TECMMLOtV 
INC  BELLEVUE  HA  H  D  KIHUfifl  FEB  B8  C-11-3I3 


The  discrepancy  between  the  values  calculated  by  Bethe’s  formula  and 

those  of  the  empirical  differential  ionization  cross-sections  can  be  explained 

as  a  result  of  energy  loss  caused  by  inner  shell  ionizations.  The  L-shell 

ionization  loss  calculated  from  an  empirical  differential  ionization  cross- 

10  ®2  4 

section  proposed  by  Peterson  and  Allen  is  «0.5  eV-A  for  a  10  eV  primary 

energy,  which  is  6%  of  the  M-shell  ionization  loss  given  in  Table  I. 

Recently,  Bretagne,  et  al.  have  also  reported  an  empirical  differential 

ionization  cross-section  for  the  L-shell  of  argon,  and  an  energy  loss  of  3  eV- 

A^  is  obtained  from  their  formula  for  a  10^  eV  primary  energy.  Combining  this 

®2 

L-shell  ionization  loss  with  the  M-shell  ionization  loss  of  5  eV-A  ,  listed  in 

®2 

Table  I  for  Bretagne,  et  al.,  yields  a  total  ionization  loss  of  ^^8  eV-A  . 

This  is  close  to  the  value  calculated  from  Peterson’s  formula  (see  Table  I), 
but  still  «20?{  lower  than  Bethe’s  loss  value. 


Bethe’s  formula  calculates  the  energy  loss  directly  and  is  generally  in 
good  agreement  with  experiments;  whereas  the  method  using  the  empirical 
formulas  for  the  differential  ionization  cross-sections  is  sm  indirect 
approach  and  is  sensitive  to  the  values  of  the  cross-sections  at  high 
secondary  electron  energies.  The  empirical  formulas  all  tend  to  underpredict 
the  differential  ionization  cross-sections  for  secondary  electron  energies 
ilOO  eV.®  This  leads  to  an  underprediction  of  the  stopping  power. 


These  disagreements  regarding  the  stopping  power  obtained  from  the 
empirical  differential  ionization  cross-sections  become  more  pronounced  for  Kr 
and  Xe.  When  the  empirical  formulas  proposed  by  Green  and  Sawada^^  are  used, 


"v  y 


shell  ionisation  rate  becomes  larger  for  high  Z  atoms.  The  total  ionization 
cross-section  for  inner  shell  ionization  is  very  small  compared  with  the  outer 
shell;  however,  its  contribution  to  the  stopping  power  is  significant  due  to 
its  high  ionization  energy.  Unfortunately,  to  our  knowledge,  differential 
ionization  cross-sections  for  Kr  and  Xe  are  not  available. 


Although  the  empirical  formulas  for  the  differential  ionization  cross- 
section  suffer  from  inaccuracies  at  high  secondary  electron  energies,  they  are 
still  useful  for  estimating  the  W-value  of  a  gas.  The  W-value  is  given  by  the 
ratio  of  the  rare  gas  ion  density  divided  by  the  stopping  energy  for  electrons 
in  the  gas.  The  empirical  formulas  tend  to  underpredict  the  stopping  power, 
but  they  also  tend  to  underpredict  the  ion  density  by  approximately  the  same 
amount.  This  is  because  the  ion  density  prediction  initially  uses  the 
electron  distribution  calculated  for  the  primary  electrons,  but  this  same 
distribution  relies  on  the  stopping  power  prediction.  Peterson,  et  al.^^  and 
Bretagne,  et  al,  ’  have  found  that  the  net  formation  efficiencies  of  ions 
and  excited  states  are  nearly  constant  for  initial  electron  energies  greater 
than  ^50  eV  [net  formation  efficiency  is  defined  here  as  the  ratio  of  the 
ionization  (excitation)  energy  divided  by  the  total  energy  loss  per  ion 
(excited  state)  formed] .  This  implies  that  the  ion  or  excited  state  formation 
efficiency  is  relatively  independent  of  the  shape  of  the  electron  distribution 
during  the  continuous  process  of  energy  cooling.  Hence,  the  ratio  of  the  ion 
density  to  the  stopping  power  using  the  empirical  formulas  is  approximately 
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equal  to  the  ratio  of  the  real  ion  density  to  the  real  stopping  power. 

This  assumption  was  validated  by  computing  the  W-value  using  the  three 

different  empirical  formulas:  1)  Peterson  and  Allen^®,  2)  Bretagne,  et  ai.® 
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without  L-shell  ionization,  and  3)  Bretagne,  et  el.  with  L-shell 
ionization.  The  obtained  W-values  for  pure  argon  are  29,  25,  and  27  eV, 
respectively.  Hence,  despite  the  large  differences  among  values  for  stopping 
power  (of  order  50^),  the  predicted  W~values  are  within  15%  of  each  other. 
Moreover,  the  formation  efficiencies  of  the  excited  species  are  also 
approximately  equal  for  all  three  cases.  In  this  paper,  it  is  assumed  that 
the  other  rare  gases  also  follow  this  general  trend. 

In  order  to  estimate  the  secondary  electron  density  in  a  pure  rare  gas, 
it  is  necessary  to  correct  the  density  calculated  using  the  energy  loss  result 
from  the  empirical  formulas  by  multiplying  the  density  by  the  ratio  of  the 
stopping  power  predicted  by  Bethe’s  formula  over  the  stopping  power  predicted 
by  the  empirical  formulas.  For  this  method  of  correction  to  be  valid  it  is 
necessary  that  the  electron  energy  distribution  be  independent  of  the 
excitation  or  deposition  rate.  This  is  true  when  the  effects  of  electron- 
electron  collisions  can  be  neglected,  which,  as  was  discussed  earlier,  is  true 
for  the  cases  studied  here. 

As  mentioned,  the  preceding  method  for  calculating  the  secondary 
electron  density  and  energy  distribution  is  applicable  to  the  case  of  a  pure 
rare  gas.  When  the  medium  consists  of  a  mixture  of  gases,  the  situation 


becomes  considerably  more  complex  because  the  energy  losses  predicted  by  both 
Bethe’s  formula  and  the  empirical  ones  depend  on  the  type  of  gas. 

The  low  total  ionization  cross-section  of  the  inner  shell  means  that  the 
density  of  inner  shell  ions  is  very  low  compared  to  outer  shell  ions. 

However,  the  inner  shell  ions  can  eject  Auger  electrons  with  high  enough 
energy  to  produce  additional  excitations  and  ionizations. Therefore,  the 
inner  shell  ionization  process  effectively  corresponds  to  the  production  of 
Auger  electrons  at  a  rate  defined  by  the  inner  shell  total  ionization  cross- 
section.  The  produced  Auger  electrons  are  cooled  through  the  formation  of 
ions  and  excited  atoms  at  probabilities  given  by  the  reaction  rates  at  the 
energy  of  the  Auger  electrons.  The  ratios  of  the  ionization  cross-sections 
for  Ar  and  Kr  at  10^  eV  (primary  electron  energy)  and  200  eV  (approximate 
Auger  electron  energy  in  Ar)  are  similar.  Bence,  Auger  electrons  produce  Ar* 
and  Kr^  at  approximately  the  same  probability  as  the  primary  electrons.  It  is 
therefore  possible  to  account  for  the  effects  of  Auger  electrons  by  increasing 
the  effective  total  cross-section  of  the  primary  electrons. 

In  this  paper,  rather  than  attempt  to  compensate  for  the  ionization 
cross-section  differences  at  the  high  energy  part  of  the  electron 
distribution,  by  assuming  an  empirical  differential  ionization  cross-section 
for  the  inner  shell,  a  simpler  approach  is  taken.  The  total  ionization  cross- 
sections  given  by  the  empirical  formulas  and  the  excitation  cross-sections  for 
each  of  the  rare  gas  species  are  multiplied  by  a  constant  factor  such  that  the 
resultant  energy  loss  agrees  with  Bethe’s  theory  for  rare  gas  mixtures. 


The  preceding  procedure  is  performed  only  on  the  differential  ionization 
cross-sections  for  the  primary  electrons.  Thereafter,  this  compensation 
procedure  does  not  have  to  be  applied  to  the  hot  secondaries  and  their  cooler 
offspring.  Instead,  the  uncorrected  empirical  formulas  are  used  to  predict 
the  energy  loss  within  a  computational  cell  of  the  medium  as  the  hot 
secondaries  produce  additional  but  cooler  electrons.  Since  the  electron 
energy  is  conserved  within  the  cell,  the  error  in  the  energy  loss  prediction 
does  not  cause  an  error  in  the  energy  balance  within  the  cell. 


The  total  ionization  cross-sections  obtained  from  the  empirical 
differential  ionization  cross-sections  for  the  outer  shell  agree  well  with 
experimental  values.  Therefore,  minimal  error  is  introduced  by  ignoring  the 

inner  shell  ionizations  when  calculating  the  ionizations  and  excitations 
occurring  in  mixtures  of  rare  gases,  except  for  the  error  caused  by  the  change 
in  the  secondary  electron  distribution  during  the  continuous  process  of  energy 
cooling.  Since,  as  was  shown  for  pure  argon,  this  latter  error  tends  to  be 
very  small  in  a  pure  rare  gas,  it  is  assumed  in  this  paper  that  this  trend  is 
also  true  for  rare  gas  mixtures.  This  last  assumption,  however,  has  not  been 
validated. 


For  the  excitation  cross-sections,  the  expressions  given  by  Peterson  and 

Allen^^  for  argon  gas  are  used.  For  neon,  krypton,  and  xenon,  the  empirical 
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cross-sections  of  Ganas  and  Green  and  the  formula  proposed  by  Vriens  and 
18 

Smeets^  are  used  for  optically  allowed  transitions.  At  sufficiently  large 
electron  energies,  the  cross-section  of  optically  allowed  excitations 
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Qaiioa^^)  approaches  the  form  : 
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where  Kj  is  a  constant.  Since  there  is  no  data  available  (except  for  argon), 
the  cross-section  for  the  optically  forbidden  np-(n+l)p  transitions  are 


assumed  to  be  of  the  form: 


1  -  e.v/e 


'forbid'^’  =  "2 


where  K2  is  a  constant,  and  is  the  threshold  energy  for  the  excitation. 

At  sufficiently  high  energy,  Eq.  (5)  approaches  Kg/e.  During  calculations  in 

pure  rare  gases,  the  values  of  Kg  for  neon,  krypton,  and  xenon  are  adjusted  so 
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as  to  yield  W-values  consistent  with  the  experimentally  determined  values.  ’ 
The  cross-section  distributions  for  inelastic  collisions  used  during  the 
calculations  are  shown  in  Figs.  1  and  2  for  Ar/Kr  and  Ne/Xe  mixtures, 
respectively.  The  electronic  classifications  and  their  energy  levels  which 
are  included  in  the  model  are  listed  in  Table  II. 


The  results  presented  in  this  paper  are  for  a  fixed  primary  electron 

energy  of  10  keV.  Of  course  in  real  laser  systems  the  e-beam  energy  is 

typically  100 's  of  keV  to  several  MeV.  Nonetheless,  calculations  of  the 

energy  dependence  to  create  ions  and  excited  states  from  that  part  of  the 

distribution  above  10  keV  yield  similar  results  even  at  1  MeV  energy 
9  10  13 

levels.  ’  ’  Thus,  the  assumption  can  be  made  that  the  ionization  and 

excitation  efficiencies  for  primary  electron  energies  above  10  keV  is  constant 
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and  equal  to  those  obtained  at  10  keV.  This  approximation  allows  the  number 
of  secondary  energy  segments  considered  in  the  model  to  be  minimized,  thereby 
reducing  the  computation  time. 

The  results  presented  in  this  paper  are  for  a  constant  excitation  rate. 

Changes  in  the  current  density  affect  the  overall  magnitude  of  the 

3 

distribution  function  [#/cm  eV]  without  affecting  the  overall  shape  of  the 
high  energy  region.  During  the  calculations  presented  next,  the  system 
represented  by  Eq. (1)  is  simultaneously  integrated  over  «250  energy  segments 
using  a  fourth-degree  Runge-Kutta  method.  A  provision  for  increasing  the 
widths  of  the  energy  intervals  with  increasing  energy  is  introduced  to  reduce 
the  number  of  segments  required  and  to  keep  the  accuracy  high  during 
calculations  at  the  lower  energy  levels. 

III.  RESULTS  AND  DISCUSSION 

A.  Distribution  of  High-Energy  Secondary  Electrons 

Depicted  in  Fig.  3  is  the  temporal  evolution  of  the  electron  energy 
distribution  calculated  for  a  3000  Torr  pure  neon  case.  The  absolute  number 
density  of  the  secondary  electrons  shown  in  this  figure  has  already  been 
adjusted  to  account  for  the  error  in  stopping  power  given  by  the  empirical 
formulas  for  the  differential  ionization  cross-sections  against  that 
calculated  by  Bethe’s  theory.  The  calculation  continued  until  the  electron 


densities  around  20  eV  and  greater  reached  constant  values.  This  final 
distribution  represents  the  steady-state  distribution.  The  electron  densities 
tend  to  reach  steady-state  values  from  the  direction  of  the  higher  energy 
side.  This  is  because  the  densities  at  each  energy  are  achieved  by  a  balance 
between  a  formation  term,  which  includes  direct  formation  by  the  primary 
electrons  and  by  electrons  at  high  energies,  and  a  loss  term  due  to  inelastic 
collisions  by  neutrals  and  electrons. 

Note  how  in  Fig.  3  the  steady-state  distributions  rapidly  decrease 
beyond  the  threshold  of  the  first  excitation  (16.7  eV)  of  neon.  Densities 
lower  than  ^17  eV  are  plotted  by  a  dashed  line  to  indicate  an  uncertainty  in 
the  calculation  for  this  part  of  the  curve.  This  is  because  the  densities  at 
these  lower  energies  are  sensitive  to  slow  electron  reactions  (e.g., 
recombination,  electron-electron  collisions,  ionization  and  excitation  of 
excited  molecules,  and  superelastic  de-excitations),  which  have  not  been 
included  in  the  calculation  yet.  Detailed  results  for  the  low-energy  electron 
distribution  are  reported  elsewhere.  When  calculating  the  distribution 
function  for  the  lower  energy  part,  the  previously  obtained  steady-state 
distribution  of  the  high  energy  part  (€>17  eV)  is  used  as  the  effective  e-beam 
source  (i.e.,  the  driving  term  in  the  Boltzmann  equation  for  the  low-energy 
electrons) .  Formation  rates  of  ionized  or  excited  species  and  of  lower  energy 
electrons  are  calculated  from  this  e-beam  source,  and  are  discussed  in  more 
detail  in  Section  III(B). 


The  normalized  steady-state  distribution  of  the  high  energy  tail  for 
3000  Torr  neon  is  shown  in  Fig.  4.  This  is  compared  with  a  Maxwellian 
distribution  at  a  characteristic  temperature  kT^  (=  2/3  <e>)  of  4.6  eV,  which 
corresponds  to  the  average  energy  of  the  steady-state  distribution.  The 
steady-state  distribution  has  a  considerably  laurger  number  of  electrons  at 
energies  >40  eV. 

The  mean  energy  required  for  the  formation  of  an  ion-electron  pair  can 

be  estimated  from  the  ion  formation  rate,  calculated  using  the  steady-state 

distribution,  ratioed  to  the  e-beam  deposition  rate.  This  value  effectively 

corresponds  to  the  W-value.  The  energies  required  for  the  formation  of 

various  excited  states  (also  referred  to  as  the  W-values  for  excited  states) 

are  obtained  in  the  same  manner.  The  formation  of  a  neon  ion-electron  pair 

requires  38  eV,  and  Ne*,  Ne**,  Ne***,  and  Ne****  require  1.6  x  10^,  1.3  x  10^, 
3  3 

1.5  X  10  ,  and  9.0  x  10  eV  for  formation,  respectively  (see  Table  III).  This 
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calculated  W-value  of  38  eV  is  very  close  to  the  reported  value  of  36.6  eV. 

The  yield  of  total  excited  states  created  during  the  formation  of  one  ion- 

electron  pair  is  wO.56.  This  value  is  significantly  higher  than  the  value  of 

21 

0.33  reported  by  Blauer,  et  al.  They  obtained  this  value  from  the  ratio  of 
excitation  and  ionization  cross-sections,  and  ignored  the  highly  excited 
states  (Ne**,  Ne***,  Ne****) .  Jauicaitis^^  obtained  a  ratio  of  0.47  by  solving 
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the  Fowler  equation^  ’  including  several  excitation  levels.  The  energies 
required  for  the  formation  of  ions  or  excited  states  calculated  from  the 
Fowler  equation  were  obtained  during  the  cooling  processes.  This  means 
Jancaitis’  calculation  used  a  slightly  different  definition  from  our  model. 


HoveYer,  it  is  clear  that  a  considerable  number  of  excited  states  at  higher 
levels  are  created  together  with  the  excited  state  at  the  first  excitation 
level.  If  collisions  with  neutral  particles  efficiently  deactivates  higher 
level  excited  states  to  the  first  or  second  excitation  levels,  this  would 
significantly  impact  the  laser  kinetics. 

The  calculated  secondary  electron  distributions  for  760  Torr  xenon  are 

shown  in  Fig.  5.  Because  of  xenon’s  larger  cross-section  for  inelastic 

collisions,  the  steady-state  distribution  is  achieved  in  a  shorter  time  than 

the  3000  Torr  neon  case,  even  at  low  pressures.  The  distribution  is,  again, 

truncated  near  the  threshold  of  xenon  excitation  and  a  large  population  of 

electrons  reside  in  the  high  energy  tail.  The  calculated  W-value  is  22  eV, 
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which  is  in  good  agreement  with  the  reported  value  of  21.7  eV.  Energies 

2  2 

required  for  the  formation  of  excited  states  are  1.3  x  10  ,  1.2  x  10  ,  and 
9.1  X  10^  eV  for  Xe*,  Xe**,  auid  Xe***,  respectively  (see  Table  III).  The 
total  yield  of  excited  states  during  the  formation  of  an  ion-electron  pair  is 
0.38. 


Results  calculated  for  neon  gas  containing  0.5%  xenon  at  a  total 
pressure  of  3000  Torr  are  given  in  Fig.  6.  The  distribution  has  been 
calculated  up  to  100  ps;  however,  the  distribution  near  the  excitation 
threshold  of  Xe  has  not  reached  steady-state  yet.  Calculations  beyond  KlOO  ps 
require  the  inclusion  of  slow,  low  energy  electron  reactions,  such  as 
electron-electron  collisions  and  superelastic  collisions  that  can  form  hot 
electrons  from  excited  states.  These  reactions  have  not  been  included  in  this 
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part  of  the  calculation.  Two  shoulders  in  the  shape  of  the  distribution 
function  can  be  seen  near  the  excitation  thresholds  of  neon  and  xenon.  The 
distributions  beyond  about  17  eV  where  neon  dominates  have  reached  a  steady- 
state  and  are  about  the  same  as  those  of  pure  neon  shown  in  Fig.  4.  Hence, 
accurate  estimations  of  the  electron  energy  distribution  lower  than  the  first 
excitation  threshold  of  neon  need  to  be  calculated  separately,  with  the 
distributions  for  the  high  energy  electrons  still  valid  using  the  preceding 
technique.  This  is  particularly  true  for  mixtures  with  low  concentrations  of 
the  minor  rare  gas.  If  the  partial  pressure  of  the  minor  rare-gas  is  high 
greater  than  a  few  percent),  the  lowest  boundary  of  the  high-energy 
electron  distribution  can  be  set  near  the  excitation  threshold  of  the  minor 
rare  gas.  An  example  of  this  is  demonstrated  by  the  dotted  curve  in  Fig.  6, 
which  is  for  5%  Xe  in  neon  at  3000  Torr. 

Electron  cooling  in  these  e-beam  pumped  plasmas  involve  the  following 
processes.  Primary  electrons  produce  neon  ions  and  hot  secondary  electrons 
(<100  eY) .  These  are  estimated  using  the  distributions  and  magnitudes  of  the 
inelastic  collision  cross-sections.  The  hot  secondary  electrons  continue  to 
create  neon  ions,  thermalize  in  energy,  and  then  create  neon  excited  states 
and  warm  electrons  (<17  eV) .  The  very  small  concentration  of  xenon  means  the 
effective  formation  of  xenon  ions  and  excited  states  begins  after  the  warm 
electrons  have  cooled  down  to  less  than  the  threshold  of  neon  excited  state 
formation.  At  this  point  the  warm  electrons  create  cool  secondary  electrons 
(<8.4  eV)  during  the  formation  of  xenon  ions  and  excited  states.  Slow  low- 
energy  reactions  then  create  cold  electrons  (of  order  few  eV)  from  the  cool 


electrons. 

Typically,  the  ion  and  excited  state  formation  of  different  gases  in 
these  mixtures  are  treated  in  kinetic  models  as  individual  pure  gases  using 
their  respective  W-values  and  excited  state  yields  obtained  from  measurements 
or  estimations  for  pure  gases.  This  is  particularly  true  in  any  of  the  rare- 
gas  halide  laser  models  excited  by  an  e-beam.  However,  the  secondary 
electrons  created  by  collisions  with  different  gases  cannot  be  distinguished 
individually  in  the  mixture.  The  electron  energy  distribution  established  by 
the  major  gas  acts  as  an  effective  primary  electron  distribution  for  the  minor 
gas.  Therefore,  the  W-values  and  excited  state  yields  obtained  for  mixtures 
must  include  the  competition  among  the  vairious  inelastic  collisions  occurring 
among  the  different  gases.  These  are  included  exactly  in  our  calculations. 

For  the  0.5%  xenon  gas  mixture,  W-values  of  39  and  3.5  x  10  eV,  for  neon-ion 
and  xenon-ion,  respectively  are  obtained.  Since  the  basis  of  the  W-value 
calculation  is  the  total  deposition  rate  of  the  e-beam,  xenon  at  low 
concentrations  will  have  extremely  high  W-values. 

Although  it  lacks  any  physical  meaning,  it  is  useful  to  calculate  the  W- 
value  normalized  by  the  fraction  of  e-beam  energy  deposited  in  each  of  the  gas 
constituents.  This  is  referred  to  as  the  "effective*  W-value  in  this  paper. 
Values  of  6.3  eV  and  38  eV  are  obtained  for  xenon  and  neon,  respectively. 

These  are  the  values  which  should  be  used  in  the  laser  kinetics  model  instead 
of  the  W-values  obtained  for  pure  rare  gases.  However,  these  values, 
especially  for  xenon,  strongly  depend  on  the  minor  gas  partial  pressure.  The 


sane  calculations  are  repeated  for  a  mixture  of  Ne/Xe  =  2850/150  Torr. 

Cooling  of  the  secondary  electrons  above  the  first  excitation  threshold  of 
xenon  becomes  accelerated  because  of  the  higher  density  of  xenon.  The 
calculated  steady-state  distribution  is  given  in  Fig.  6  by  the  dotted  curve. 

It  is  very  similar  to  the  distribution  for  pure  xenon  despite  the  low  xenon 
pressure.  For  this  mixture,  the  normalized  W-value  with  the  partitioned 
deposition  energy  for  xenon  increases  to  13  eV,  but  this  is  still  9  eV  smaller 
than  the  value  for  pure  xenon.  Despite  the  differences  between  the  W-values, 
this  does  not  have  a  strong  influence  on  the  rare-gas  halide  excimer  formation 
because  the  dominant  channel  of  xenon  ion  formation  in  the  Ne/Xe  mixtures  at 
low  xenon  concentration  is  Penning  ionization  by  neon  excited  states  and 
charge  transfer  from  neon  dimer  ions. 


The  energies  required  for  the  formation  of  each  excited  state  in  a  0.5% 
Xe  mixture  are  6.4  x  10^,  4.6  x  10^,  and  6.5  x  10^  eV  for  Xe*,  Xe**,  and 
Xe***,  respectively;  and  1.7  x  10^,  1.4  x  10^,  1.5  x  10^,  and  9.2  x  10^  eV  for 
Ne*,  Ne**,  Ne***,  and  Ne****,  respectively.  These  are  listed  in  Table  III 
together  with  the  results  for  pure  rare  gases.  Efficiencies  calculated  by  the 
ratio  of  the  W-value  to  the  ionization  or  excitation  threshold  energy  are  also 
listed.  The  results  of  calculating  the  total  excited  state  yields  are: 
}](Xe*/Xe^)  =1.4,  and  J^(Ne*/Ne^)  =  5.3  x  10  While  the  addition  of  a  small 
quantity  of  xenon  to  neon  has  only  a  slight  affect  on  the  total  excited  state 
yield  of  neon,  the  total  number  of  xenon  excited  states  created  is  now  larger 
than  the  ion  yield.  This  is  because  the  xenon  ionization  must  compete  with 
excitation  of  neon  excited  states  due  to  the  overlap  of  their  cross-sections 
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in  the  secondary  electron  spectrum  as  shown  in  Fig.  2.  This  results  in  a 
decrease  of  xenon  ion  formation  rate.  Using  these  results,  the  total  number 
of  species  in  a  0.5%  Xe  mixture  created  by  the  e-beam  source,  which  act  as  the 
precursors  for  rare  gas  halide  formation,  can  be  calculated  from: 


39  eV  (LOSS)  1  Ne""  +  0.76  +  0.11  Xe"  +  0.15 


[6] 


where  ^®total  bhe  total  population  in  the  various 

excitation  levels.  Hence,  an  e-beam  loss  of  30  eV  yields  approximately  2.0 
precursors.  Increasing  the  xenon  concentration  to  5%,  the  total  excited  state 
yield  of  xenon  relative  to  the  ion  formation  decreases  to  0.01  because  the 
formation  efficiency  of  xenon  ion  improves.  In  this  case,  an  e-beam  loss  of 
45  eV  results  in  the  formation  of  2.5  precursors.  It  should  be  mentioned  that 
increasing  the  xenon  concentration  tends  to  degrade  the  laser  performance  due 
to  formation  of  triatomics. 


It  is  difficult  to  accurately  estimate  the  W-values  of  the  minor  rare 
gas  excited  states,  in  this  case  xenon,  due  to  the  fact  that  the  distribution 
function  near  the  first  threshold  of  xenon  cannot  be  accurately  calculated  by 
considering  only  the  fast  reactions  with  hot  electrons.  (Although  the  error 
in  the  values  listed  in  this  paper  is  negligible  (<1%),  the  degree  of  error 
increases  with  decreasing  concentration  of  the  minor  rare  gas.)  The  root  of 
the  problem  lies  in  the  definition  of  the  W-value  itself.  As  described 
earlier,  the  cooling  process  that  the  hot  secondary  electrons  undergo  does  not 
occur  in  well-defined  discrete  step,  but  is  a  continuous  process  with  no  clear 
boundary  between  the  fast  reactions  occurring  with  hot  electrons  and  the  slow 


reactions  occurring  with  cool  electrons.  As  an  example,  for  the  minor  rare 
gas,  the  excited  state  formation  is  affected  by  slow  reactions,  such  as 
electron-electron  and  superelastic  collisions,  which  can  cause  a 
redistribution  of  the  hot  secondary  electrons  whose  energy  are  larger  thaui  the 
excitation  threshold.  In  this  manner,  the  W-value  of  excited  state  formation 
and  their  relative  yield  are  a  function  of  the  e-beam  deposition  rate, 
particularly  for  the  minor  rare  gas. 

Since  the  total  electron  and  excited  state  densities,  and  the  rate  of 
electron-electron  collision  and  superelastic  collision  are  affected  by  the 
deposition  rate,  it  should  be  again  emphasized  that  the  W-values  and  excited 
state  yields  for  the  minor  gas  in  the  mixture  cannot  be  accurately  determined 
and  are  obviously  different  from  the  values  for  a  pure  gas.  The  overall 
magnitude  of  the  values  calculated  in  this  paper,  however,  are  still  valid. 
This  observation  regarding  W-values  suggests  that  the  usual  method  for 
calculating  the  formation  of  ions  and  excited  states  for  the  minor  gases  in  a 
mixture  by  simply  using  W-values  and  relative  yields  is  limited  by  an  inherent 
inaccuracy. 

Similar  calculations  are  repeated  for  1020  Torr  pure  argon,  760  Torr 
pure  krypton  (Fig.  7),  and  various  mixtures  of  Ar/Kr  (Fig.  8).  Table  IV  lists 
the  calculated  effective  W-values  and  excited  state  yields.  The  overall 
characteristics  of  the  distribution  functions  are  very  similar  to  those  for 
neon  and  xenon  gases.  The  major  difference  with  the  neon  and  xenon  mixture  is 
that  all  the  threshold  energies  for  ionization  and  excitation  of  argon  and 
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krypton  aure  very  close  to  one  another.  They  are  located  between  about  10  auid 
16  eV.  Hence,  only  one  shoulder  is  seen  in  the  distribution  functions  of  the 
Ar/Kr  mixtures.  Since  not  only  the  cross-section  for  krypton  ionisation,  but 
also  the  excitation  of  krypton  share  the  same  secondary  electron  spectrum  with 
the  excitation  and  ionization  cross-sections  for  argon,  the  excited  state 
yield  of  Kr  obtained  from  a  mixture  of  Ar/Kr  =  1020/114  Torr  is  only  slightly 
higher  than  that  obtained  in  pure  krypton.  The  number  of  precursors  created 
by  the  e-beam  source  for  this  particular  mixture  is  1.7  with  an  e-beam  loss  of 
29  eV. 


KrF  laser  performance  has  been  investigated  with  mixtures  at  high 
25-29 

krypton  concentrations.  It  would  be  very  useful  from  a  laser  kinetics 

modeling  viewpoint  to  analyze  the  W-value  and  excited  state  yield  as  a 
function  of  krypton  concentration.  Figures  9  and  10  show  calculated  W-values, 
normalized  by  the  partitioned  deposition  energy  for  Ar  and  Kr,  and  excitation 
yields  of  individual  states  relative  to  the  ion  formation  as  a  function  of  the 
Ar  and  Kr  concentration,  respectively.  The  nature  of  the  competition  between 
the  elastic  collision  reactions  which  share  the  same  secondary  electron  energy 
spectrum  change  with  mixture  composition;  hence,  both  the  W-value  and  excited 
state  yield  are  not  unique  constants,  but  depend  on  the  mixture  composition. 
Again,  the  accuracy  of  the  estimations  for  the  excited  state  yields  drops  as 
the  minor  gas  concentration  decreases . 


m 
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B.  Fornation  Rates  of  Low-energy  Electrons  and  Various  Species 

As  described  earlier,  a  full  Boltzmann  code  calculation  is  necessary  to 
calculate  the  low  energy  distribution  of  the  secondary  electrons  in  order  to 
obtain  accurate  rate  constants  for  the  secondary  electron  reactions.  This  is 
especially  true  with  regard  to  dissociative  attachment  by  halogen  gases,  which 
can  affect  the  low  energy  part  of  the  electron  distribution.  The  effective 
role  of  the  high-energy  electron  distribution  function  is  as  the  ionization 
and  excitation  source  (referred  to  in  this  paper  as  the  e-beam  source)  in  the 
Boltzmann  equation.  Calculated  are  the  formation  rates  and  distribution  of 
the  low-energy  secondary  electrons  and  the  formation  rates  of  various  species 
using  the  calculated  steady-state  high-energy  electron  distribution  and  the 
primary  electrons  in  the  e-beam.  The  calculation  discussed  in  this  paper  used 
a  high-energy  electron  distribution  calculated  for  mixture  containing  0.5%  Xe 
in  neon  at  a  total  pressure  of  3000  Torr. 

During  the  high-energy  distribution  calculation,  the  effect  of  the 
halogen  gas  is  ignored  due  to  its  negligible  contribution;  however,  the 
formation  rates  of  the  halogen  ion  and  the  excited  halogen  molecule  are 
calculated.  The  approach  used  in  this  paper  divides  the  entire  secondary 
electron  distribution  into  a  high-energy  part  and  a  low-energy  part  with  the 
division  point  at  16.7  eV  corresponding  to  the  excitation  threshold  of  neon. 
This  means  the  formation  rates  of  the  low-energy  electrons  and  species  that 
are  described  next  are  the  rates  due  to  the  high-energy  electrons  above 
16.7  eV  and  by  the  primary  electrons. 
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To  calculate  the  distribution  of  the  low-energy  electrons  produced  by 
ionization  of  excited  states  and  halogen  molecules,  an  empirical  formula  for 
the  differential  ionization  cross-sections  must  be  used  that  also  gives  good 
agreement  with  the  total  ionization  cross-sections.  '  The  same  basic 

function  proposed  by  Bretagne,  et  ai.^  for  argon-based  excited  states  is 
utilized: 


where  a(ep,e)  is  the  differential  ionization  cross-section  (in  A)  for  a 

primary  and  secondary  electron  of  energy  and  c,  respectively;  is  the 

ionization  energy;  and  Eg  is  the  ionization  energy  of  hydrogen.  K,  a,  p,  and 

7  are  fitting  parameters.  It  is  found  that  the  values  a='f=3,  p=0.9,  and 

K=2.49  give  fairly  good  agreement  between  the  integrated  total  ionization 

cross-section  and  the  reported  values^^  for  ionizations  of  Ne*(3s),  Ne**(3p), 

Xe*(6s),  and  Xe*'''(6p).  These  same  parameters  are  applied  to  Ne^,  and 

* 

NeXe  .  K=0.5  with  the  same  a,  P,  and  7  given  previously  yields  reasonable 

33 

agreement  for  HCl  ionization. 

The  calculated  results  of  the  characteristic  time  constant  [sec  for 

3 

various  species  are  listed  in  Table  V  for  an  excitation  rate  of  1.0  MW/cm  . 

-1  -3 

The  actual  formation  rates  [sec  cm  ]  are  obtainable  by  multiplying  the 
densities  of  the  parental  species  by  these  time  constants.  These  formation 
rates  are  proportional  to  the  e-beam  deposition  rate.  Therefore,  once  the 
high-energy  distribution  and  the  formation  rates  of  species  has  been 
calculated,  these  values  can  be  used  during  the  e-beam  excitation  by 


proportionally  adjusting  the  rates  to  the  instantaneous  e-beam  energy 
deposition  rate.  It  should  be  emphasized  that  during  these  calculations,  the 
W-values  and  excited  state  yields  are  not  the  primary  parameters  needed  to 
predict  the  secondary  electron  energy  distributions. 

19 

As  shown  elsewhere,  calculation  of  the  low-energy  electron 
distribution  requires  a  time-dependent  Boltzmann  equation,  which  includes 
electron-electron  collisions,  elastic  collisions,  and  all  low-energy 

reactions.  This  is  solved  in  a  manner  similar  to  that  reported  by  Bretagne, 
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et  atl.  for  pure  argon  gas  and  in  conjunction  with  a  detailed  heavy  particle 
kinetics  code.^ 

The  preceding  detailed  treatment  is  paraunount  when  modeling  the 
secondary  electrons  in  e-beam  produced  laser  gas  mixtures.  Not  only  does  the 
non-Maxwellian  distribution  of  the  high-energy  electrons  affect  the 
predictions,  but  also  the  effect  of  local  electron  loss  in  the  energy 
spectrum,  caused  by  dissociative  attachment  reaction  with  halogens,  can  cause 
significant  differences  in  the  transient  concentrations  of  electrons  and 
excited  species,  halogen  consumption,  and  so  on,  from  those  calculated  by 
models  that  assume  a  Maxwellian  distribution.  This  is  especially  true  when 
the  e-beam  deposition  rate  is  low,  and  the  electrons,  now  at  relatively  low 
density,  cannot  effectively  create  a  Maxwellian  distribution. 


IV.  CONCLUSIONS 


The  distribution  functions  of  high-energy  electrons  in  e-beam  excited 
rare  gas  mixtures  are  calculated  for  the  first  time.  The  steady-state 
distributions  for  energies  higher  than  the  first  excitation  threshold  of  the 
major  rare  gas  are  quickly  established  and  a  relatively  large  number  of 
electrons  are  populated  in  the  high-energy  tail  as  compared  with  a  Maxwellian 
distribution.  Competition  between  the  various  elastic  collisions  of  the 
different  rare  gases  in  the  mixture  result  in  a  dependence  of  the  W-values  and 
excited  state  yields  on  the  mixture  composition.  Differences  in  the 
predictions  for  the  W-values  and  excited  state  yields  of  the  minor  rare  gases 
can  be  particularly  large  when  the  values  for  a  pure  rare  gas  are  used  rather 
than  those  calculated  for  a  mixture. 
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TABLE  II.  Electronic  classifications  and  their  energy  levels  used  in  the  electron  distribution 
calculations . 


TABLE  III.  Calculated  W-values  and  yields  of  ions  and  metastables  for  neon 
and  xenon  based  gas  mixtures.  The  efficiencies  are  defined  by  the  ratio  of 
the  W-value  to  the  ionization  or  excitation  threshold  energy.  The  number  in 
parenthesis  corresponds  to  the  power  of  10. 


SPECIES 


W(eV) 


YIELD 


EFFICIENCY («) 


3000  Torr 

Ne"" 

3.8(1) 

1 

Ne* 

1.7(2) 

2.4(-l) 

Ne*" 

1.3(2) 

2.9(-l) 

Ne*** 

1.5(3) 

2. 6 (-2) 

Ne**** 

9.0(3) 

4.0 (-3) 

760  Torr 

Xe* 

2.2(1) 

1 

Xe* 

1.3(2) 

1.8(-1) 

Xe** 

1.2(2) 

1.8(-1) 

Xe*** 

9.1(2) 

2.0(-2) 

1  =  2985/15  Torr 

Ne"" 

3.9(1) 

1 

Ne* 

1.7(2) 

2.2(-l) 

Ne** 

1.4(2) 

2.8(-l) 

Ne*** 

1.5(3) 

2. 6 (-2) 

Ne**** 

9.2(3) 

4.0(-3) 

Xe"" 

3.5(2) 

l.l(-l)  [1] 

Xe* 

6.4(2) 

6.0(-2)  [5.5(-l)] 

Xe** 

4.6(2) 

8.5(-2)  [7.7(-l)] 

Xe*** 

6.5(3) 

6.0(-2)  [5.0(-2)] 

TABLE  lY.  Calculated  W-values  and  yields  of  ions  and  metastables  for  argon 
and  krypton  based  gas  mixtures.  The  efficiencies  are  defined  by  the  ratio  of 
the  W-value  to  the  ionization  or  excitation  threshold  energy.  The  number  in 
parenthesis  corresponds  to  the  power  of  10. 


SPECIES 

W(eV) 

YIELD 

EFFICIENCY (R) 

Ar  =  1020  Torr 

Ar"" 

2.5(1) 

1 

6.3(1) 

Ar* 

1.8(2) 

1.3(-1) 

6.1 

Ar** 

2.7(2) 

9. 5 (-2) 

4.9 

At*** 

2.9(2) 

8. 8 (-2) 

4.9 

Ar**** 

6.6(2) 

3. 8 (-2) 

2.3 

Ar 

2.7(2) 

9. 5 (-2) 

5.8 

Kr  =  760  Torr 

Kr"" 

2.1(1) 

1 

4.8(1) 

Kr* 

1.1(2) 

2.0(-l) 

1.1(1) 

Kr** 

2.0(2) 

l.l(-l) 

6.1 

Kr*** 

1.5(3) 

1.4(-2) 

8.8(-l) 

^  ♦♦♦♦ 

Kr 

4.7(3) 

4.0(-3) 

3.0(-l) 

v.v.  .r. .r. . . 
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TABLE  IV.  (cont.) 


1020/114  Torr 

At"" 

2.9(1) 

1 

At* 

2.2(2) 

1.3(-1) 

Ar** 

3.0(2) 

9. 6 (-2) 

Ar*** 

3.2(2) 

8. 9 (-2) 

At**** 

7.4(2) 

3.9 (-2) 

At 

3.0(2) 

9. 5 (-2) 

Kr* 

1.6(2) 

1.8(-1)  [1] 

Kr* 

6.7(2) 

4.3(-2)  [2.4(-l)] 

Kr** 

1.3(3) 

2.2(-2)  [1.2(-1)] 

Kr*** 

1.2(4) 

2.0(-3)  [1.4(-2)] 

Kr**** 

3.8(4) 

8.0(-4)  [4.0(-3)] 

TABLE  V.  Calculated  formation  rates  of  various  species  created  by  the  e-beam 
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source  (e>16.7  eV)  and  by  10  eV  primary  electrons  for  a  mixture  of  Ne/Xe/HCl 
=  2980/15/4.8  Torr  at  an  excitation  rate  of  1  MW/cm^.  The  number  in 
parenthesis  corresponds  to  the  power  of  10. 


REACTION 


Xe  Xe 


TIME  CONSTANT  (sec"^) 


FIGUKB  CAPTIONS 


FIG.  1.  Distributions  of  ionization  and  excitation  cross-sections.  Dashed 
curves  are  for  krypton-based  ionization  and  excitations.  Solid 
curves  are  argon-based  ionization  and  excitations. 

FIG.  2.  Distributions  of  ionization  and  excitation  cross-sections.  Dashed 

curves  are  xenon-based  ionization  and  excitations.  Solid  curves  are 
neon-based  ionization  and  excitations. 

FIG.  3.  Temporal  evolution  of  the  high-energy  electron  distribution 
calculated  for  3000  Torr  pure  neon  at  an  excitation  rate  of 
0.82  MW/cm^. 

FIG.  4.  High  energy  tail  of  the  steady-state  electron  distribution  calculated 
for  3000  Torr  pure  neon.  Dashed  curve  is  a  Maxwellian  distribution 
at  a  characteristic  temperature  kT^  =  4.6  eV. 

FIG.  5.  Temporal  evolution  of  the  high-energy  electron  distribution 
calculated  for  760  Torr  pure  xenon  at  an  excitation  rate  of 
0.41  MW/cm^. 

FIG.  6.  Temporal  evolution  of  the  high-energy  electron  distribution 

calculated  for  a  rare  gas  mixture  of  Ne/Xe  =  2985/15  Torr  at  an 
excitation  rate  of  1  MW/cm  .  Dotted  curve  is  the  distribution 


obtained  at  t  =  100  ps  for  a  Ne/Xe  =  2850/150  Torr  mixture  at  an 

3 

excitation  rate  of  1.1  MW/cm  . 

FIG.  7.  High-energy  electron  distribution  at  49  ps  (approximately  steady 

state)  calculated  for  1020  Torr  pure  argon  and  760  Torr  pure  krypton 

3 

at  excitation  rates  of  0.27  and  0.34  HW/cm  ,  respectively. 

FIG.  8.  Steady  state  high-energy  electron  distribution  calculated  for  rare 
gas  mixtures  of  Ar/Kr  =  1020/51,  1020/114,  and  1020/510  Torr,  at  an 
excitation  rate  of  0.52,  0.57,  and  0.88  HW/cm^,  respectively. 

FIG.  9.  Plots  of  W-values  normalized  by  the  partitioned  deposition  energy  for 
the  argon  ion,  and  yields  of  excited  states  of  argon  relative  to  the 
argon  ion  formation  in  Ar/Kr  mixtures  as  a  function  of  argon 
concentration. 

FIG.  10.  Plots  of  W-values  normalized  by  the  partitioned  deposition  energy  for 
the  krypton  ion,  and  yields  of  excited  states  of  krypton  relative  to 
the  krypton  ion  formation  in  Ar/Kr  mixtures  as  a  function  of  krypton 


concentration . 
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Abstract 

The  electron  energy  distributions  for  energies  lower  than  17  eV  and 
their  time  dependent  evolution  are  calculated  for  electron-beam  (e-beam) 
excited  Ne/Xe/HCl  mixtures.  A  time  dependent  Boltzmann  equation  including  all 
interactions  between  electrons  and  ground  or  excited  state  species  is  solved 
together  with  a  detailed  full  kinetics  model  for  XeCl  lasers.  Effects  of 
electron-electron  collision,  HCl  concentration,  and  excitation  rate  on  the 
steady-state  electron  energy  distribution  are  examined. 
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High  pressure  gases  excited  by  intense  electron-beams  (e-beams)  have 
been  used  to  produce  high  densities  of  excited  or  ionized  species  leading  to 
high  poirer  laser  emission,  such  as  in  rare  gas  dimer  and  rare  gas  halide 
lasers.  In  such  a  plasma,  most  of  the  ions  and  excited  species  are  formed 
either  directly  by  the  primary  or  the  high  energy  secondary  electrons.  At 
high  electron  densities,  the  effects  of  electron  impact  excitation  and 
quenching  can  be  important  even  for  relatively  low  energy  electrons  (lower 
than  the  ionization  energies  for  the  rare  gas  constituents) . 


The  mechanism  of  electric  discharge  excitation  employed  in  discharge 

pumped  lasers  is  fairly  well  understood,  and  theories  for  the  electron  energy 
1*3 

distribution  ,  which  relates  the  discharge  parameters  to  the  excitation 

rates,  are  available.  However,  a  similar  detailed  electron  energy 

distribution  theory  has  been  lacking  for  e-beam  excited  lasers.  Models 

published  to  date  typically  assume  that  electrons  are  in  a  Maxwellian  energy 

distribution  with  either  a  constant  or  a  time  dependent  average  electron 

4*8 

temperature;  this  is  particularly  true  for  rare  gas  halide  lasers^ 


Eliott  aind  Greene®,  and  Bretagne  et  al.^®’^^  have  reported  electron 

energy  distribution  functions  for  pure  xenon  and  argon  gases,  respectively. 

Their  calculated  distribution  functions,  using  a  detailed  Boltzmann  equation 

12 

solver  that  includes  the  Rockwood  formalism  for  electron-electron  and 
elastic  collisions,  show  significant  departures  from  Meucwellian 
distributions.  Similar  calculations  performed  for  mixtures  have  been 

reported  by  Rozenberg  et  al.  In  their  mixtures,  the  dissociative  electron 


M'S'* 
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attachment  reaction  of  F2  plays  an  important  role  in  influencing  the 
distribution  of  electrons  lower  than  the  excitation  threshold  energy  of  argon 
gas  (11.8  eV) .  Fluorine  tends  to  increase  the  average  electron  energy  as  a 
result  of  removal  of  low  energy  electrons  by  attachment  reactions. 

In  rare  gas  mixtures  containing  a  trace  of  HCl,  the  relatively  large 
vibrational  excitation  cross-sections  of  ECl  and  the  large  dissociative 
attachment  rate  constants  of  vibrationally  excited  HCl  can  significantly 
affect  on  the  electron  distribution  function.  A  detailed  calculation  of  the 
electron  distribution  function  is  paramount  in  order  to  obtain  the  correct 
rates  for  reactions  between  electrons  and  ground  or  excited  state  species. 
These  rates  are  necessary  for  predicting  the  densities  of  various  transient 
species  in  laser  gas  plasmas,  such  as  the  secondary  electrons,  the  halogen 
molecules,  and  the  excited  rare  gas  atoms. 

In  this  letter,  calculated  results  of  the  low  energy  electron 
distributions  in  e-beam  excited  Ne/Xe/HCl  laser  mixtures  are  described.  The 
numerical  technique,  which  employs  a  Boltzmann  equation  that  includes  all  the 
collisional  processes  in  which  the  electrons  are  involved,  is  the  same  as  the 
one  applied  by  Bretagne  et  al.^^  to  pure  argon  gas. 

The  high  energy  electron  distribution  function  calculated  for  a  Ne/Xe 
mixture  using  a  reduced  Boltzmann  equation,  which  was  reported  in  our  previous 
paper^^,  is  used  as  the  driving  term  for  the  full  Boltzmann  equation  of  the 
low  energy  electron  distribution.  Since  electrons  higher  than  the  excitation 


threshold  of  neon  (16.7  eV)  establish  a  steady-state  distribution  in  a  very 
short  time  (less  than  100  ps) ,  the  high  energy  distribution  function  can  be 
treated  as  static  in  overall  shape,  but  changing  in  magnitude  proportional  to 
the  temporal  change  of  the  e-beam  excitation  rate.  It  should  be  emphasized 
that  W-values  (defined  as  the  energy  required  to  form  an  electron-ion  pair) 
for  rare  gas  ion  formation,  and  yields  of  rare  gas  excited  species  relative  to 
the  ion  formation,  which  are  typically  used  in  models  for  e-beam  excited  gas 
lasers,  are  not  used  in  our  calculations.  In  our  previous  paper  ,  it  was 
concluded  that  the  dependences  of  these  values  is  sensitive  to  the  mixture 
ratio  of  the  rare  gases  (e.g.  Ar/Kr  or  Ne/Xe)  and  that  the  definitions  of 
these  parameters  can  be  ambiguous,  particularly  when  applied  to  the  minor  rare 
gas  in  these  mixtures.  Therefore,  by  using  a  high-energy  electron 
distribution  function  rather  than  these  macroscopic  values,  the  production  of 
ions  and  excited  species  can  be  more  accurately  predicted. 


As  reported  by  Bretagne  et  al.^^,  the  evolution  of  the  distribution 
function  f(e,t),  for  energy  e  (lower  than  the  excitation  threshold  energy  of 
neon)  is  expressed  as  follows: 
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where  the  first  term  on  the  right-hand  side  corresponds  to  electron-electron 
collisions;  the  second  to  elastic  electron-ground  state  species  collisions; 
the  third  to  inelastic  collisions  including  ionizations,  excitations, 
superelastic  de-excitations,  attachments,  recombinations,  and  Penning  effects; 
and  the  last  term  to  the  production  of  secondary  electrons  with  energy  e  by 


the  primary  electrons  Cp  of  the  e-beam. 


Electron-electron  and  electron-neutral  collisions  are  modeled  using  the 

Rockwood^^,  and  Elliott  and  Greene®  formulas  as  modified  by  Bretagne  et  al.^^ 

to  adapt  to  non-equal  segments  for  the  energy  intervals.  Ionization  terms  are 

calculated  using  differential  ionization  cross-sections  reported  by  Green  and 
15 

Sawada  for  ground  state  rare  gases  and  using  empirical  formulas  obtained  for 

HCl  and  excited  rare  gas  atoms^^.  The  set  of  cross-sections  given  in  our 

previous  paper^^  are  used  to  calculate  the  excitations  of  the  rare  gases. 

16 

Cross-sections  reported  by  Davis^  are  used  for  HCl  ionizations  and 

17 

excitations,  and  those  reported  by  Domcke  and  kundel  are  used  for 
dissociative  attachment  reactions  with  ground  and  vibratioaally  excited  HCl. 


During  the  calculation,  the  initial  secondary  electron  production  rates 

for  each  of  the  energy  segments  lower  than  16.7  eV  are  calculated  from  the  e- 

beam  electron  energy  and  the  static  high-energy  electron  distributions.  These 

calculations  are  performed  only  once  at  the  beginning.  Afterwards,  the  rates 

are  adjusted  in  direct  proportion  to  the  time-varying  e-beam  excitation  rate 

and  the  density  of  parental  species  during  e-beam  excitation.  The 

differential  equations  [i.e.  Eq.  (1)]  obtained  for  the  electron  populations  in 

each  of  the  energy  segments  lower  than  16.7  eV  are  simultaneously  solved  using 

a  fourth-degree  Runge  Kutta  integration  method.  An  important  difference 

9 

between  our  calculations  and  the  ones  done  by  Elliott  and  Greene  ,  or  Bretagne 

11  6 
et  al .  is  that  a  detailed  heavy  particle  kinetics  model”  is  simultaneously 

solved  together  with  the  Boltzmann  equation.  This  permits  a  more  accurate 


accounting  of  the  transient  ions  and  excited  species  densities  that  are 
involved  in  the  Boltzmann  equation. 

Running  both  models  in  a  true  continuous  parallel  mode  is 

computationally  prohibitive  on  the  computer  system  (VAX  11/780)  used  for  this 

work.  Therefore,  the  following  algorithm  is  used  to  conserve  on  computation 

time;  1)  During  the  kinetics  calculation,  the  electron  energy  distribution  is 

frozen,  thereby  making  the  reaction  rate  constants  involving  electrons 

constant.  2)  When  the  electron  density  changes  by  5%,  the  Boltzmann 

calculation  is  called.  During  the  Boltzmann  calculation,  the  kinetics 

calculation  is  stopped,  and  all  the  densities  of  species  are  frozen.  The 

Boltzmann  equation  calculation  is  continued  until  the  energy  distribution 

_2 

function  reaches  steady-state  using  a  convergence  factor  of  10  (defined  as 
the  ratio  of  the  previous  calculated  value  to  the  current  calculated  value) . 

3)  The  kinetics  calculation  with  new  rate  constants  for  the  electron  reactions 
restarts  and  the  process  repeats  itself . 

Figure  1  shows  the  time  evolution  of  the  low-energy  electron 
distribution  function  calculated  for  a  XeCl  laser  mixture  of  Ne/Xe/HCl  = 

18 

2980/15/4.8  Torr.  An  e-beam  current  waveform  from  the  Tahoma  laser  device 

3 

at  an  excitation  rate  of  250  kW/cm  and  a  pulse  width  of  500  ns  full  width  at 
half  maximum  (FWHM)  is  used.  Early  in  time,  the  distribution  function  is 
truncated  near  the  first  excitation  threshold  of  xenon  (8.4  eV) ,  which  was 
also  observed  during  the  high  energy  electron  distribution  calculation^^,  and 
indicates  that  ionization  and  excitation  processes  of  the  ground  state  rare 


gas  quickly  remove  electrons  higher  than  the  rare  gas  energy  threshold. 
Structure  seen  on  the  distribution  curve  at  energies  less  than  10  eV  are 
caused  by  vibrational  excitation  of  HCl. 

At  later  times,  as  the  electron  number  density  increases,  the  electron- 
electron  collisions  tend  to  redistribute  the  low  energy  electrons  towards  a 
Maxwellian  distribution  at  an  average  electron  energy  of  i^l.S  eV  during  the 
time  interval  50-600  ns.  Superelastic  de-excitation  reactions  also  become 
important  as  the  densities  of  excited  species  increases  and  contribute  towards 
a  smoothing  out  of  the  distribution  function  caused  by  heating  of  the  lower 
energy  electrons.  At  this  HCl  concentration,  dissociative  attachment  and 
vibrational  excitation  reactions  are  slower  than  the  redistribution 
processes. 

There  are  two  important  points  to  be  noted.  First,  the  induction  time 
to  establish  a  steady  state  distribution  is  ^50  ns  for  the  condition  studied. 
This  time  depends  on  the  risetime  of  the  current  pulse  and  the  excitation 
rate.  Second,  the  part  of  the  distribution  function  higher  than  the  first 
excitation  threshold  of  xenon  shows  significant  departure  from  a  Maxwellian 
distribution,  even  with  the  redistribution  effects  caused  by  electron-electron 
processes  and  superelastic  de-excitation  reactions.  Therefore,  assuming  a 
Maxwellian  distribution  for  the  secondary  electrons  can  cause  a  significant 
error  during  the  calculation  of  the  rate  constant  of  electron  reactions,  which 
will  result  in  errors  in  the  transient  species  densities. 


In  order  to  confirm  the  effect  of  electron-electron  collisions,  Pig.  2 
shows  the  results  of  the  steady-state  electron  distribution  as  a  function  of 
the  e-beam  excitation  rate  calculated  including  the  electron-electron 
collision  process  [Fig.  2(a)]  and  without  it  [Fig.  2(b)].  The  gas  mixture  is 
the  same  as  in  Fig.  1.  It  is  clear  that  even  at  a  low  excitation  rate  of  2 
kW/cm  ,  the  electron-electron  collision  process  is  very  important.  The 
structures  seen  on  the  curves  in  Fig.  2(b)  exactly  correspond  to  the 
vibrational  excitation  of  HCl.  From  Fig.  2(a),  however,  it  is  observed  that 
the  distribution  functions  at  low  excitation  rates  do  slightly  depart  from 
Maxwellian  distributions,  although  they  are  smoothed  out.  At  this  particular 
HCl  concentration,  the  required  excitation  rate  to  establish  a  Maxwellian 
distribution  in  the  low  energy  region  is  approximately  50  kW/cm  . 

Figure  3  shows  steady-state  distribution  functions  calculated  for 

various  HCl  concentrations.  The  e-beam  excitation  rate  is  kept  constant  at 
3 

200  kW/cm  .  Since  the  distribution  function  in  the  low  energy  region  is 
established  as  a  result  of  the  balance  between  the  electron-electron 
collisions  and  vibrational  excitation  of  HCl,  HCl  excitations  become  important 
at  high  HCl  concentrations  and  quench  electrons  higher  than  the  vibrational 
excitation  threshold. 

Under  the  more  specialized  situation  of  relatively  high  excitation  and 
relatively  low  initial  HCl  concentration,  an  alternative  simpler  model  that 
assumes  a  quasi-Maxwellian  distribution  is  possible.  Under  these  conditions, 
electron-electron  collisions  can  effectively  establish  a  Maxwellian 


distribution  for  the  low  energy  electrons.  In  this  simplified  approach,  the 
distribution  function  of  the  low  energy  region  would  be  treated  as  Maxwellian 
whose  characteristic  temperature  is  defined  by  solving  an  energy  balance 
equation  for  the  low  energy  electrons  during  excitation.  It  should  be 
emphasized  that  even  for  this  specialized  case,  only  the  low  energy  electron 
distribution  can  be  approximated  as  Maxwellian.  The  high  energy  regions  would 
still  be  treated  as  a  static  non-Maxwellian  distribution  as  calculated  from 
the  reduced  Boltzmann  equation^^.  Hence,  this  simplified  approach  still 
differs  from  the  models  which  assume  that  the  entire  distribution  is 
Maxwellian. 

The  modeling  process  described  in  this  letter,  in  conjunction  with  that 
reported  previously  for  the  high  energy  electron  distribution,  is  a  more 
accurate  alternative  to  the  usual  modeling  methods  that  use  stopping  power,  W- 
values  of  pure  rare  gases,  excitation  yields,  and  a  Maxwellian  assumption  for 
the  entire  electron  energy  distribution.  This  new  approach  should  improve  the 
performance  of  the  model  predictions.  It  is  also  applicable  to  other  kinds  of 
lasers  excited  by  e-beams. 

This  work  was  supported  by  the  Office  of  Naval  Research,  Contract 


No.  N0014-85-C-0843. 
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FIGURE  CAPTIONS 


Fig.  1.  Temporal  evolution  of  the  low-energy  electron  distribution  function 

calculated  for  a  mixture  of  Ne/Xe/HCl  =  2980/15/4.8  Torr  by  solving 

the  time-dependent  Boltzmann  equation  together  with  a  full  kinetics 

3 

model.  The  excitation  rate  is  f3250  kW/cm  and  the  e-beam  pulse  width 
is  500  ns  (FWHM) .  The  distribution  functions  between  100  ns  and 
560  ns  sure  approximately  same  as  the  one  at  100  ns . 

Fig.  2.  Steady  state  low-energy  electron  energy  distribution  functions  as  a 
function  of  e-beam  excitation  rate  for  a  mixture  of  Ne/Xe/HCl  = 
2980/15/4.8  Torr.  (a)  Electron-electron  collisions  included.  (b) 
Electron-electron  collisions  excluded. 

Fig.  3.  Steady  state  low-energy  electron  energy  distribution  functions  as  a 

function  of  HCl  concentration.  The  mixture  contains  0.5%  Xe  with  the 

3 

balance  Ne  to  3000  Torr.  The  e-beam  excitation  rate  is  200  kW/cm  . 
Electron-electron  collisions  are  included. 
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